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SYNOPSIS 

i'he first part of the present study which deals with 
soLie q_uaiiti ba Give aspects of flotation kinetics has been sub— 
dividea asJ (l) a phenomenological model for flotation kinetics; 

(2) an algorithm for estimation of flotation rate distribution; 

(3) the effect of process variables on the flotation rate dis~ 
cribucions; and (4) the effect of process variables on the flota- 
tion rats distribution in a continuous operation. The second par 
deals ?/ith (5) the optimal synthesis and design of a flotation 
circuit. 

^ • 4., jjlgJlog§Ilol..Pbi cal Model for Plotat io n Ki net ics: 

A phenomenological model for flotation kinetics has 
been derived ?bich provides an explicit relationship between 
the apparent rate constant and the specific rate constant for 
particle-bubble attachment subprocess. In first instance a 
general forniulation of the kinetics hci s been constructed in terms 
of a set of non-linear differential equations. Mathematical 
conditions for the free flotation caseumder vbich flotation 
process can be described by a simple first order rate equation 
are stipulated, A simplified method for analysis of kinetics 
in terms of grade of the concentrate has been proposed in which 
flotation kinetics is concurrently analyzed for total mass 
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recovery anct mass of valuable mineral. 


2 . 


-Bate Pi s t r i but lo n B v N n m e r i ca 1 


Ihe 

function in 


problem of estimation of the rate distribution 
appaxent rare constant in first order flotation 


rate equation is equivalent to the inversion of the Laplace 
transform. A reliable method for estimating the flotation 


rate constant distribution iia cumulative mode by the numerical 
inversion of the Laplace transform has been developed. 'Ihe- pro- 
posed algorithm has been tested for S5cnthetic data generated 
from different distribution fmictions. The effects of errors 
in the data, number of data points, initial guess in the optimi- 
zation step and non-floatables on the accuracy of the proposed 
algorithm have also been studied. 

3 • L he Sffec t oi Pro cess Ya rlables on the flotatio n Rate 
■ List r i but i 0 n s ; 

The computed cumulative distributions in apparent rate 
const, '.nt 5 as a whole, have been examined for appropriate rela- 
tionships between cumulative distribution and aeration rate, 
particle size and pulp density. The distribution function is 
found to remain unalte.red, but the argxmaent of the function, 
the apparent rate constant, is rescaled with change in the rate 
of aeration and particle size. The scaling feictor increases / 
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rapidlj- with the rate of aeratior aM 
max iiiTujn for interrneciiate particle size 


exhibits a characteristic 
The calculated distri- 


butions are independent of initial pulp density when dilute, but 
show marked divergence for concentrated pulp. A partial proof 
for ohe validity of the model proposed in part 1 has been given. 


4 . 


ge ed Rat on,. tjpe^ljTi^^^ 
a Co ntinuou s Cel l : 


An algorithm for estimating the floteition rate distribution 
in a continuously operated cell has been developed and using the 


tabulated data given in the literature, the effects of particle 
size and feed rate on the flotation rate distribution have been 
exaiDined. The distributions are found to be independent of 
slurry flow rate and the effect of particle size is similar to 
that in the semi— batci] process. The estimated rate distributions 
have been employed to predict the kinetics of cleaner cell whose 
feed comprises of the rougher concentrate. Similar results have 
bean obtained for the gangue component. 


5 • and Desi /ni of Jlo^t^apn,, : 

A preliminarj- investiga-cion of the feasibj.lity of synthe- 
sizing and designing optimal and suboptima 1 flotation circuits 
by optimiza'cion of a number of appropriate objective functions 
of recovery ana grade and profit has been carried out. In this 


approach, the optimal structure aid. design parameters are 
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e^^trac 1:6(3. simultaneously from a generalizeci circuit by direct 
search. Pirst, a flotation circuit with two cells and feed 
comprising of two species of kno\TO overall rate constant is 
synthesized and designed and the results are compared with the 
real optiraum circuit for this case, obtained by using the trial 
and error technique. The approach is then extended to a flota- 
tion circuit with four cells and for batteries of cells for a 
feed comprising of valuable, gangue and middlings. 



CHAPTER I 


lETRODUCTION 

Trotn flotation of particulate solids, suspended in 
v.-ater, by means of air buobles is perhaps the Diost important 
technique available at present for separation of valuable 
minerals from associated gangue impurities. The commercial 
exploitation of froth flotation can be traced back to first 
decade of the twentieth century.'’ Probably no other process 
in the history of mineral beneficiation has been responsible 
for such increased mineral production as froth flotation. At 
present, many low grade and finely grained complex ores, which 
were previously considered worthless, are mined only because 
they can be economically enriched to concentrates by the 
flotation process. 

In viev7 of the importance of the flotation process, 
attempts are continuously being made to obtain better under- 
standing of the underlying kinetic mechanisms and physical 
phenomena which may lead to: 

(1) Development of the process strategy for ores which, 
at present, cannot be beneficiated by flotation. 

(2) Modification of the existing processes, so that as the 
mineral rich ores are depleted these processes may be adopted 
to lo'w quality ores. 
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(3) Iiaprovement in the efficiency of the existing processes 
for taprovement in the quality of the product, the plant 
throughput and the operating cost. 


Once 


improvement 


tne basic process has been evolved, any f-urther 
in its efficiency is in general possible only by 


optimization and automatic control of the process. Ihis, 
however, requires a reasonably accurate mathematical kinet 
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model ^^lich should be able io predict the performance of the 
unit opexation under varying operating conditions. Because 
of the large tonnage of ores treated, even a marginal improve- 
ment in the process efficiency can have significant economical 
JJnpact. Thus, the importance of a quantitative understanding 
of flotation as a rate process is obvious. 

At present, it is not possible to derive the flotation 
kineuic model from first principles, although some attempts 
have been made recently in this direction. 1 feasible and 
useful approach is to develop a kinetic model which is in 
agreement ?/itb overall flotation kinetics and then to study 
systematically the effects of Vu,ricus process variables on 
the specific rate parameters appearing in the model. However, 
the nuiuDer o± variables which have an appreciable effect on 
the flota,-cion behaviour is very large. An exhaustive list of 
the variables has been given by Sutherland and Wark^, 

but as suggested by Arbiter and Harris^, this problem can be 
somewhat simplified by grouping together the varia.bles 
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pertaining to (l) ore and mineral properties; (2) nature of 
reagi^nt tieaument; and (3) aggregate of flotation machine 
characteristics. Since the flotation kinetics depends on 
nonlinear interactions oi the variable Sj this classification 
is^ strictly speaking, not mutually exclusive. Group 1 includes 
Particle sise^ shape anu density, and mineral composition and 
distribution etc. whereas variab-les related to adsorption or 
reaccion oi chemical reagents, e. g. collector, frother and 
modifier etc. at the siorface of mineral particle and air bubbles 
cons'citute group 2. In group 3, the most important variables 
are the number of bubbles and their size distribution, and the 
complex hydrodynamic environment prevailing in the cell. 


Although a large number of studies on flotation kinetics 
have been reported, a majority of these ‘are experimental and 
empirical in nature and a comprehensive model, which can 
interpret the existing kinetic data satisfactorily and further 
predict the performance of the process under varying operating 
conditions, is ye g to be evolved. In the present study an 
attemvpt has been made to study some aspects of flotation kine- 
oico leading to development of a more 'Tieanlngful kinetic model 
based oi> the dist3?ibuted particle-bubble attachment rate 
constant. The resulting model has been reduced to distributed 
apparent flotation kinetics, v/hich has been highly successful 
for the analysis of the existing kinetic data, under free 
flotation constraints. A method for numerical inversion of the 
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Laplace transform has been proposed which is suitable for 
e^ Li..u<., ti^,, of the apparent flotation rate distribution. The 
effects of some of the process variables on flotation kinetics 
ueem otjdied and some of the results are interpreted on 
the basis of the new kinetic model. 1‘inally, optimal synthesis 
and deoj-gn of t^ic industrial flotation circuits has been 

carried out. These objectives are discussed in more detail0 
in Chapter 3. 



CI-IAPTER 2 


I'lTERATURS REVIEW 

There exists considerable literature on flotation 
processes, for convenience, it can be divided into five 
broad catetiories: literature pertaining to i) physical 
cheiiiistxy ox flotation; 2) kinetics of the process; 3) design 
of flotation cell and various other technological aspects; 

4) optimization of the flotation circuit and process vari-- 
ables; and 5) automatic control of the process. It is beyond 
the scope of this study to review all these areas; only, 
i berature perociining to secono. and fourth categories is 

discussed under the titles el llotatiorj kluetios and optUol- 
zation. 

A. I'lOTATIOh ICIKETICS: 

'The mabheiaatical description of flotation as a rate 
process vjas first given by Gracia-Zuniga^ ?/ho, on the basis 
of ejbperimeiitai dota in a semi-batch process, reported that 

the material remaining in the pulp decays exponentially with 
time t , i.e. 


h^3(t) 


exp [-Ea t] 


( 2 . 1 ) 


, I ' ■ 

wheru. Mg(t) is fraction of material remaining in suspension 
at time t and Ka is apparent or overall flotation rate constani 
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This empirical eq-aation means that a plot' of In (Bil’(t)) versus 
t should he a straight line with slope equal to minus Ka. Since 
then Several kinetic equations and models have been proposed. 

In this revie?/, these models have been discussed in more or less 
the same order in which they were evolved and this is reflected 
in the degree ox sophistication in models which tends to increase 
as ?/e go from one model to the next. 


s Based on Ch emical Ki netic Analogy* 

Starting from the premise that the rate of flotation 
can be expressed as a first order chemical reaction, Beloglazov^ 
proposed a differential equation for a semi-batch flotation 
process, similar to that of G-racia-Zuniga^, >Yhich can be written 
as 


dT' 


Ka C(t) 


( 2 . 2 ) 


where C is the concentration of solid particulate species in 
pulp, ihe first order kinei-ic form of the rate expression in 
equaoion (2.r) prooably led drunder and Kadur*^, and later, 
Plaxsin, ivlassen ana Berger^ to adopt the law of xiass action 
in chemical kinetics for their kinetic equations for flotation 
process. Superficially, the chemical kinetic analogy seems 
quite reasonable if one considers the formation of bubble 
particles aggregate as the rate controlling step, similar to 
the chemical reactions governed by interaction between atoms; 
molecules or ions. According to this view point, a general 
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equation for the kinetics of the semi-hatch process can he 
postulated as^'. 



where, as defined e..rlier, G refers to the concentration of 
flotation material, is the concentration of a particular 
constituent i in the pulp, k' is specific flotation rate 
constant, n^ is the order of reaction v/ith respect to i-th 
species and n^ is the number of species in the pulp taking 
part in aggregation step. If one assumes that there are only 
two constituents in the pulp, namely, solid particles and air 
bubbles, and if the air supply is constant, merger of all 
constant tertas ?\dth k reduces equation (2,3) to 


d_C^(t) 

dt"" 



( 2 . 4 ) 


where n is the overall order of reaction. Schuhmann^ pointed 
out that for n = 1, the rate constant in equation (2.4) has the 
same dimensional significance as the specific reaction rate 
constant for a first order chemical reaction. Mika"^*^ has, ' 
however 5 ? ,correc1;ly'. poiiited out- that 'the analogy .balweeri ■ 
flotation and chemical reaction rate equations arises from the 
use of bhe terms, such as rate constant and order of reaction, 
and these parameters, as such, have no physical significance 
in flotation phenomena except when n = 
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There has heen some controversy in the literature 
concerning the value of n in equ 3 ,tion (2.4). Gracia-Zuniga^ 
Beloglazov®, Grunder and ICadur^, Schuhraann^, Sutherland'^'’, 
iiiOrris nrovm and Sraith'’"', Jowett and Safvi"’^ 


Volkova 


12 


Modi and i'uerstenau/'^, Imaizumi and Inoue 
Young postulated n to ht 
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unity, whereas Irbiter^-^, Bennett 


and Chi and 
19 


Chapaan and Dell 2°, Hum21, Mitrofanovas. 23 and 


Salman2'^ proposed a value of tv». Bogdanov and co™rkers25.26 
Horst ^ ■ , deBruyn and Modi ^9^ and Volin and Svami^O assigned 
various other values. Plaksin et al^ asserted that it would 
he oncorrect to assume a constant value for n during the entire 
period of flotation because conditions and forms of the minera- 
lized hubbies undergo change during the process. Zrokhin^'’ and 
Kanovskii^^ divided the whole process into several stages and 
described each stage by a different equation. Examination of 


systems that exhibit non-first order kinetics led Pogorely^^, 
and Tomlinson and Fleming^'^ to introduce the concept of free 
and inhibited flotation. In the case of free flotation, 
whicn requires dilute pulp and high aeration rate, the bubbles 
o.re only sparsely coated with mineral particles and at no 
stage there is a deficiency of the bubble surface. low it is 
likely tnat the conditions and forms of mineralized bubbles 
remain constant and, therefore, the first order kinetic equa- 
tion describes the process during the entire period of 
flotation. In the case of inhibited flotation, bubble rising 
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through the pulp becomes effectively saturated with the 
particles and reaches the surface with a load of particles 
determined either by its surface capacity to accommodate 
particles or its lifting power. A generalized equation 


?/hich includes both these cases was given by Pogorely^^ 


_dq 

dt 


Sq'V L 1-exp (. 


Ka_C T 


')] 


(2.5) 


where q is mass of solids in froth product, is the ratio 
of total available bubble surface to bubble surface occupied 
bj. oOlid particles in unit time, y is the amount of parti- 
culate solids on unit bubble surface area, and 1 is the average 
holding time ox bubbles in pulp. For low concentration of 
solids in pulp, equation (2.5) reduces to. a first order process 

while in the case of high concentration, the equation tends to 
zero order. 


■^*2 on— f l oatabl e Compo nent s : 

An assumption implicit in the foregone discussion is 
that the feed solids are homogeneous. I'he experimental data 
when plotted according to n-th order rate law, however, fre- 
quenoly failed to verify the proposed rate expression. Morris^^ 
modified the first order rate equation by incorporating an 
additional parameter C^. as follows; 

= - i^a (C - (2.6) 


dt 
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He did not 'provide any physical interpretation of 0^ hut 
computed it from experimental data. Bushell^^s however ^ 
recognized 0^,^ as the tern^inal concentration j representing 
the material that would not float even after an infinite 


time , 


observed that 0 


nonlinear function of parti- 


cle size. 

' 7 . ‘‘7 

Cooper"'' , 


i>imilar observations have also been reported by 
On the other hands the work of Jowett and Safvi'''^ 


strongly indicabes that is merely an empirical parameters 
whose actual behaviour does not correspond to the true terminal 


concentre cion andj moreover, instead of being a strictly minera' - 
logical Variable it apparently represents the hydrodynamic 
nature ox the test cell as well as the flotation chemistry. 


5 Probabili stic and Stoc h astic Kine tic Models: 

According to Schuhmann^ the kinetics of first order 
homogeneous reaction is applicable to flotation , and the ovcra_ 
probability of flotation of a solid particle may be expressed 
in terms of probabilities of collision (P ), adhesion (P ) and 
particles being retained by, bubble (P^). Thus, if individual 
probabilities are independent of each other, the overall proba- 
bility (P^) of flotation is' 


P„ = P P P 
0 car 


( 2 . 7 ) 


Similar approach v/as adopted by Bogdanov et al^^’^^ A 

modification to equation (2.7) has bee ii proposed by Tomlinson 
and Pleraing ^ 
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p p p 
c a e f 


( 2 . 8 ) 


where P^ is probability of a particle being retained and 
lifted through the liquid and P^. refers to the probability of 
a particle being retained in the froth. Strictly speakings 
the hypothesis of statistically independent events, iniplicit 
in equations (2.7) and (2.8) is not valid. For example the 
probability of particle- bubble attachment depends on the num- 
ber of particles already present on the surface. For a battery 
of continuously operated cells arranged in seriesj Kelsall^^ 
showed that the mass of material, Mrp, leaving in tailings 
after passing in n^ cells in series is given by 


% ( 1 - 


(2.9) 


where is mass flow rate of particulate solids in feed. 

Kelsall also p3xiposed a graphical method for evaluation of P . 

0 

... . 42 , 

lavis has applied this approach to industrial plants floating 
lead and zinc concentrates. However , as pointed out by Woodburn'^'- 
the main objection to this model is that it fails to take into 
account the increase in residence time per cell down a conse- 
quitive cell bank which results in gradual increase in proba- 
bility of flototion. 


Bodzinoy'^'^ has been among the first authors to develop 
a stochastic model for flotation. He showed that the flotation 
can be described as a so called 'simple death process' which 
reduces to the first order kinetics in the mean. Recently, 
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Kat2 and Shinnar^^ have considered flotation as a random pheno- 
m.ena belonging to the family of the Markov processes. 

^ fwo ghase^^ Mod el : 

A kinecic model of the flotation process 5 which recog— 
niaes that the cell contents are partitioned between pulp and 
froth and that there exists a mutual traffic of material 
between these two phase Sj was proposed by Arbiter and Harris^ 
and further extended by Harris and Rimmer^"^. The basic assump- 
tion in this model is that both pulp and froth phases are per- 
feebly mixed 5 and first order kinetics is applicable for indi- 
vidual transfer step. This model was experimentally verified 
by Harris, Jowett and G-hosh^® for continuous flotation of coal 
at steady state. Bogdanov et al^^ and Bull^^ have also supportc ’ 
this approach. later, Ball, Kapur and Buerstenau^ introduced 
distributed flotation rate constant into this model. 

Several authors such as l^Iiemi^^, Woodbum and loveday^^, 
and Imaizumi and Inoue^^ pointed out that for sufficiently fast 
rate of froth removal the rate of return of material from froth 
to pulp is negligible and, therefore, under such conditions 
two phase model is not necessary. This point, however, has 
been criticized by Maksimov and Khainman^^ who observed that 
even for high concentration yields and with the use of reagents 
which ensured a stable froth, the probability of particles 
remaining in the froth was much smaller than unity, and the 
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rate of return of particles increased rapidly mth particle 
. size. Bushell. assumed that the rate of re-entry of particles 
to the pulp ?/as proportional to the rate of transfer of material 
from pulp to froth. One of the hasie assumptions of the two 
phase model, that the froth phase is perfectly mixed is really 
questionable. Gooper's-^"^ concept of horizontally homogeneous, 
but vertically distributed froth phase seems physically more 
reasonable. 

5 Kin _ed :Ju i.|Iodel3 Based on Mechanis t ic Approa ch: 

Models included in this category are based on one or 
more sub-mechanisms involved in the process. Although these 
models lead to result far too complex to be of any industrial 
utility, they, nontheless, attempt to provide valuable insight 
into the likely mechanism of the process. 

Mineral particles are transferred from pulp to froth 
through any of the following four mechanisms : (1) attachment 
of the mineral particles with the air bubbles in the pulp 
follov/ed by the transfci' of bubble to the froth; (2) selective 
precipitation of the gas dissolved in water on the surface of 
mineral particles; (3) mechanical straying of the mineral particles 
into the froth column; and (4)., attachment of mineral particles to 
the bottom of the froth column by the agitation of the pulp. It 
has been observed by many investigators that the third and fourth 
mechanisms of mineral transfer from pulp to froth are not signi- 
ficant in many systems and the major fraction of the material 
is transferred in the form of bubble-particle aggregates.^ 
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56 

Taggart-" proposed the hjrpothesis of huhlDle-particle 
aggregate formation due to precipi'cation of gas on the hydro- 
phobic surface of minera.l particles from supersaturated gas 
solLition. 'This, however , -was criticized by Bogdanov and 
4 ilanovski''*'^ , and Gaudin^^^ who postulated that the direct 
encounter betvi^en pai'ticles and bubbles is more plausible for 
the formation oi bubble particle aggregate. The direct encoun- 
ter h3/potnesis has been supported by Malozemoff and Ramsey^^, 
Plante and Sutherland^® , Spedden and Plannan^^, wbelan and 

Brown^^, Kirchberg and Topfer^^, and Imaizumi and Inoue^'^. 

64 

Klassen' ‘ proposed an intermediate mechanisia in which a large 
bubble coalesces vrith a micro bubble precipitated and attached 
to the mineral surface. Ihe conditions under which gas preci- 
pitation on the mineral surface becomes important have been 
reviewed by Smith et al^^, and ERassen^^. It should be pointed 
out that the phenomenon of mineralization of bubbles by preci- 
pitation of gases raises a nuivber of fundamental questions, 
the ansv^er to which may be of great importance in flotation 
practice. Por example f in the flotation of very fine parti- 
cles, owing to their small mass, particles are normally swept 
away with bubble stream lines and undergo collision with great 
difficulty. In the case of coarse particles also, formation of 
several clusters of small bubbles by precipitation of gas 
may enhance the rate of flotation. , 
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-^he foriDa'cio]! oi bubble iDineical. parbicle aggregates 
by direcb encounters which is the most frequent mechanism under 
noriflal flotation conditions,, involves several sequential stages, 
viz.^ adsorption of reagents on the surface of pa.rticles and 
bubbles, collision between bubbles and particles, attachment, 
partial detachment of particles from the bubbles due to gravi- 
tational and Inertial forces, and froth un stability etc. The 
kinetic models developed on the basis of mechanistic approach 
incorporate one or more of these steps, 

'first such model was proposed by Sutherland'* . He 
assumed that water is non viscous and incompressible, the 
bubbles and particles are both rigid spheres and incompressible, 
and the inertia of particle is negligible, and calculated the 
probability of interception from streamline flow around spheres* 
The collision radius, r^, of the bubble for a particle of 
diameter 1, as obtained by him, is . 

2 5 ® ^ ( 2 . 10 ) 

where a is the diameter of the air bubble. Photographic 
observations by .Spedden and Hannan^"*, however, showed that the 
larger particle deforms the bubble appreciably during collision. 
This deficiency 'ms removed by Philippoff and Evans, who 
independently postulated a bounce theory in which the role of 
inertial forces is emphasized, further modifications to the 
theory have been proposed by Whelan and Brown^^. On the 
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other hand the empirical correlations reported by iangmuir^^, 
Herne \ Rana and \/ong'^\ on collision of particles seem to be 
consistent with the results of Sutherland. Using a similar 
approach, Bogdanov and his collaborators^®"'^'^ have calculated 
the probabilities of collision and adhesion. Recently, Flint 
and liowarth ~ evaluated the collision efficiency of a particle 
with a bubble from considerations of particle trajectory as a 
function of particle size and hjrdro dynamic conditions prevail- 
ing in the cell. Although 'this stud 3 r provides a qualitative 
understanding of the collision phenomena, its extension in 
quantitative sense, to a system of swarms of particles and 
bubbles in turbulent medium is not available at present. 


A rate equation based on analogy between adsorption and 
mineralization of air bubbles "vas derived by Matveenko'^^. A 
detailed analysis of bubble particle interaction for small and 
rnediuiTi size particles (< 100 /-i ) , based on physico-chemical 

r-7 A 

aspects, has been made by -Derjaguin and Bukhin'^. The mecha- 
nism of bubble particle attachment has also been studied by 


Sieeles 


^ 5. '6 


Morris 5 ZaidenbersC"^ and Glembotski 


80 


et al 


81 


A more detailed and structurally complex phenomenolo- 


gical model for flotation kinetics has been proposed by Mika 

S2 

and Fuerstenau wherein each elementary rate process is 


formulated in terms of linear first ox'der differential equations 
for individual discrete size particles and bubble species. 
Although this model clearly identifies and elucidates the basic 
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stup,j in interphase- transfer of particles, its practical 
util j. rcicains to be demonstrated because it involves a 
nuino .,,1 oj. unknovn parameters \';hich are at present not amena- 
ble to theoretical or experimental estimates. Recently, 
Vi'oodburn, king and Golborn®^. have proposed a kinetic model 
for a continuous flotation cell in viiich the rate constant 
ic derived as a function of particle size, its mineralogy and 
the tailings flov/ rate. The final expression is product of 
three probabilities pertaining to events of Collision P 

c ^ 

adhesion P^ and disengagement P^ as iollo?/s: 


P^ = B exp (- ) 


2 


a 


( 2 . 11 ) 
( 2 . 12 } 


■D 

^d 



max 


for 1 < L. 


max 


= 1 for L>L (2.1'5'i 

<^/here B is a pi'oportionality constant, S is a hydrodynamic 

paramei-er, L is particle size, x. refers to contribution of 

surfac:. hydration to the probability of adhesion, and I is 

" ■ max ■ 

maximum particle size Miich will remain attached to a bubble. 
It should be noted that in the derivation of these probabi- 
lities, a number of simplifications and crude approximations 
had to be made and at best, these represent a first tentative 
attempt in this direction. 
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^ Models with D .i e d 

iia t e Cons tant ^ 


App arent Flotati on 


In the literature cited above, it is assumed that the 
rate constant is single valued. It is, however, more likely 
that bhe feed material comprises of several different species 
with different values of the rate constant associated with each 
of these species as well as with surface and composition hetero 
geneity within a species. Thus, it is more rational to asso- 


ciate a distributed rather than a single valued rate constant 
with che feed material, The over simplified probabilistic 
model of kelsall^”^ reflected this fact in that there are at 
least two species of floatable particles, an easily recoverable 
frac-Gion having a high probability of flotation, and a low 
probability fraction which floats with difficulty. 


I^Oo^^sly 5 end Imaizumi and Inoue^^ wei'c among the 
first authors to formally introduce the concept of variable 
and distributed rate constant. The latter postulated that the 
rate of removal of solid from pulp to froth is a time invariant 
firs'c order process but the particulate species are distributed 
with respect to the overall specific flotation rate constant. 
Thus, if lig (iia, t) is mass of solid particulate species with 
attribute overall rate constant Ka, remaining in the cell at 
time t, the rate of removal of solids of this specification 
from the pulp is 
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fflvIgdCa, t) 

^ , ±) , 0 < la < oo (2.14) 

Solution to this eq_uation is 

Mg(Sa, t) = Mg(Ka) exp [-Ka t] (2.15) 

where , 

MgdCa) = M3(iCa, 0) ' (2. 16) 

Therefore 5 the total material reniaining in the pulp '.at time 
t is, 

CO 

i'-gltj =/ Mg(Ka) exp [-Ka. t] dSa (2.17) 

0 

If Mg(t = 0) is normalized, M^^(Ka) "becomes a true density 
function s and. 5 explicitly, Mg(Ka) dKa represents the mass frac- 
tion of particles with overall rate constant in range Ka to 
Ka + dKa in the feed material. The main difficulty in this 
model lies in evaluation of i«fg(Ka) given noisy experimental 
data Ivlg(t). If Mg(Ka) is a discrete density function, equation 
(2.17) becomes, 

Mg(t) = Mg(Ka^) exp [-Ka^j t] + Mg(Ka2) exp [-Ka2 t] 

+ ... +B5g(Ka^ ) exp [-Ka t] (2.1S) 

s s' 

where .. 

Mg (Kai) ' = , 

' X— T 


1 


(2. 19) 
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Imaiz-urai and Inoue^^ proposed a method for evaluating Ka^ 
and corresponding Mg(Ka^) by considering asymptotes on the 
log (Jiig(t)) versus t plots. Similar approach has also been 
taKien by Hiioer— Panu Unfortunately, the number of species 
^s 1 ^ 0 1 kiioTO a priori and must be assumed arbitrarily. 

Comple ue identification of the disc-rete distribution from noisy 
data is a very difficult task.. A more commonl,y used method 
of determining the density function consists of assuming a 
funccional form for Mg(Ka) in equation (2.17), and then deter- 
mining the parameters of the distribution from the experimental 
data, thus, Woodburn and- Loveday^^, and Inoue and Imaizumi®^ 
assumed a gamma distribution and evaluated the parameters by 
graphical methods, loveday , however, observed that a better 
fit to experimental data is obtained with incomplete gamma 
distribution. Ball, ICapur and Puerstenau'^'^ also used gamma 
distribution for Ivlg(lia.) and carried out the analysis of the 
flotation kinetics using the moments of the dis’fcribution. 

.QP 

Harris and Chakravarti'^- suggested a bimodal gamma distribution 
for Ivig(iCa). Kinetic equation by Huber-Panu®^, however, starts 

with a rectangular distribution of rate constant. When M„(Ka) 

s 

is assumed to be an incomplete or a bimodal gaimma distribution, 
areater flexibility and enhanced agreement is to he expected 
because now M.,(I^a) is expressed in terms of three or more 

■O' , 

adjustable parameters instead of two in the case of gamma 
distribution. Hecently, Black and Paulkner^^ have arbitrarily 
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assigned 'values to n and Ka- and evaluated M (lCa. . ) in equation 
(2.18) by linear regression techniques. This is essentially a 
curve fitting e3:ercise and provides no insight into the structure 
of the distribution. A large number of different equations can 
fit the data quite accurately and therefore, evaluation of rate 
constant distribution from equation (2.13) in this manner is 
rather meaningless. Methods for analyzing kinetic data assuming 
distributed rate constant models have also been discussed by 

Woodburn, King, Buchalter and Piper^”', Huber-Panu^^, Tille and 

9 1 9A 

Panou and Outtris . 

riarris and Ghakravarti '' correctly pointed out that in 
equation (2.17), M (t) is the lapl 3 .ce transform of M_(Ka) and, 
therefore, Mg(Ka) is inverse of the Laplace transform. Hence, 
the prcoblem of evaluating flotation rate distribution, Mg(Ka), 
is the problem of inversion of the Laplace transform. In the 
case of experimental data Mg(t), this inverse 'must be computed 
by numerical means. Unfortunately, m-imerical inversion of the 
Laplace transform is not an easy task and success is not always 
assured, especially, in presence of experimental errors. This 
is because Laplace inverse is a very unstable operator which 
implies that arbitrarily small variations in M„(t) may cause 
arbitrarily large changes in M^(Ka) as highlighted by Krylov and 

S'- 

OK 95 

Skoblya , and, Bellman Kalaba and Lockett , Although in the 
literatuire a nuraber of algorithms are available for numerical 
inversion of the exact data, they invariably become ineffective 
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in the case of noise corrupted data. Some of these algorithms 
are discussed in Appendix A. In Chapter 5, we shall present a 
simple and robust method for the numerical inversion of the 
Laplace transform. 

it should be understood that although Mg(lCa) is a fre- 
quency function, it is not an exclusively internal coordinate or 
attribute of the particles. In other words, the apparent rate 
constant, ICa, depends, as pointed out by Harris and Chakravarti®^, 
on numerous factors only a few of which, such as sise, degree of 
liberation and distribution of valuable mineral in a particle,, 
its surface characteristics and conditioning etc, are specific 
attributes of the particles. .Other factors that contribute to 
the structure of the apparent flotation rate constant are aera- 
tion rate, population, size, and residence time distributions 
01 bubbles, as well as, complex hydrodynamics and physico- 
chemic.,,1 conditdons prevailing in the pulp. I'la incorporates in 
some complicated fashion a number of elementary interphase 
transfer steps such as particle bubble attachment and detachment 
rate, and escape of particle ladden bubbles from pulp to froth. 

iforic of Laidenberg, Lisvoskii and Burovoi^'^ is notevvorthy 
in this connection. They have attempted an analysis of flota- 
tion kinetics using a multivariate distribution of the particle 
entities in particle size, degree of liberation in individual 
size fraction, and floatability for a given size and mineral 
composition. However, the main drawback of their model is that 
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invariably neithc-r the conditional distribution in mineralogy 
for a i^^iven size nor tne floatability as a function of particle 
size and mineralogy is knovn. Hov/ev or, recent results of 
"hoouDUTij, King and Colbom'~^ have shovn that even a limited 
e-'^tenbioxi oi firsfc order treatment on the liiaos suggested by 
Eaidenberg et al offers significant advantages for simulation 
and control of the flotation circuit. Recently, King^® has 
furcneo. extended anct applied this mu3-tivariate distribution 
model for analysis of kinetics for a flotation circuit. 

flotation Cell: 


liost of the models discussed till nov/ v^ere developed for 
semi-batch flotation operation. Although these can be easily 
extended to steady state continuous case, the analysis based on 
semi-batch system, as such, has been criticized by many investi- 
ga'uors on the ground that both the mineral and reagent concentra- 
tions undergo continuou.s change with tiiiie, thereby, possibly 
enhancing tiie complexity of the system. Further^ most of the 


information has to be collacted in a r. 
lead to erroneous results. Ihe merits 


ther short time v/hicb may 
and demerits of both 


semi-batch and continuous operations for development of flotation 
kinetic models have been discussed by ilrbiter and Harris^, 
Woodburn et al^”*, Brown and Smith’’^, and deBruyii and Modi^^. 
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The basic difference betv/een sa-i- batch and continuous 
operations is tnat in the latter case the residence time of all 
elements in the system is distributed. If in a continuous 
steady state cell with spatially homogeneous pulp the residence 
tiixJe distribution of solids is S(t) it, i.e.^the fraction of 
solids in tailings exit that has spent time t to t + dt in 
pulp 5 tne tailings isass flov/ rate of particles of size 1 with 
overall rats constant ICa. is^^ 

CO 

M^(Ka, W = f I4p(Ka,L) B(t) exp f-ICa t] dt (2.20) 
o 

where (lia^L) is mass flov/ rate in feed of particles of size 
1 and overall rate constant Ea. In distributed rate model the 
overall flow rate of particles of size L in tailings, Mj,(l), 
and concentrate streams, Iti^Cli) , are, respectively, 

oo CO- 

Mr^^(l) = f S Mp(Ea,l) E(t) exp [-Ka t] dt dKa (2,21) 

0 0 

and 

, oo , 

M^(I) = / M (Jfa,l) dKa - (2.22) 

, O ' ' 

Many investigators have postulated that a continuously operated 
sub-aeration flotation cell is close to a perfect mixer^^’®^’^"*’"^''”^''^- 
and therefore, E(t) is an exponential distribution function. 

However, Niemi showed that although the residence time of 
liquid phase can be assumed to be that of a perfect mixer, in 
the case of solid particles, particularly of large particles, 
the assumption of perfect mixing requires some modifications. 
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He further reported that even if solid particles were apparently 
ideally mixed, Laeir mean residence time was different from 
that of the liquid phase. 

Bogdanov, Hainman and Maximov^® have attempted to 
explain the observations of Bogdanov, Kisivalter and Khainman'^^'^, 
and Zakhvatkin et that the time of flotation required for 

a given recovery is reduced ?/ith increase in throughput. They 
deduced that with increase in throughput the flotation time is 
improved due to changes in the hydrod3rnamic regime of the pulp. 
However, analysis of the data of Yfoodburn et al in Chapter 7 
of the present work shows that the throughput rate has negli- 
gible influence on the apparent flotation rate constant. Bogdanov 
Hainman and Maximov^'^ derived a rate equation ' for an ideal case 
where there is no mixing and cell behaves as a plug flow 
reactor. In such a case flotation is predicted to be indepen- 
dent of the flow rate, loveday^'^ reported that reduction in 
the mixing coefficient is likely to give improved grades at a 
constant recovery level, results which are also well known in 
the theory of continuous polymer reactors. 


In a significant contribution which includes valuable 
and extensive experimental date, Woodburn et al'’^'^ have shown 
that the residence time distribution is adequately characterized 
by the gamma distribution 
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wnerc ch© pai’aiUG’Gers a and p are 
rate and particle size and massj 
variance of the .distribution by- 


functions of slurry flow 
and aro related to mean and 



( 2 . 24 ) 



(2.25) 


ihe additional parameter a? vvhich was found to be quite small, 
takes into account the deviation from perfect mixing behaviour 
(a = O). Substitution of equations (2.25) and (2.21) into 
(2.22) leads to the xollov/ing expression for recovery.' 




(L) = 


/ 

o 


IvIj,(Ka 


, 1 ) [1 


,a+1 


(K3. + 0 


a+1 


] dKa 


( 2 . 26 ) 


M|Q£Li^^ocess 'triable s on^^ 

S.® Particle size is one of the most ■ 
important parameters and has been studied extensively, among 
others, by iCidd and Wall,’^'^^ Sutherland^, Chuan and Zimmerman’’®^, 

•Tu- • 35 . ■ *ir,p 

morns , uaudin, Schulimanii and Schlechteri -Bo,<'danov and Khain- 
man Sigeles and Leviush"^®-^, Brown and Smith deBruyn 
and Modi Bennett et al I'hrokhin”'"^ , Mackenzie and Matheson^"*^, 
Huber-Panu and coworkers’’''’^"'''’^, Tomlinson and Pleming^"^, 

Imaizumi and Ihoue^^, Heinrich and Walter Mika'*^, Ek’’'^^, 

Tevari and Biswas'*'’®, Suwanasing and Salman ^^,lbodburn et al 
and Clement et al”’”’^. Prom' these studies, it may be concluded 
that regardless of the mineral, collector or flotation machine. 
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there is a definite trend with regard to variation in the 
recovery as a function of particle size, and there exists an 
intermediate size or size range at which the recovery is 
optimum. This size depends on the rxiineral to be floated and 
the hydrodynamic conditions in the flotation machine. Both 
fine and coarse particles float with difficulty, and in the 

latter case practically no flotation occurs above a limiting 
size. 

The explanation for low recoveries of very fine parti- 
cles is attributed to low probabilities of collision because 
these have a tendency to follow the stream lines. On the other 
hand, the probability of detachment increases with particle 
size because of disruptive and inertial forces that come into 
play, ihese opposing influences of the particle size give 
rise to a critical intermediate size at which the flotation 
rate is optimum. Combining the probabilities of collision, 
adhesion and detachment given in equations (2.11), (2.12), 
and (2, 13) ffoodburn et al^^ have, recently, proposed an 
expression for the overall rate constant IQa as a function of 
particle size, L, and hydrodynamic parameter, 0, as follows: 

ICa = k(G) r®^p(- -I)] [l-( 1 )1-5j ( 2 . 2 ,) 

Ii ■ max 

where k(G) is a function of the degree of hydration of a parti- 
cle of grade G and it depends on the mineralogy of the 'oarticles 
and its surface treatment with the aurfactants. 
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(ii) J 3 X-i^£a t io_n^ I'he effect of aeration rate on 

flotation kinetics has been studied by Huber-Panu et al 


114,120 


Inoue and Imaisurai®^, fevmri and Biswas^ 


lynch 5 Edwards 


Burden and Bennett et al^°. The first three groups of 
investigators reported that the overall flotation rate increased 
uitij uLe enhanced rate ou aeration. Mathematical model proposed 
hy king also supports this observation. Edwards, however, 
observed that althouga the flotation kinetics was very sensi- 
tive to the rate of aeration, the efficiency of separation 
remained more or less unchanged. Bennett et al pointed out 
that the enhanced aeration rate could lead to increase in 
flotation rate only if the distribution of bubble size does 


not change because the maximum flotation rate is achieved with 
largest nuiiiber of small bubbles = lynch, however, did not find 
a significant correlation between flotation kinetics and 
aeration rate. 

(iii) Pulp D ensit y : The effect of pulp density on 

flotation rate has been investigated by Fewton and Ipsen”^^^, 
uutsalyuk"'^' ' 5 iJ.assen'^^'^, Burciona^^, Imaizumi and Inoue'^^, 
deBruyn and Modi^^, and loveday®^. They observed that except 
for highly concentrated suspensions the average or instantane- 
ous flotation rate is independent of pulp density, v/hich is 
in conformity with the first order kinetics in free flotation 
regime. In the case of concentrated pulps the flotation 
depends primarily on the capacity of the carrier air phase to 
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lift the Particles from pulp into froth rather than on the 
particle “bubble attachment rate constant, and therefore, flota- 
tion rate constant is no longer independent of the pulp density. 
Now the specific rate constant in terra of per unit mass should 
decrease with increasing concentration of solid particles in 
the .pulp 5 as shov/n in Chapter 6 of this study, 

( iv ) Bf fe^c tp_of^ he.aSejSi.§.* 

Collectors: Both the structure of the collector molecules 

and its concentration affect the flotation rate by influencing 

the particle-bubble attachment rate. In general, the adhesion 

of particle to bubble improves with increase in the length of 

hydrocarbon chain, in the collector molecule. The effect of 

collector concentration on the flotation kinetics has been 

studied by Beloglazov^, Schuhmann^, Inoue and Iniaizumi®^, and 
0*1 

lynch . These studies shov/ that both the recovery and flota- 
tion rate increase with concentration of collector in the pulp. 
Hoivever, this is usuall;/ aoconrplished with loss of selectivity. 

In order to maxi .ii?e recovery, most industrial flotation circuits 
are operated at .saturated colloctor dosage. 

Brothers: In addition to providing stability to air bubble, 

the function of the frother is to decrease the bubble size 

' 1 ' '■ " 

and bubble velocity. Studies by Grunder et al , and’ Bennett 
et al^*^ have shown that increase in the amount of frother tends 

to increase the number of bubbles by reducing their average size, 

" ' ■ 121 

and the flotation rate is enhanced, lynch , however. 
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att-ribu ced the xi'icreased id-otation rate with frother addition 
to enhanced rate of water transfer from pulp to froth. 

It may he pointed out here that the effect of reagents 
on flotation kinetics has not been discussed here in detail 
because in order to do .justice to this 'topic it will be 
necessary to treat individual mineral S3?'stems5 which cover 
a Vast area. 

B. OPIEIIZAIIOd OF FIOTAIION OPERATION; 

Once the kinetic behaviour of a single flotation cell 
is Icnomc in terms of its process variables 5 it is then possible 
in theory 5 to optimize the process efficiency and automatically 
control the process'* . The problem of optimization, 
considered so far essentially consists of defining a suitable 
objective function, followed bj'- determination of the values 
of these process variables such that the objective function is 
optimal. The problems associated with the formulation of the 
meaningful objective function, avppropriate to flotation process, 
have been concisely stated by Nfron and Moir'*'^^. Jov/ett and 

-I ■7'Z, 

Ghosh optimi.zed the performance of a flotation circuit 
consisting of a sequence of cells with series tadlings through- 
put and parallel froth withdrawal. Their profit fimction 
(objective) was defined as the difference between the selling 
price of the concentrate and the production cost which included 
the cost of the raw material and the cost of operation. 
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Sergee\7 and Smimov also used a similar objective function 
in which, however, the number of cells connected in series 
was predefined. E.ecently, V'/oodburn et al^"* have proposed a 
simpler objective function where the recovery of concentrates 
was maximized such that the grade of the product was always 
better than a minimum acceptable grade value. 

In industrial flotation practice due to various reasons 
the separation of valuable mineral from the associated gangue 
is Seldom complete and in order to improve the process effi- 
ciency it is a noimal practice to employ a number bf flotation 
cells or banks of flotation cells. The basic drawback of all 
the above mentioned approaches is that the configuration, i.e., 
bhe number of cells and/or their linkages is assumed a priori. 
In general, a flotation circuit is comprised of many flotation 
cells, performing rougher, cleaner and scavanger functions 
and it is q.uite tedious to enumerate all possible meaningful 
configurations. A particular configuration which is optimal 
for a given type of ore under certain process conditions may 
not be optimal for other ores. In general, there is no 
justification for assuming a configuration ana then optimizing 
the process variables, since the assumed configuration may not 
be the best possible. Thus, the problem of optimization 
should be conceived of comprising the follov/ing two steps: 

(l) synthesis of flotation circuit , that is, hovr the various 
cells are interconnected such that the optimal process 
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efficiency is embedded in the derived circuit, and (2) deter- 
mination Ox the Values of the process parameters of this 
configuration. 

nach flotation cell can be considered as a separate 
reactor unit. Recently, systematic strategies have been 
developed for synthesis of the chemical reactor units for 
optimum efficiency which are relevant to the flotation prob- 
lem. ‘I'hese are briefly reviewed in the following section. 

Synth e sis of R e actor Units: 

Except for the isolated treatments of the heat exchange 
networks and multicomponent distillation separation sequences, 
the development of formal systematic approaches to synthesis 
in the design of chemical processes dates back to 1968 only. 

The subject is still in considerable flux. The chief diffi- 
culty encountered during the process structux’e synthesis lies 
in the large number of possible configurations for the typical 
industrial processing tasks. Examination in turn of each of 
these configurations is not only tedious but occassionally it 
is not possible to enumerate all possible configurations. 
Conventional optimization techniques designed to hardle dis- 
crete Variables (representing equipment connections) are 
generally applicable to fairly small size problems only. 

135 

Rudd proposed an approach to process system synthesis 
which is based on decomposition of the main problem into a 
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seq.uc;nce oi subproblecis until the level of availaole technology 

is reacned. 1‘his approach has been applied by ilishida et al'^^*^s 

and Kobayaski c^t al ^ for optimal synthesis of heat exchangers ^ 

i 3 8 

by kenzies et al for synthesis of energy recovery networks , 
and by Umeda and Ichikav.'a^ for synthesis of generalized systems 
Masso and Kudd^^*^ further modified the decomposition approach 
by inoroducing the use of heuristic problem solving techniques. 
Ihe combined heuristic design decomposition approach to process 
analysis has been applied for synthesis of chemical processes 
by Siirola 5 Powers and Rudd''^''^ and Rishimura and Pliraizumi 
etc. It mayshowever, be pointed out that the final configura- 
tion depends heavily on heuristic rules which are used at the 
Various decomposition stages. As Masso and Rudd correctly 
pointed out, really satisfactory rules are yet to be developed. 
Iherefore , solution to practical problems by this technique 
is fairly limited at this stage. 


A synthesis technique, knoTOT as 'Evolutionary^ synthesis', 
has been suggested by King, Gantz and Barnes and Ichikawa 
and Ran ’’ . Ihis technique is based on synthesis of nev/ process 
by modification of previously generated techniques. A somewhat 
more theoretical approach to the evolutionary sjTithesis was 
taken by McGalliard and Wes t er berg”* who applied Lasdon 
dual feasible decomposition method for process optimization. 

The strategy for optimal design of chemical processes from the 
vie?/ point of statistical design theory has been studied by 
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Watanake, I'Tishiiiiura and Matsu'bara”'^^. The branch and boimd 
techniq.ue 8 of problem solving have also been shown to provide 
guidance iii solving the design problems. As shovni by lee 5 
Masso and Kudd , the design problem is branched to simpler 
design problems v/hich bound the original problem. The main 
difficulty 3 however 5 comes in formulating the alternate prob- 
lemj solution of which bounds the solution of the original 
design problem. Umeda, Hirai and Ichikawa*'^’" have attempted 
to apply v/ell known optimization and mathematical programming 
techniques to the synthesis of chemical processes. In their 
approach all possible configurations are imbedded into a 
generalized configuration by defining all the interconnections 
which raight exist betw/een various reactor units. Split frac- 
tions with Values zero to one are assigned to all such inter- 
connedtions and the solution to the s^/nthe sis problem is obtained 
by determining optimal split fraction (O or 1) for each linkage 
along with the optimal design paraxxie-ters. -Values of one or 
zero indicate whether a pair of equipment are connected or not. 
Shis ■ approach has been criticized by Hendry 5 Paidd and Seader'^^*^ 
on the grounds that the parametric space for most of the indus- 
trial situations becomes so large that it is difficult to 
ensure a global optimum in the resulting mix integer optimi- 
zation problem. Recently, this' problem has been taken care 
of to some extent in some systems by a multilevel' technique 
which partitions the integrated problem, into smaller 
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independent subproblems which can be solved with ease, lo 
coordinate the solvitions of subproblems and determine the 
overall system an adjustment problem is required to be 


carried out. 



CtLAPTSR 3 


OUTLIiiE OR WORK 


The present work can be divided into t?;o parts. In the 
first part some aspects of flotation kinetics have been studied 
whereas the second part deals with the synthesis and design 
of a flotation circuit. 

Part I; RLOTATIOH KBiSTICS: 

Shis can be further subdivided as follov/s: 

( 1 ) .4 .P henp ineno logical M odel for Rlota tion Kine tics: 

From the literature review given in the previous chapter 
it is evident that the kinetic models based on distributed 
apparent flotation rate constant lack justification for the 
first-order kinetics wdth respect to the species with a given 
rate constant. 1 phenomenological model for flotation kinetics 
has been derived which is perhaps more meaningful in so far 
as the model provides an explicit relationship between the 
apparent rate constant and the specific rate constant for 
particle bubble attachment suhprocess which, from the flota- 
tion mechanism point of view, is a more pertinent attribute. 
Concurrently the effect of bubble population is incorporated. 
Ivlatherocitical conditions for the free flotation case under which 
flotation process can be described by a simple first order 
rate equation given by equation (2.17) are stipulated. 
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^ i i j'lotatloi'i Ra t e -Djstribut ion_ _b y Hign erical 
Inv er s_iqii_ the__L a. plac e Trans form : 

The problem of determining the density function M_(Ka) 

O 

in equation (2.17) has been solved, so far by assuming a 
priori a specific freq.uency function (usually, the gamma 
function) and then evaluating the parameters of the distribu- 
tion from experimental data by some graphical techniau;e. These 
graphical methods are quite tedious and the results depend, 
to some extent, on the personal judgment of the computer. 
Although recently, Woodburn et al"^ have used numerical methods 
for estimating the parameters, a more fundamental objection, 
ho?/ever, lies in the rather arbitrary assumption of the form 

of the distribution function for M (iCa). Harris and Chakra- 

s 

86 

varti' have correctly pointed out that in equation (2.17) 

M_(t) is the Laplace transform, of I\i (Ka), therefore, M^(Ka) 
is inverse of the Laplac-. transform. A numerical method for 
estimating the r; ta constant distribution in cumulative mode 
by inversion of the Laplace transform has been developed. The 
method has been tested for synthetic data generated from 
different distrihiitions. Effects of errors in data, number 
of data points, initial guesses in the optimization step and 
non-floatables have been studied. This method of inversion of 
the Laplace transform has been found to be especially suitable 
for noise corrupted data which is inevitable in flotation 
experiments. 
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( iii ) TAe^l£ecM.pf_.Aerati0]a Rate, Part icle Size and ?u1p 

-'fc.he Flotat ior i Rate Sistributioiis: 

In the literature cited above, the effects of process 
variables on flotation kinetics were determined either in 
terms of the rate consrants in some empirical kinetic expres- 
sions or b.v monitoring the instantaneous flotation rates. 

Q/C 

Perhaps the only exception is the study of Inoue and Imaizumi 
who studied the effects of the aeration rate on the distribution 
as a v-hole. In this part of the v/ork, using the numerical 
inversion technique, developed in part (ii), cuQiulative 
apparent flotation rate distributions have been evaluated 
for kinetic data available in literature. Appropriate rela- 
tionships between the estiraa.ted distribution functions and the 
aeration rate, particle size and pulp density, respectively, 
are examined and a partial proof for the validity of the model, 
developed in part (i) is provided- 

( iv) f,he__ Effjcjts^, PA Particle Size^ ana Fe ed, , Aat.GL.pn. .thg . 

Plotc tion Rate Pistrib uti ons^ AA PaII • 

Pirst a reliable method, sjjr.ildr to the nimerical 
inversion of the Laplace transform, for computation of Kj,(ICa,L) 
in expression for mass flow rate of particles of size L in 
concentrate stream in the steady state, given in equation (2,26), 
is developed. This method is again not restricted to any a 
priori assumption concerning the mathematical form of the 
distribution function. Using the tabulated data of 
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WoodlDurn et al the effects of particle size and feed rate on 
rate distribution have been exaiained. The estimated rate dis- 
tribution is then employed to predict the kinetics of a cleaner 
cell whose feed comprises of the rougher concentrate. The 
analysis has been extended to gangue component for which data 
is also available. 

Part II: OPTIlVlii MB SlTBOPTnvUIi SYITTHESIS AMD PUlSia-K OP 
PLOTATION CIRCUITS: 

The integral synthesis approach proposed by Umeda , Hirai 
and Ichika?/a 5 has been applied for S3mthesis and design of 
an optinal or suboptimal multicell flotation circuit, using 
appropriate objective functions of recovery, grade and profit. 
For illustration purpose a flotation circuit with two cells and 
feed comprisiLig of two species is considered. The validity of 
the method is checked by compariiag the results with the real 
optimal circuit obtained by examining each configuration separa- 
tely. The approach is then extended to a more realistic flota- 
tion circuit with four cells and for batteries of flotation 
cells with a three species feed comprising valuable, gangue 
and middlings. 



PART I 


CHAPTER 4 

A PIiSEOMMOlOGIGAL MOREL POR PLOTATION KINETICS 

The model developed in the present work takes into 
accoiint the two distinct states of particles in pulp, namely, 
free particles suspended in water phase and particles attached 
to air buhbles. In first instance a general formulation 
explicitly incorporating the important particle-- bubble attach- 
ment subprocess is constructed which is then reduced to the 
expression in equation (2.17) under clearly defined conditions. 
In this manner it v/ill be possible to identify the relation- 
ship that exists between the apparent rate constant Ea and the 
specific rate constant, K, for particle bubble attachment 
subprocess, in terras of a weighted mean of the bubble popula- 
tion size distraibution. 

MATiiH.iATIC.lu MOiOEL TOR ILOTATION KHSTIC3: 

The model is based on a number of restrictions and 
assumptions as follows, some of which are only approximately 
true, but nevertheless, have frequently been invoked in the 
analysis of flotation kinetics. 

1. The bubble population and size distribution in pulp is 
time invariant-^ i.e. steady state exists. The residence time 

RX Q7 Q-| 

of the bubbles conforms to that of a perfect mixer''^ s - 5 
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'Thusj the escape probability of a bubble from pulp to froth 
is a function of air hold up and aeration rate only (both 
assumed to be constant), independent of the pulp density and 
particle load on the bubble. The latter stipulation is perhaps 
strictly valid for the free flotation case only^"^'. The gas 
phase dispersion does not change^*'' ,i.e bubbles neither 
coalesce nor break-up during their sojourn in the pulp. 

2. The net rate of attaoimient of particles to bubbles in 
pulp is a first order rate process with respect to mass of 
particles in suspension, but the particles are distributed in 

1^5 O 

specific attachment rate constant which is assumed to be 

152 

independent of bubble size and population in the pulp . 

The rate constant for detachment of particles from bubble is 
a function - albeit in some inverse sense - of the attachment 
rate constant as well as the particle size. In pax’ticular, 
attachment and detachment rate constants are independent of 
the aeration' rate in the highly turbulent regir'e. 

3. A'uxiliary interphase particle transfer process, such as 
direct entrainment from pulp to froth and its reverse, are of 
second order significance and can be neglected . Inclusion 
oi these processes will add to complexity of the kinetlc^ ^ 
formulation, but v/ill not materially alter the basic frame work 
of the model. 
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4, The model presented here refers to a semi-ha.tch process 
with fast ±rot:h removal rate, hov/ever, it can he extended to 
the continuous mode of operation, including netY/orlc of cells, 
with appropriate modifications for particle residence time and 

O'!? 

tailings flow rate as suggested hy Woodhurn et ,sl 


let l_a-st/i:j he specifications of an air hubhle of size 
a in the pulp at time t from start of the flotation process 
which was inti*oduced in the pulp at timeT'(t >T ), and 
m(K5l,a,t 5I:') be mass of particles of size L with specific 
attachment rai-e- constant K on the bubbles of specification 
[a,t,T^. Tbe iotal mass of particles adhering to the bubble 
[ajtjTj is 


O 00 


m(a,t;r) = I / m(K,L,a,t,T) dK dL 


(4.1) 


‘1 


where I ^ and Ig smallest and largest particle 

sizes, respectively. It is recognized that the capacity of a 
bubble to accommodate more particles on its surface is an 
exceedingly complicated function of its size, and, the raass 
and size distribution of particles already present on the 
surface vbich would naturally change wit)-- time. This difficulty 
legis to the concept of an instantaneous potential bubble 
capacity mCa,t,r), defined such that the rate of attachment of 
particles of attribute [K,!] on the bubble [a,t,‘i^ is 

m(a,t,r) - ^ . s 

^ (K,l) = aU) [1 -:r--;-w3 K 


( 4 . 2 ) 



where r_^(lv5L) is the rate of attachmert of particles [KjL] 

on a bubhle Mg(K,L,t) is mass of particles £^5!)] in 

free suspension at time t, and A(a) is a function of the 

bubble size and hydrodynamic conditions, including the 

intensity of mixing, of tee pulp* Bquetion (4»2) is similar 
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to the relationship suggested by Sastr^^ and 1‘uerstenau 
v/ith the difference that the bubble size is incorporated in the 
form of a general, but unknov^m function A(a) and m is a function 
of time. Similarly, the rate of detachment of particles [K,1] 
from bubble [ajt,'!] to free suspension state is proportional 
to the mass of particles on the surface as follov/s: 


r^(K,L) = D(K,L) m(K,3j,a,t ,T) 


(4.3) 


where r^(iv,L) is the rate of detachment of particles [K,lJ 
and DCKjL) is the detacMent rate constant. I'he net rate of 
attachmeiit of particles [K,Lj on the bubble £a,r,’r] is 

dm ( u , Xj , a , 5 ) 
dt 


?(K,Jj 5a , t , 


= rjK,L) - r^(2,l) 


^ m(a,t,T) 

A(a) [1- ;t--’ ] K Mg(K,L,t) 

ra(a,t,i") 


-D(K,L) m(K;,L,a,t,T) 


(4.4) 


and the net rate of attachmont of all particulate entities 


on bubble £a,t,T is, 
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Substitution of equation (4.4) into (4*5) gives 

m (a t ^ 2 

r(a5t,T) = A(a) [l- -j f " f g dK dl 

in(p.,t,T) ' Is ^ o 


-f f BdC^L) mdC.L.ajt^r) dil dL (4.6) 

L. 0 


Let iqCajtjX’) "be the nucnlDer of bubbles of specification [ad-sfl 
in the pulp. The mass of particles on these bubbles is 


( Kj L jU 5 1 j'j’) — I'i •j^( a 5 1 , ) ni( Kjl 5 a 5 1 5 it) ( 4 * 7 ) 


Therefore 5 the mass of particles on bubbles of all sizes 

and ages, i.e.jOn the air phase in the pulp, is 

CO -j^ 

M^(K,l,t) = f f lL^(a,t,T) m(K,l,a,t,f) dT da (4.8) 

00 


Moreover, 

m(lI,L,a,t,t) =0 (4.9) 


Differentiation of eo^uation (4.8) v/ith respect to time followed 
by substitution of equation (4.4) results in 


dl.l, (K,I,t) «= t 

J J 4k (a:,t,T) r(K,l,a dr da 

dt o 0 

t dliT.(a,t5T) ^ ^ ■ 

f m(K.L,a,t,r) dT da (4.10) 

' O o’- V . : dt ' , . 

Nov/ the net rate of change of the mamber of bubbles [a,t,l^ is 
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dxI^Ca ,t ,T) 


d-b 


X I'T-j^(a ,t 5 1 ) +dI^(a 5 t 5 T) 6 (t— T) 


0 < T < t 


(4.11) 


where .X is inverse, of the mean bubble residence time T. Hence 






( 4 . 12 ) 


virhere Q is feed rate of air and H is hold up of air in the 
perfectly mixed cell. ^^(ajtjX') 6 (t~T) is the rate of number of 
bubbles of size a fed at time t, and 6 (t-f) is the delta function. 
But 

N^(a,t,r) 6(t-f) = 0 , r f t 

= N^(a),T= ^ (4.13) 

Oonibining equations (4.10) and ( 4 . 11). 


diM^(ll,L,t) 

dt 


00 "f: 


= / / H^(a,t 5 T') r(K,Lsa jtjT) dT do 


0 o 


CT t 


/ %(svS:.,t) m(K,l,a,t,r) dT da 


o o 


t 


+ / / H^Ca jt^T) 6(t_r) m(K,L,a^ 

o ' 0 

( 4 . 14 ) 

In view of equations ( 4 . 9 ) and ( 4 . 13 ) the last double integral 
term vanishes. Further , substitution 01 equation ( 4 *S) into 
equation ( 4 . 14 ) results in 
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di>'L ( r.jij jt) °= , t 

= - XlL (K,Ljt) + / / IT (a,t,T) 

o o 


~dt 


rCKjLjajtjT) dT da 


(4.15) 


The double integral t-enp. in this equation will be recognized 
as the net rate of change of niass of solid particles [KsIj] 
in suspension, vis.. 


dIg(K,l,t) 


/ / K,(a,t,T) r(K,l,a,t,r) dT da (4.16) 


dt 0 0' 

Substitution for rCnjL^ajtsT) from equation (4.4) leads to 

CO t 


dlX(K,l,t) 

b = _iai (K,L,t) / f A(a) N^(a,tj‘^) dX da 

“ 0 0 


dt 


avi (K,L,t) / / 


t iii(a,t,'r) 

0 0 m(a,t,7') 


A(a) dT da ^ S f D(K,L) 


J 

0 0 


ra(K ,1 ,a , t ,T) I,^( a , t ,r ) dT da 


(4.17) 


Substitution of e-juption (4.S) gives 




-Eg(I.,I.,t) izi + D(X,l)lg(K,L,t) 

(4.18) 

where 0 is a weighted mean bubble size defined ao 


dt 


0 = f S A(a) IT (a,t,T) dT da 


0 0 

'■■■ ' OO, ■ 

/ A(a) li^(a) da 
o 


(4. 19) 
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?/here H-.p(a) is the steady state nianber of btibbles of size a 
in the pulp and is a time dependent weighted mean 


bubble size defined as 

c 

= I 


oo 


t in ( a j 1 5 f ) 


j : r — A (a) 

O lii { a 1 1 atf} 


"b 


(a,t,T) da (4.2C) 


By definition 


m( 


a . 5 ? 


:) 


( 4 . 21 ) 


m(a jt jX") 

It then follows that 

'-V (t) < 0 for t > 0 

Combining equations ( 4 . 16 ) and ( 4 . 1S) with (4.15) 

-r j.\ ni 


( 4 . 22 ) 


dt 


Mg(K,l,t) 0 [ 1 - — ] 

-- D(K,L) li-j^(E5i?i) — ^ 


(4.23) 


Since froth remo-al is assumed to be vary fast and mass of the 
solids is conserved 


A.. 

'(it 


r M^(K,l,t) + M, (K,I,t) + M^(K 5 L,t)] = 0 

S ' M . 


(4.24) 


v/here M^(K,l,t) is mass of particles [K,!] in the concentrate. 
Substitution of equations (4.1 8 ) and (4.25) i^ico (4.24) reoultu in 


dlvT(K,L,t) 

0 --, (K,l,t) 

dt ° 


(4.25) 



Equations ( 4 . 18)5 (4.23) and (4.25) then formally describe the 
flotation kinetics. However, even approximate characterization 
of ■'^(t) by experimental means or from theoretical considera- 
tions on the basis of our current understanding of the physical 
mechanisms of the flotation subprocesses, presents formidable 
problems. But it is of passing interest to note that some 
simplification in the model is possible, if it is assumed that 
m(a,t,r) is time invariant, i. e., 

m(a,t 5 T) = m(a) (4.26) 

and both A(a) and £(&) increase v/ith increasing bubble size 
in such a manner that 

i-ial ^ f (4.27) 

m(a) 

where i is some unknoTO parameter independent of bubble size. 
¥/ith these assumptions equation ( 4 . 20) can be revritten as 

M^(t) = i f m(a,t,T) ll^(c’,t,T) dT da 

■ 0 0 

= '5 M. (t) (4-28) 

i D 

Substitution Of equation (4.26) into (4.13) and (4.23) leads to 
follov;ing two equations ^ 
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and 

(lWL(K,L,t) __ < 

-“'-dT = 0 [1- 

- D(K,L) ~ XM^(K,L,t) (4.30) 

and as before 

m (KjLjt) 

^ (4.31) 

The term in sq_uare bracket in equations (4.29) and ( 4 . 30) 
explicitly, but only approximately, accounts for the fact 
that the capacity of bubble phase to accommodate more parti- 
cles is limited b 3 ^ the amount of particles already present in 
its Surface. 

Fre e flotat ion C ase: 

The potential application of equations (4.18), (4-23) 
and ( 4 . 25 ) is in the important free flotation case which is 
frequently encountered in practice'^^’^*'. Free flotation, 
which implies that the bubbles are not over loaded, should 
result when one or more of the following conditions are met. 
The aeration rate t- is very large, \'.hich happens to be 
invari£ 5 bly true, and f(IC,L) is relatively significant in t 
relationship to K. These stipulations concurrently justify 
the quasi-steady state assumption for M^(K,L, t) , that is 

(aivi^ac,L,t) 
dt 


0 as M^(K,I.,t) 0 


(4.32) 
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Moreover, from the definition of free flotation, it immediately 
follows that m/m' - 0= Iherefore, in equation (4.20) 


(t) -> 0 


(4.53) 


Suhstitution of equations (4.32) and (4.33) into (4.23) 


re suit s in 


M^(K,L,t) 


0 K Mg(K,L,t) 

D(K,L) + A 


(4.34) 


Substitution of equations (4.33) and (4.34) into (4*18) and 
(4.25) and rearrangement of te3?ms results in following two. 


equations 


and 


ciig(^L,t) 

”*dt 


dlii^(K,L,t) 

“”"“"dT'"” 


-ICa. Mg(K,l,t) 


= ‘ ICa II (K,L,t) 


where the apparent flotation rate constant is 
Ka 


d(k jli) + ,X 


Integration of equation (4, 35) give; 


(Iv,I,t) = M (K,l) exp f-lCa tj 

S ■ D . . . ■ *« 


where 


= M^(K,I,0) 


(4.55) 


(4.5s) 


( 4 . 37 ) 


(4.38) 


MgdC,!) 
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'Total amount of particulate species of size 1, remaining in 


the pulp -is 


M^(l,t) = f M (K,L) exp [-ila t] dK 


(4.39) 


Mg(l jt) 


f M (Ka,l) exp [-lia t] dKi 


(4.40) 


are a change of scale from K to Ka has been made so that 


, (Ka.l) 


Mg(H(Ka) ,1) 


vfhere l(Ka) is a solution to equation (4.37) i.e, 
K = H(]Ia) 


(4.41) 


(4.42) 


Eqtiation (4.40) is evidentlj/- equivalent to equation (2.17) 
for single size species, and equation (4.37) describes the 
relationship that exists between the apparent flotation rate 
constant Ka and the first order attachment rate constant K. 
This relationship is important not only for the insight it 
provides into the structui'e of apparent rate constant in terms 
of 0,Xj K and D but, in theory, it also permits evaluation 
of within a constant of proportionality from 

Mg(ICa,L) since the latter can be computed from experimental 
dlta. Eurthor re'duction in equation (4.57) is possible when 
free flotation results from stipulation of very large aeration 
rate. Thus, when A » 3)(K,L), a simple, ^j.elationsb^^^xis|s 

i \ 

between ICa. and K \ 

\ . . K •■■"lllJ 


Ka = 0 K 


^cc. 
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This equation should be useful for deteraiining the effect of 
aeration rate on the flotation rate distribution. Moreoverj 
integrating equation (4.40) over entire particle size range 
gives 

DO 

= / MJli,t) dl 

0 ® 

CX) OD 

= / / M^(Ka,L) exp [-Ka t], dKa df 

D 

0 0 
CO 

= / M^(Ka) exp [-Ka t] dKa (4*44) 

s 

0 

Thus 5 we find that in case of dispersed size particulate feed, 
first order flotation kinetics with distributed apparent 
flotation rate constant is valid only if the mean residence 
time of the bubbles is very small. 

BjLPjqyjrj^ anji^. 

P.ecover.y and grade are the two most ioaportanb aspects 
of the flotation process. Although, the analysis presented 
thus far deals with recovery only as given by 

l^(L,t) = / M^(Xa,l,t) dKa (4.45) 

it should be noted that the flotation experiments conducted 
for detemiination of the distribution ivig(Ka jl) in equation 
(4.40) by some numerical technique (a suitable algorithm will 
be presented in the next chapter) can he extended to compute 
the distribution ? (Ka,!) also, where ?g(Ka,3j) dKa is mass of 
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valuable rnineral embedded, in the feed of size L which floats 

with apparent rate constant in range ICa to Ka + dKa. Thus, if 

?i/e conduct flotation experiments and deterinine, in addition 

to the mass of valuable mineral V (L/c) remaining in 

s 

cell at time t, Laplace inversion of 


Vg ( L , t ) = f Yg ( Ka ,1 ) exp [-Ka t| diva. (4.46) 

0 


will give Y^(Ka,ij). The cumulative grade of the material in 
pulp is 


or 


. , V (L,t) 

0 < G (L,t) = ^ 

■ = led-t) ■■ 


GgCL.t) Mg(L,t) = Yg(l,t) 


(4.47) 


(4.48) 


Application 


of Borel ' s theorem 
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for inversion of the Laplace 


transform gives 

Yg(Ka,L) = 


ICa 

/' G„(k 4,L) M rKa~la,L) dsi 

O ' o 


(4.49) 


Gg(lva. jL) will be icecogniized as the impulse response function 
in sample space and it relates the tvro feed distributions 
I»i^(Ka,L) and Y (Ka,L). G„(Ka,L) may be calculated by the 

O . O ^ V 

Ys(L,t) 

Laplace inversion of experimental data or directly 

' ■} R2! 

from equation (4.49) by well Icnown techniques. ^ 
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DISCUSSION MD CONCLUSIONS:- 


It should be stressed that the approach adopted in 
our analysis is not unique. One can formulate alternate 
kinetic model in terras of multivariate distributions which 
includes the grade of the individual entities, as done by 
Zaidenberg et al . Our procedure, however, seems to possess 
inherent advantage of simplicity and practical utility. The 


analysis of flotation kinetics can be carried out by computing 
tv/o parallel mathematical streams, one for total mass recovery 
and other for mass of valuable mineral keeping in mind the fact 
that these two computations refer to the same particulate 


entities whose common attribute is Ka or E but the dependent 
variable is particle mass or its valuable mineral content. In 
this connection, consider a raultivariate distribution as 
follows: 

M^(S,L,ff) dS dL dU (4.50) 

which is the mass of particles with attribute surface affinity 
or reactivity including its hydrophobicity ranging from S to 
S + dS, particle size from L to L + dL and grade (fraction of 
valuable mineral) in range G to & + dG. Hera S is an index 
of particle surface characteristies such as fractional area 
of the valuable mineral on surface, adsorption of flotation 
reagents and geometry of .particles (in particular high energy 



loc.ctions) vmicu, in tircn, may or may not depend 
on L and G. Now the Particle-bubble attachment rate constant, 
Iv, can oe expressed as a function of S, 1 and G 

" K(S,L,G) 

Accually L is perhaps a function of particle mass also which 
IS, however, taken care of bjr inclusion of terns 1 and G. It 
is now possible to convert the density function M^(S,I,G) into 
an aloemate trivariate density function It is 

assimied thar equation ( 4 . 51 ) can be solved for S as a function 

of n,l,G and, one to one transformation is possible. Let the 
solution to equation (4.51) be 

It then follows that a trivariate distribution in attribute 
K, L and G is 






,1 5 >_r 


) [ 


iS.siuD. 

oK 


] dK dL 


dG 


which can be written in simpler notation as 
Li2^k.5L5Gj dlL dL dG 
'X’herefore, 


(4.53) 


(4.54): 
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These two distributions are then the points of departure for 
the kine-cic model described here and result in tractable 
mathematical forms and concurrently provide the desired 
in forma I'^fiuh respect to botn grade and recovery for a 
given particle size. 

The icinetic model presented here is based on rather 
simplified formulations of particle bubble attaclmient detach- 
ment rates in equation (4.4), Detailed investigations of 
particle-bubble aggregation, adherence and detachment sub- 
processes lead to more complicated relationships which, 
nontheless, continue to suffer from various assumptions and 
extreme simplifications. In fact, rigorous derivation of 
these rates from first principles and explicit solution to the 
resulting equations is not yet possible. 

The following results emerge from the kinetic analysis: 

1. The particle-bubble attachment (interception and 
adherence lumped tc,gether)rate constant IC is related to the 
apparent rate constant by equation (4.37). 

2. The a:’ole of bubble population can be incorporated by 
a v/eighted mean 0 as in equations (4.18), and ( 4 . 23 ). 

5. I'irst order flotation kinetics with respect to distri- 
,./ ■ ' state 

buted K is possible if qua si- steady/as simaption in equation 

(4.30) is made. 
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4. J'irst order kinetics for all size particles require 
that the residence time of bubbles be very sEall. 

5. Ti'ie problem of grade can be simplified considerably 
by t?/o parallel analyses in total mass and mass of valuable 
mineral which are related to cumulative grade by equation (4.4?) 




CHAPTER 5 


ESTEiATE OP THE EEOTATIOH RATS DISTRIBUTIONS 

As po 111 bed out earlier also, one of the most important 
problems in the ilotarion kinetics analj^'sis is that of evalua- 
tion of bhe flotation rate constant distribution in equation 
(2. 17) which is reproduced here. 

C50 

Mg(t) = f M (Ka) exp [-Kg. tj dKa (5.l) 

0 

In this chapter a relatively straight forward method for 

reliable computation of rate distribution from M^(t) data. 

s 

which exploits the structural features of the problem, is 
presented. The inversion is carried out in terra of cumulative 
distribution rather than the frequency function. The algorithm 
consists of discretizing the appropriately transformed integral 
in equation (5.l) by the Lobatto quadrature and then 
carrying orit optimization of the resulting constrained variables, 
first, the general method is developed and then its applica- 
tions under different conditions of numerically generated data 
with suitable amount of errors are presented. It may be 
pointed out that the validity of any estimation method can be 
established only if one can extract back the distribution 
that is used to generate the numerical data. The closeness of 
fit between experimental and computed M^Ct) is an unreliable, 
even misleading, criterion for this purpose becaiise the Laplace 



59 


transform is inherently a very stable operator and consequently 

ligCb) is quite insensitive to reasonable variations in M (Ka). 

s 

NUIaEElICAl 1/EE. SION 01 THE LAPLACE TRAEoEOBia: 


After nuj.aerous attempts it was established that success 
is more readily assured if inversion is carried out in the 


cumulativ^e distribution mode. This is because the cumulative 
form oi che distribution function is inherently smoother and 
less sensitive to the influence of errors. Moreover, this 
mode provides a number of constraints wiiich considerably 
shorten the region -of search in the optimization step and 


prevent wild fluctuations vdiich tend to occur in the hitherto 
proposed methods for numerical inversion of the Laplace trans- 
form in presence of errors. Let R(Ka) be the normalized cumu- 
lative distribution function, hence 


R(iia.) = / UCld.) dKi 

Ka ^ 


where 


and 


R(0) 


Mg(Ka) 


S diva 

O 


DR(Ka) 

aKa 


Mg(t) 


f exp [-Ka t] dICa 


(5.2) 


(5.3) 


(5.4) 


Substitution' of equation ( 5.4) into ( 5, 1) leads to 

o-vn +-1 (5.5) 
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Integi:ating by parts we obtain 

oo 

M(t) = J R(Ka) exp [_Ka. t] dKa 


(5.6) 


where 


M(t) 




t 


(5.7) 


Equation (5.6) is the cumulative distribution analog of 

equa-Gion (5.I) and it is required to extract Rdfe) by inversion 

of the laplace transform. By definition of the Laplace trans- 
form 


M (.-i,.) ^ f Qxp f- M—tn 

r\ S 


dKa 


(5.8) 


v/here s is an adjustable scaling factor, (s >0). Define a 
transformation variable x as • 


X 


= 2 exp [- ^ ] - 


(5.9) 


ilhere 


2: ore , 


Ka = s 


C ill- ] 


(5. 10) 


and 


L^ ] dx 


(5.11) 


Substitution of equations (5. IO) and (5.11) into (5,8) 


give s 




dx 


(5.12) 



61 


The limits oi the integral are now suite ole for a n-point 
loba b L-o puadrarure formula for discretization of the integral 
term, Therefore., 


s 


.1. 1 r ., ILll -t iC- ll 1 

„t\L“ r,{a-l) J 

dL ^ 


n-l 

+ 2 w. f(x. )i (5.15) 

i=2 ^ ^ i 


where w^. and x. are i-th weight and abscissa, respectively, 
of the liObatto quadrature formula (both are tabulated in the 
literature and 


^(x. ) = P(x. ) [x.4-l] 


t-1 


(5.14) 


where 


Mien 


and 


P(Xj_) = R(s In [x™j'J) 

1 

Ka = 0 


P( 1) = 


(5.15) 


(5.16) 

(5.17) 


Sirjailarl.y, when Ka = 

x^ = -1 ( 5. 18) 

and 

P(-l) = 0 (5.19) 

On substituting equations (5.17) and (5. 19) , equation (5.13) 
is simplified as 


1 {±) 


n 


r^rrr 


1 


n-1 


2 1=2 


E w- P(x . ) [x. -fll 


t.-1 


(5.2G) 
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Tn equai--'-’ - ■ ( 1:5,20)5 the only unknowns are i = 2 , 3 ,... 

n_1, I3q.ua ■t ion (5.20) cj 


R(Kaj^) and Ka.; a 


an be solved for '£'(x^) , and hence for 


5 using 


some oP'tiKsi- 


function c i io in e 


are related by equation (5.9) 
tion techniques. A least squares objective 

following form v^as found to be quite suitable 


for this purp 


O SP 


E 

□ =1 


t . 


30 


n(n~T)‘ 


n~l t.. 

Z w, f(x )[x.+l] 3- 
i=2 ■' ■ ^ 


(5.21) 

where li is bhe number of data points5Mg(t . ) 5 j = 132,..., H. 

D 

I'he problem is now reduced to minimization of the v. pective 
function wibh sippropriate choice of l'(x^). 

The unlcnovm variables of the objective function in 
equation (5.21) are subjected to a number of constraints which 
arise from t-he following properties of the cumulative distribu- 
tion function i 


(i) !l(lva.; ) _> O3 0 C 5. 


( 5 . 22 ) 


(ii) R(la) is a iioiiotonioally non-increasing function of 


Ka 3 i.8. • 

■1 y rX his- ) X 3.2) — — y rCkp o 

when 

0 < Ka < Ife-g * • • S:-^m ’ * * - ” 

(iii) R(o) = 1 and R(“) = 0 


(5.23) 


(5.24) 
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ihe Loba cto quadrature formula was chosen because it is 
exact at the tw extremes and, therefore , the last constraint 
is automatically satisfied (equations (5.17) and (5.19 )). 
fhe second constraint, which is key to the success of the 


meGhoo. , drastically narrows doTO the region of search to that 
withiii Ghe conj-ine of two neighbouring points and prevents ?/ild 
oscillations. -This feature can be readily incorporated in many 
optimization algorithms such as iinivariant search method, 
similcir to I'riedman and Savage"^^"^ method using Fibonacci 
search technique ^ Rosenbrock''^-% Random Search 
methods etc. The optimization techniques used in this, as well 
as in folio w'ing chapters are discussed in Appendix B. Essential"! 


same results were obtained v/hen the i.nitial guess was changed 
or different optimization methods were employed. It v/as conclu 
ded that this step of the inversion algorithm does not present 
any special difficulties. 


Eight point Lobatto quadrature formula was employed to 
approximate the integral in equation (5.12) , the values of 
weights and abscissas are given in Table 5- 1- In some instances 
this did not cover adequately the abscissa range of interest 
especially when function R(Ka) varied sharply over a narrow 
interval of Ka. This difficulty was readily taken care of by 
Suitable adjustment in the scaling factor s in equation (5.9) 
in order to generate more points in the region of interest. 

The optimization was usually initiated with s . equal to unity, 
and subsequently the scaling factor was changed if found necess^- 



T'able 5. 

l: Abscissas and 

Weights of 

8 Point Lobatto 


and Legendre 

Quadrature I’ormulae. 

Lobafto 

Quadrature 

Legendre Quadrature 

Abscissa 

Weight 

Abscissa 

Weight 

1.00000000 

0.03571428 

0.96028936 

0. 10122854 

0.87174015 

0.21070422 

0.79666678 

0.22238103 

0.59170018 

0. 34112270 

0.52553241 

0.31370666 

0.20929922 

0.41245880 

0. 18343464 

0.36268378 

-0.20929922 

0.41245880 

-0. 18343464 

0. 36268378 

-0.59170018 

0. 34112270 

-0.52553241 

0.31370666 

-0.87174015 

0.21070422 

-0.79666678 

0.22238103 

-1.00000000 

0.03571428 

-0.96028986 

0. 10122854 
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EYALITATION OE THE MTEOD: 

1 * proposed ine 1 ;hod was tes'ted with 

numerical data Mg(t) generated from the follomng test dis- 
trihut ions. 

(i) triangular distribution 

(ii) Rectangular distribution 

(iii) Exponential distribution 

(iv) G-amma distribution 

(vj Bimodal (mixture of two) gamma distribution. 

table 5*2 gives the assumed values of the parameters 

in these distributions and the explicit forms of function M (t) 

s 

obtained by substitution of these distributions in equation 
(5. 1 ) 5 'Which were then used to generate essentially error 
free numerical data divided in 20 points at well spaced 
intervals in range 0.01 < M^(t) < 1 . figures 5.1 and 5-2 show 
that the cumulative distribution estimated by the proposed 
inversion technique is in reasonably good agreement with the 
test distribution in every case. Very accurate inversion of 
triangular and rectangular distribution is not feasible because 
of abrupt changes in the slopes of these frequency fimctions. 
Intuitively distributions with such sharp variations are not 
expected to be encountered in case of actual minerals. 

Eig. 5.2 also illustrates the use of sealing factor for 
generation of additional points on the curve. 



Table 5,2: functional Form of M^dCa), Mg(t) and the Values of the 

Parameters Used for Synthetic Generation of Flotation Data 
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2 # g|‘£g,g,'fc,_qf_.,S^ro rs in thft Data; 

its mentioned earlier, inversion of the Laplace transform 
is, very sensitive to the accuracy of the data., fhe effect of 
errors nes been scudied for exponential and gamma distributions 
by rounding off Mg(t) data (O < ^gCtj < 1 ) beyond 1,2 and 5 
decimal places, using 16 data points at nvell spaced time 
intervals. Percentage errors, thus introduced in the M (t) 
data is shown in Pig. 5*3 for the gamraa distribution for all 
three sets of data. It will -be seen that with increasing flota- 
tion time the data has become progressively degraded. This 
trend is expected to be in line with actual experimental results 
as the rate of flotation decreases continuously with time 
and errors accumulate. Pig. 5.4 shows that, as anticipated, 
the deviation of calculated distribution from the corresponding 
theoretical curve increases with increasing error content in 
Mg(t) but the method is reasonably successful if the data is 
accurate to twro decimal places. It is felt that this degree 
of accuracy is attainable in a reasonably carefully conducted 
experiment. Purther discussion of the Laplace inversion v/ill 
be illustrated with two decimal place data only, generated from 
the gamma distribution. A plot of log (M^(t)) versus t for 
this distribution is shovm in Pig. 3.5. Pig- 5.4 also shows 
that inversion algorithm is robust in that the test data which 
were read off directly from the graph in Pig. 5.5 co.uld also be 
inverted with acceptable precision* The data points shown in 
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Fig. 5.5 represent M^Ct). values back calculated by substi- 
tuting the computed gamma function in Fig. 5.4 in equation (5-20) 
It will be seen tnat excellent agreement is obtained. ' 

3 . 


It IS obvious that the computational efficiency of 
■the method in optimization step will depend on the choice of the 
initial guess vector. Fig. 5.6 shows that the widely different 
starting guesses essentially converged to the optimum value and, 
therefore, the success of the method does not seem to depend 
on the choice of the base point. But the search time was found 
to be strongly correlated with the starting vector. 

Ml e.cJL„Pl. Jnmber of lata F eints : 

'i‘he effect of nionber of data points used in inverting 
the Laplace transform is shOTOi in Fig. 5.7 where the distribu- 
tions were estimated using 75,11 and 20 well spaced data points 
in range of O.O1 < M^Ct) < 1. it will be seen that the results 
are quite satisfactory even when only 7 data points were 
employed. Fowever, it should be emphasized that if relatively 
small n-umber of data points are used then it is necessary to 

ensure that the data should be v/ell distributed over the entire 
range. . . 

^ p>f.. Non-Floa tqble 3 : ' , 

I'he synthetic data employed, thus far, was restricted 
to floatable species only. In practice one frequently 






Compute 
Data poi 






77 


encounters ores with significant amount of non- floatable 
component. To examine the applicability of the inversion 
method for this case, s:inthetlc data tath different percentages 
of non-floetables r«re generated, let Mj(t) bo the fraction 
of floatable material remaining in the cell at thvie t and 1^, be 
the fraction of non- floatable in the original feed. The frac- 
tion Of total material remaining in the cell at time t is given 


Where 


and 


Mg(t) - M^(t) [l-4ij + 


M.(0) 

O 


(5.25) 


( 5 . 26 ) 


Mf(0) =1^(0) [1-Mj (5^27) ■ 

^ig(t) data were generated with 10 to 50 percent non 
iloatables. Substitution for M^(t) in terms of the frequency 
function Mg(Ka) in equation (5.25) leads to 

oo.' ■ ■ 

ni^(tf = / M^(lCa.) exp [-Ka t "] dKa + Ir, (5.28) 

O c. 

Substitution of equation (5.4) into (5.28) gi¥es 

^ ' ^ exp [-Ka t] dlia + (5.29:) 

Integrating by parts 

'DO ■ " ■ ■ 

M(t) = [i-MJ / R(I{a) exp [-Ka. t] dKa (530.) 


or 


M(t) = / R'(Ka) exp [~Ka t| dKa 

lot,: I....' M: • • 


(5.3f) 
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where * 

R'(Ka) - [i-Mj R(Ka) ( 5 . 3 . 

Although equation (5.51) is similar to equation (5,6), it 
should he noted that now R'(0) is equal to an untoiov/n [1-MJ . 
IhuSj one of the end points corresponding to = 'j iu equa- 
tion ( 5 . 17 ) is no longer preassigneds and therefore, the 
Lohat to quadrature cannot he used. Using the transformation 
in equation (5.9), equation (5.51) is transfomred into 


M (~y~) 


t N o 1 


R'(s In dx (5.35) 


2 " -1 

and. disci'ebized using a 8 point Legendre quadrature formula 
Now the required objective function corresponding to equation 
( 5 . 21 ) become s 


E -i 

3 =n 


8 


s 


I I 


2 3 i=1 

I (R^, Rg) 

subject to constraints 




. 1-1 


where 


0 < < r-g ... < Pg < 1 


s'(x.) = r’(s In 


(5.34) 


(5.35) 


= P, 


(5.36) 


and and x^ are i-tb weight and abscissa, respectively, of 
the 8 point Legendre quadrature formula tabulated in Table 5.1. 
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OptiiniZdtion results are showr in Pig. 5.8 ard it will be seen 
that for four cases with = 0.1, 0,2, O .3 and 0.5, the calcu- 
lated distribut xons ars quite close to the theoretical curves. 

t , 

oince R (O) is equal to , this method automatically gives 

the amount of noh-floatables present in the feed. It should be 
meritioned fcna l- it is also XDossible to employ Radau' s quadrature 
formula for discretizing the integral v/hich is exact at one 
of the end points, namely, ?Aen Ka - R'(Ka) =0. 

BISCUSSIOI MD COl'TClU SIGNS: 

The proposed method of estimating the flotation rate 
■distribution by inversion of the laplace transform has been 
tested for its reliability using synthetic data based on a 
number of test distribution functions. It is noteworthy that 
the method has been shown to actually extract, back the distri- 
bution function from noisy data, provided the errors are random 
and unbiased; which has not been demonstrated for earlier methods 
using graphical or numerical techniques. Moreover, in our 
approach no a priori assum.ption is made -regarding the form of 
the distribution function. This' has ir,iportant consecuence as 
sho'wi in Chapter 7. The inversion algorithm is reliable and 
robust as it can tolerate considerable error in the data which 
is inevitable in any experimental study. The non-f loatables can 
be estimated also. Estimationof the cumulative distribution 
rather than the frequency function does not present any diffi- 
culty in the sense that the kinetic equations can be reformulated 
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completely m the cumulative mode. With the availability 
of the algorithm it is now possilDle to analyse the structure of 
the apparent rate const.-nt in temns of particle size, aeration 
rate, pulp densioy, flow rate etc. Thus a broad new area for 
flotation iiivescigation becomes accessible for systematic 
investigations 3 some aspects of viiich are dealt in next tvro 
chapters. 



CHAPTER 6 


the EEESCTS 0? A3JI10H RATE, EARTIOIE SIZE MI, EDIT 
On TaE i'EOTATIOH li/iTS EiSTSIE’JIIOOIS 


As insiiliioiisd. 6f: T’liB'" „-p -3. 

s-cudies in past on 

the effeets of proceee variables on flotation klnetloa have 
been canniei out in tei.3 of none ovenofe or instantaneous 
flotation rates oi pararaeters in eisplrical klnetio expressions 
and not in tems of rate constant distribution as a -rhole. The 
reason for this oan be attributed to the fact that the notion 
of the rate distributed particulate species is coaparatiuely 
new and, moreover, till now no satisfactory and reliable 
method was available for estimation of the apparent flotation 
rate distribution Mg(Ka) in equation (5.1), given noisy 
experimental data M^(t). In this chapter, using some of the 
nuaaerous experimental data available in the literature, the 
changes that occur in the flotation rate distribution, as a 
whole, as a function of aeration rate, particle sise and 

pulp density have been studied. The procedure comprises of 

the -following two steps. 


J'lgCt) data is employed for computation of the flotation 
rate distribution by numerical inversion of the iaplace trans- 
fonn in equation (5.1) using the algorithm developed in the 
preceding chapter. The computed 'cumulative apparent flotation 
rate distribution, R(K;a), is then examined for appropriate 
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relationships hetvi/een R(lCa) and the aeration rate, particle 
sine and pulp density. In this manner it has been possible to 
develop useful correlations for the structure of the apparent 
flotation rate constant 'Ka in the distributed parameter 
flotation model and to check the validity of equation (4*43) 
for the effect of the aeration rate in the phenomenological 
model developed in Chapter 4. 

EFI'BOI Of AERATION RATE; 

Eig. 6. 1 shov/s normalized cumulative distributions in 

the apparent flotation rate constant Ka for four aeration 

118 

rates. These were computed from Tewari and Biswas data 
for four different aeration rates for -140 + 200 mesh (A. S.T.M.) 
calcite by numerical inversion of the Laplace transform. The . 
data employed for computation of these distributions was read 
off from the graph and is tabulated in Table 6.1. It v/ill be 
seen that distributions shift to the right and the flotability 
increases rapidly with increasing aeration rate. 

Under the condition /- >> D(K,L) , the relationship between 
apparent flotation rate constant, Ea, and the first order 
particle to bubble attachment specific rate constant K is 
given by equation (4«43). fhe parameter 0 in this equation, 
an unspecified function of weighted mean bubble size and 
hydrodynamic conditions prevailing in the cell, depends pri- 
marily on the rate of aeration. The distributions in iig. 6, lean 
now be used to test the validity of this simple relationship. 








85 


I'aole 6.1: fraction of Calcite Remaining in the 


Cell 

Rates 



Cumulative Times 


for Jour Aeration 


Time 
(Sec. ) 

Aeration Rate in litre/rain 


2.4 

2. 2 

1.9 

1.5 

0 

1 

1 

1 

1 

6 

■ 0.805 

0.942 

0.980 

0.993 

15 

0.422 

0.824 

0.945 

0,97s 

30 

0.252 

0.603 

0.376 

0.953 

45 

0.181 

0.450 

0,795 

0.916 

60 

0.141 

0.352 

0.705 

0.852 

75 

0. 107 

0.288 

0.624 

0.778 

90 

0 . 0 S 1 

0. 248 

0.544 

0.705 

■ 105 

0.053 

0.212 

0.485 

0. 64s 

120 

0.040 

0, 185 

0. 440 

0.606 

135 

0.027 

0. 166 

0.410 

0.562 

150 


0. 146 

0.384 

0.530 : 

165 

- . 

0.133 

0.367 

0.500 

180 



0.119 

0.350 

0.478 
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By dexinicion 5 the ciHiiulative rate distribution in 
attachment rate constant K can be written as 


R(K) = r M (K') dK’ 


K 


Substituting for K from equation (4.43) 


S(K) 


I? ( iSS 
^ ^ 0 ' 

i i ) « 


iva 


Equation (4# 39) can be rewritten as 

CO 

= / M^(K) exp [-Ka t] dK 


( 6 . 1 ) 


( 6 . 2 ) 


(6.3) 


Comparison of this equation with (2.17) leads to following 
equation 


,(K) dK 


M^(Ka) dICa. 


Substitiztion of equation (4.43) results i 


in 


(6.4) 


Mg(Ka) = 

Gombining equations (4. 43) and ( 6. 5) with ( 6, 2) 


R (^) 


/ M„(Ka> dKa 
Ka ® 


Erom equations (5.2), (6.2) and (6,6) it folj'oY/s tha' 
R(Ka) = R (^) = 1{K) 


( 6 . 5 ) 


( 6 . 6 ) 


(6,7) 
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Since R(iv) is the cianulative rate distribution in Kj it is 
independent of the aeration rate. In general., the valiie of 
the function 0 is not known, hut it can be readily estinated 
v/ithin an uncertaiirne of a constant of proportionality as 
follows. Denote the median apparent flotation rate in distri- 
bution R(ua) by Ila(C.5). Rx’oiri equation (S.7) 

0.5 = R(ka(0.5)) = R 

= R(I[(0.5)) (6.S) 

Hence 

MLOji-il. _ (6.9) 


But IC(0.5) is a constant for a given assembly of conditioned 
particles. Hence, 0 is directly proportional to lLa(0.5). 
Substitution of equation (6.9) in (6,7) gives 




( 6 . 10 ) 


■The constant H(0.5) may be absorbed in function R to give E 
and equation (6,10) becomes 



( 6 . 11 ) 


These relationships along with equation (6.9) imply that for 
every aeration rate there exists a parameter 0, which is 
proportional to Ka(0.5), such that if the abscissa scale in 
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plots of R(lLa) versus Ka is divided by the corresponding 
value of ICa(0.5) the resulting curves for different aeration 
rates should collapse onto a single similarity distribution 
R equal to R(I<.) which, in turn, is independent of the aera- 
tion rate. Rig. 6.2 sho’vs that within limitations of noise in 
expez'iiiiental data and computational error in the inversion of 
the Laplace transform, the similarity distributions R* as a 
function of dimensionless similarity variable Ka/Ka(0.5) have 
indeed collapsed even though' Ka( 0 . 5 ) varied nearly 18 times from 
0.14 to 2.5 min"”*. It may be concluded that the simple rela- 
tionship given by equation ( 4 . 43 ) is apparently valid. Rig. 6. 3 
illustrates the rapidly increasing strong dependence of 
Ka(0.5)s (i.e. of 0 ) on the aeration rate. Thus, the rate of 
aeration can be an important control variable in operation of 
flotation circuit - 9 '- However, in some flotation system^ 

such strong correlation between flotation kinetics and aeration 
rate has not been found, from the data on bubble size distributio" 
in pulp it is also possible to extract some approximate informa- 
tion on the nature of the function of 0 since in view of equa- 
tion ( 6 , 9 ) s there might exist a moment of the bubble size 
distribution which may be proportional to Ka(0. 5). 



Aeration Rate, Ka(0.5) 

lit. /hr min"^ 




cumulative flotation rate 
as^-;tunctior> ot dimensionless f!otat:o 





Median app^ rate, which is proportional to 

^ function of aeration rate. 



EI'E'SC'T Oi’ PARTICLE SIZE: 


Riga, b.45 5 . 5 , 6.6 and 6.7 sho’vv the apparent flotation 
ra L-e di&tributiona computed from Tev/ari and Bisv/as data for 
four size fractions of caleite ^ '*'^5 Tomlinson and Eleming data 
for five size frac cions oi apatite^^ , and Volin and Swami data 
for five suze fractions of hematite and for six size fractions 
of magnetite particles"'^;, respectively. The data employed for 
computation c± tnese rate constant distributions is reported 
in Tables 6.2, 6.5? 6.4 and 6 . 5 . In first three cases the 


originally reported data was in graphical form v/hereas in the 
latter tvra cases it v/as in tabulated form. In some cases of 
tabulated data where number of data points was not sufficient 
graphical interpolation v/as carried out. Interpolated points 
have been tabulated along v/ith the reported data points in 
Table 6.4 and 6 . 5 , We shall now investigate if the approach 
taken in the previous section can be extended for analysis of 
the particle size effect also. 

It is assumed that the contribution of the particle 
size in particle to bubble attachment rate constant K may 
be separated in the following manner 

K = r.ih)K° , (6.12) 

where 'l(L) is a function of particle sise only and K°, the 
reduced rate constant? depends on grade and surface hydro- 
phobicity which is assumed to be independent of particle size. 
Combining equation (6.12) with (4.43) 
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Ea = 0 %L)t 


rO 


(6.13) 


iilieia the aeration rate is held constant this equation may he 
revn^itten for convenieiice as 


Ka = n(L) 1 




(6.14) 


and 


1(1) = 0 >Kl) 


(6.15) 

jiollowixia the argumentsj used in derivation of equations (6.1) 
to (6.11)5 it is possible to derive a relationship between 
cumulative rate distribution in terms of Ea and cumulative 
ra ls distribution in E*^, R°(K°). Again by definition 


0 


R"'(E°) = / M° (E°) dK' 

K° 


where M° is the frequency function in E°, 
equation (6.14) into (6.16) gives 


( 6 . 16 ) 

Substitution of 


O/'KaA r 1 t„tO / Ea 


Let 


R^(^^-) 


Mg(Ka) 


/ 

lia 


^ ) dEa 
1 ^ r|' 




Substitution of equation (6, IS) into (6.17) leads to 




o /ICaN ^ j K (£a) dEa 




= R(Ka) 


( 6.17 


( 6 . 18 ) 


(6.19) 
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Table 6.2 Praction of Oalclte aoBiainlng iii the Cell at 
G-umulative Times for Four Particle Size 
Fractions'*'^®. 


Time 

(Sec.) 


Particle Sise in Mesh 


-2704-525 

-200+270 

-70+100 

-50+70 

0 

1 

1 

1 

1 

6 

0.722 

0.775 

0.829 

0.927 

12 

0.451 

0.564 

0.658 

0.855 

18 

o 

• 

ro 

0.388 

^ 0.447 

0.756 

24 

0.288 

0.338 

0.356 

0.684 

30 

0.252 

0,293 

0.290 

0.619 

45 

0.180 

0.226 

0.211 

0.514 

60 

0.144 

0. ISO 

C.163 

0. 441 

75 

0. 126 

0, 148 

0.132 

0.395 

90 ■ 

0.117 

G. 141 

0.118 

0.353 

105 

0. 108 

0.134 

0. 108 

0.537 

120 

0.094 

0. 127 

O 

m 

a 

'O-D . 

0. 314 

135 

0.036 

.0. 124 i’ 0 

■ 0,082 ; ; 

0.503 

150 

0.076 

0.113 

0.066: • \ : 

, ' 0, 283 
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Table 

6 * • 

Praction of Apatite I 

Remaining in 

the Cell at 



Cumulative 

Times for 

Piva Particle Size 



Practions^' 

'i 

r 

• 





Particle size 

iu microns 



— 





(Sec. ) 

134 

84 

67 

51 

31 

0 

1 

1 

1 

1 

1 

10 

0.853 

0.910 

0.945 

0.967 

0.980 

15 

0.780 

0.865 

0.915 

0.950 


20 

0.680 

0.804 

0.880 

0.925 

0# Sb2 

30 

0.430 

0.680 

0.810 

0.383 

0.940 

40 

0.350 

0.595 

0.745 

0,845 

0.925 

50 

0.265 

0.515 

0.690 

0.805 ' 

0.910 

60 

0.220 

0.450 

0.630 

0.775 

0.890 

70 

0. 180 

0. 380 

0.580 

0.740 

0.S77 

80 

0. 155 

0.335 

0.530 

0.710 

0.865 

90 

0.090 

0.300 

0.500 

0.680 

0.850 

100 

0.080 

0-265 

0.470 

0. 650 

0.835 

110 

0.070 

0.235 

0.440 

0. 627 

0.820 . 

120 

0.050 

0.220 

0. 413 

0.605 

0 . SCO 

130 


0. 190 

0.385 

0. 58a 

C.7S5 

140 


0.175 

0. 370 

0.563 

0.775 

150 


0. 160 

0. 350 

0.545 

0.760 

ISO 


0. 140 

0. 330 

0.530 

0.750 

170 


0. 125 

0.315 

0.515. 

0.735 

1J» 


0. 115 

0.300 

0.500 

0.725 

190 


0. 100 

0.290 

0.435 

0,720 

200 


0.085 

0.275 

0.470 

0.705 

210 


■ mm, 

0,260 

0.460 

0.690 
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Tsbl© o»4* Px^cGio];! of Hgiqb L xfo Roriiciiii ing in 
0x111111131; iv 0 I'iniss fop 3?ivG pQpficlG 
Fractions^^l 


Tirne 


Particle 

size in 

E3 Sll 

(Sec.) 

35-48 

65-80 

80-100 

100-115 

0 

1 

1 

i 

1 

5 

0.901 

0.365 

0.055 

0.901 

10 

0.825 

0.738 

, 0.722 

0.805 

15 

0.756 

0.606 

0.605 

0.732 

20 

0. 670 

0.528 

0.505 

0.63s 

25 

- 

- 

0.337 

- 

30 

0.546 

0.366 

0.2S8 

0.492 

40 

- 

0.243 

0.170 

_ 

45 

0 . 565 

0. 185 

0. 153 

0.312 

60 

0.244 

0.096 

0.057 

■ 0. 146 

75 

0 . 1 35 

0.054 

0.024 

0.079 

90 

0.091 

0.031 

0.013 ■ 

0.042 

105 

0.065 

0.022 

-- 

0.022 

120 

rTi 

0 

• 

0 



0.013 

135 

0.038 



0.010 

150 

0.023 



0.007 


uho Gsxl 3 "b 
Size 


115-150 

1 

0.925 

0.850 

0.779 

0.705 

0.593 

0.402 
0.263 
0. 163 
0, 100 
0,051 
0.023 
0.020 
0.014 



100 



e 6.5*. 

Fraction 

at G-umula 

Fractions 

of iiagnetite Remaining in the Cell 

tive limes for Six Particle Size 
,30^ 

Time 
( Sec. ) 



Particle size in Biesh 


55- 4 S 

48—60 

60-65 

80-100 

100--1 15 

150-170 

U 

1 

1 

1 

1 

1 

1 

5 

0.960 

0.928 

0.880 

0.864 

0.900 

0.920 

10 

- 

0.850 

0.775 

0.740 

0.790 

0.822 

15 

0.880 

0.771 

0, 654 

0.61S 

0.673 

0.730 

20 

- 

- 

0.’580 

0.520 

0.562 

0.638 

25 

- 

0. 630 

0.500 

0.431 

0.465 

0.526 

50 

0.782 

0.551 

0.424 

0. 300 

0.372 

0.466 

40 

- 

- 

0. 279 

0.170 

0.235 

0.342 

45 

0. 666 

0.330 

0.212 

0. 107 

0. 175 

0.293 

60 

0.542 

0.159 

0.081 

0.025 

0.094 

0.204 

75 

0.430 

0.040 

0.012 : 

- 

0,041 

0. 143 

90 

0.309 

0.011 

- 

_ 

0.014 

0.091 

105 

0. 168 

o.ooe 

_ ■ 

_ 

_ 

0.056 

120 

0.080 

wrw, 



■ 

0.018 
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From equations (6.16) amd (6.19) It follows that 
R(Ea) = = hO(kO) ■ 

' 'i 

i^'or median apparent flotation rate constant, equation (6.20) 
can be written as 


0.5 


Hence 


= R(Ka(0.5)) 




= K°(0.5) 


( 6 . 22 ) 


Since If ( 0 . 5 ) is independent of particle size for a given 
assembly of conditioned particles, it follo?/s then that ^ 
is directly proportional to Ka(0.5). Substitution in equation 
( 6, 20) leads to 

R(ICa) = R°(K°) ( 6 . 23 ) 

Absorbing constant IC°(0.5) in and denoting the resulting 
distribution by 

R(lia) = = R°(K°) (6.24) 

Accordingly, plots of R(ka) as a function of ICa/Ka(0. 5) should 
collapse onto a single similarity distribution R^, equal to 
R°(K°), which is independent of the particle size. Figs. 6.8, 
6.9, 6.10 and 6.11 show that the distributions in Figs. 6.4, 
6 . 556.6 and 6 . 7 , respectively have indeed collapsed, within 
error limits, and the validity of equation (6*14) is thus 




- S Ka/Ka(0.5) 

f<^rn of cumulative flotation rate distribution for 

^ 1 fitter t ion of ^^Imertsiontess ftotation rate constant Ka/Ko(0.5 
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conf iicjne d. -^ig. 6, 12 shows the plots 
proportional to , as a function of 
The curves for Volin and Swami^^ data 


of ICa( 0 , 5)5 which is 
the mean particle size, 
for hematite and magne- 


tite exhibit characteristic peak in the intermediate particle 
size, in conformity with the well known experimental fact 
that too large or small particles donot float well. The curves 
corresponding to Tewari and Siswas data for calcite, and 
Tomlinson and Fleming data for apatite, how^ever, donot exhibit 
such maximuin flotation rate at some intermediate particle size, 
ihis may be attributed to the fact that perhaps the range of 
particle size examined was not sufficiently extended in these 
experiments. As pointed out earlier the particle size at which 
the flotation rate is maximum would strongly depend on the 
mineral species, its density and surface characteristics, as 
well as, the hydrodynamic conditions prevailing in the cell. 

Thus, it is quits likely that for Tewari and Biswas data flota- 
tion rate will be maximnm at a Particle size range < 50/^ whereas 
for Tomlinson and Fleming data it is maximum at a particle size 
range > 134 a . 

EFFECT OF PULP LSTSITY: 

Flotation kinetic data of Imaizumi and Inoue^^ for 
quartz particles in size range 74—105 microns with five 
initial pulp densities, namely, 7.7,, 23 , 69., 208 and 460 gm/litre 
was employed for computation of a.(Ka). The data which was read 
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off from the graphs is tabulated in fable. 6,6. fig. S. 13 
shows that for three dilute pulps, the apparent flotation 
rate aistributions are, within limits of errors, bunched 
toge-cher and independent of the pulp densities. Li case of 
the t.nK. concentrated pulps the free .flotation hypothesis, 
whicn lorms the basis for derivation of equation (2.17), is 
not likely to be satisfied. The overall rate of flotation now 


depends priiriarily on the capacity of the carrier, air phase to 
accommodate and lift the particles from pulp into froth. 
Consequently , the exact nature of the flotation rate distribution 


has only a secondary influence on the flotation rate and this 
is pernaps reflected qualitatively in the steep slopes of these 
two distributions in fig. 6.13. The pulp density limit for 
free flotation v/ould depend, apart from the strong influence of 
the rate of aeration, on the volume of solid phase rather than 
its mass and it would differ for different minerals. Loveday®"^ 
had shov.n that upto 45.7 gm/litre for pyrite-silica mixture 
the instantaneous average flotation rate is independent of 
initial pulp density, for the first time, however, fig. 6. 13 
exhibits this invariance in direct teims of the calculated 
flotation rate distribution. 

Discussion ill D CON elusions; 

flotation data available in literatiore has been used 
to compute the apparent flotation rate distributions, these 
distributions and corresponding apparent flotation rate 
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Table 6 . 6 : Fraction of Quartz Remaining in Cell at Cumulative 



'■ 2 imes 

for Five 

Initial 

Pulp Densities®^. 

Time 


Pulp density in gm/litre 


(Min. ) 

7.7 

23 

69 

20S 

460 

0 

1 

1 

1 

1 

1 

0.5 

0 . 678 

0.654 

0.717 

0 . 850 

0.967 

1.0 

0.452 

0.371 

0.494 

0 . 666 

0.889 

1.5 

0.262 

0.262 

0 . 294 

0.518 

0.833 

2.0 

0.202 

0.153 

0.20 

0.352 

0.750 

2.5 

0. 155 

0.127 

0.106 

0 . 218 

0.667 

3.0 

0.102 

0.088 

0.073 

0.166 

0.589 

3.5 

0.077 

0.065 

0.037 

0.119 

0.500 

4-0 

0.059 

0.050 

0.027 

0.100 

0.433 

4.5 

0.046 

0.035 

0.021 

0.066 

0.367 

5.0 

0.034 

0.026 

0 . 01 S 

0.055 

0.317 

5.5 

0.027 

0.019 

0.016 

0.041 

0.283 

6.0 

0.023 

0.016 

- 

0.030 

0.250 

7.0 


~ . 

- 

. - 

0 . 167 

8,0 

•mm ■ ' 



0.012 

0.117 



Pulp Density , 
gm / lit 



i Apparent Flotation Rate Constant, Ka 

flotation rate distributions for quartz 
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constcuaus have been analyzed as a function of aeration rate, 
paruicle size and initial pulp density, it will be noted 
that the assumed structure of Ka in equation (6.14) is 
somev/hat similar to the model proposed by Woodbuxn, King and 
Colborn . In present notation their expression for overall 
rabe constant, gi\ 7 en by equation ( 2 . 2 ?) can be rewritten as 

Ka = 0(0, 1 ) K° (6.25) 

where 0 is a specific function of particle size and a hydro- 
dynamic parameter o defined in teims of bubble size and rela- 
tive particle bubble velocity. Kor a given cell, in first 
instance, both £ and e are unloioym functions of the rate of 
aeration. Although equation ( 6 . 25 ) bears formal resemblence 
to equation (6.14), the function 0 suggested by Woodburn et al 
does not permit separation of particle size effect as in equatio 
( 6 . 12 ). 1 'i/hen only one variable is changed, either aeration 
rate or particle size, both models require self similar distri- 
bution in scaled apparent flotation rates as shown in Pigs. 6 . 2 , 
6.85 6.9, 6.10 and 6.11. Unfortunately, due to Inclc of suitable 
lata it has not been possible to discriminate betvreeu these two 
structures of the apparent flotation rate constant. 

■ ' ■'Q'Z ■ 

libodburn et al ^ have shown that for modelling of conti- 
nuous flotation cells for design and control of industrial cir- 
cuits it is necessary to incorporate explicitly the dependence 
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of appare-n-G flo-tation rate constant on particle size since the 
residence tine distribution depends on size. In the present 
approach the distribution as a whole is scrutinized for the 


effects of both aeration rate and particle size. The princi- 
pal conclusion is that the nature of the function R(Ea) does 


no G change wich these ^ariaoles, only the argument of the 
function 5 Ka^ is altered or rescaled as in eq_\^tions (4.43) 
and (6.14)5 and the scaling parameters may be determined from 
the appropriate data. This is a significant result in so far 
as it greatly simplifies the description of the distributed par- 
ticulate species, specially Mien incorporated into a model 
of the flotation circuit. Determination of the apparent flota- 
tion rate distributions and the scaling, parameter for particle 
size effect from continuous steady state data, which is more 
meaningful from the practical point of view has been carried 
out in a similar manner, details of vrtiich are given in the 
following chapter. 

Pig. 6.13 shows that the free flotation hypothesis for 
dilute pulps is justified in context of distributed apparent 
flotation rate constant and, using the inversion algorithm, it 
is possi'.ble to determine the limiting vconcentrat ion of solids 
for free flotation kinetic.? from suitable experimental data. 



CHAPTER 7 


THS EJjPEGTo OP PARTICIiS SIZE MI) PEED RATE OI'I TIHE PIOTATION 
RATE DIS2RIBUTI0IS Bf A COKl'IKUOUS CELL 


Analysis of flotation Icinetic data, for a semi-batch 
flotation oper.ation has been discussed in the previous chapters. 
The approach applied therein can be easily extended for the 
analysis of kinetic data for continuouslj/ operated cell. In a 
continuous steady state operation the mass flov/ rate of solid 
particulate species in tailings and concentrate streamss as a 
function of the residence time distribution, is given by 
equations (2.21) and (2.22). Based on experimental evidence 
Woodburn et al proposed that the residence time distribution 
of solid particles in the pulp can be assumed to be a gamma 
distribution (equation (2.23)) vdth mean and variance of the 
distribution, as given by equations (2.24)-j and (2.23), respec- 
tively. Therefore, equations (2.21) and (2.22) can be 
rewritten as 

■,K+1 


and 




M„(L) 




( Ka 'f p) 


a+1-^ 


( 7 . 1 ) 


/ MT,(Ea,I) [l- 




a+1 


(Ka+p) 


a+1 


] dKa 


(7.2) 


In a 25 Gu ft continuously operated Denver Cell, using an 
impulse of irradiated apatite rich ts^cer, but excluding the 
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flObaoion reagenta, woodburn et al®^ have computed the parame- 
tera of the tesiaence time distrihutiou for a number of particle 
size groups and sIuxtj^ flow rates. The impulse tests were also 
conducted uuuer xlotation conditions to generate data suitable 
for computation of Based on earlier v/orks, Woodburn 

et al initially/ assumed that Mp(Ka5l) is given by the gamma 
distribution 5 but found that 5 contrary to expectations, the 
parameters of the distribution for a given size group were not 
invariant with flow rate. 

In this chapter, it is shown that this dependence on flow/ 
rate is perhaps more apparent than real and seemingly arises 
from a restricted choice of the gamma function for characteriza- 
tion of the flotation rate distribution. First a reliable method 
for computation of IvLp(Ka. ,1) in the steady state expression for 
recovery in equation (7.2) is presented, which is not restricted 
to any a priori assumption concerning the mathematical form of 
the distribution function. Using tabulated data of ’ffiDodburn 
et al it will be then possible to demonstrate that the computed 
distribution Mj,(Ea,L) is independent of throughput , and the 
apparent flotation rate constant, Ka, is related to particle 
size L by equation (6,14)» Ihe estimated flotation rate 
distribution is then employed to predict fairly accurately the 
flotation kinetics of a continuous cleaner cell whose feed 
comprises the rougher concentr/ita. ffhe analysis is then 
extandfid to tho data, of Woodbum et al for gangue component. 



^ OBJJliCTIVE I'UNGTIOl FOR ESTBliTIOI 


OF FLO'lATIOiT RATI 


DISTRIBDTID.il 


Identification of the flotation rate distrib'ation requires ' 


laiowledge of the 
parameters a and 


S'fcScCiy 8“Cat;0 flo' 
p 5 for lurarv feed 


• M;j(l) or Mq(L) , as v/ell as 
s.Lurry flo?/ rates. Discretization- 


of the inte.gral term in equation (7.2), followed by either solution- 
of the resulting algebraic e.quations or application of the 'least 


squares technique should, in . theory, provide an estimate of the 
distribution. It turns out, ho?/ever, that because of the ill- 
conditioned nature of the matrix in the discretized form, this 
estimate may prove to be unreliable, frequently physically 


meaningless, in presence of even a modest level of errors v/hich 
is inevitable in any experiraental data. A numerical method 
similar to the one used for inversion of the Laplace transform 
in Chapter 5 was developed as follows: 


let Rj,(Za,L) be the cumulative retained distribution function 

which is related to the freaueneq/ function by 

6R^XKa,L) 


lSj(Ka,I.) = 

Substitution in equation (7.2) results in 


(7.3) 


Zq(u) - j - 
0 

Integration by parts leads to 

IUq ( J-i ) = 


OR (Ka,l) _ 

[1- 




a+1 


(Ka+p) 


a+1 


] dKa (7.4) 


By (Ka ,1) 


C a+l ) 


. (7.5) 



where for corvehienoe Rj(Iia,I) has been nontallzed by Slviding 
throughont by slurry mass flow rates of particles of size 1 
in feed, iransforrn the variable ICa to x by 

Ka = 1 c\ 

+X-' S ( 7» D / 

.fherej as defined earlier, s is a scale parameter (s y O) 

¥/hichi can be arbitraril’* adjusted to generate estisate of 

Rj,(ta„,L) in any region of interest in range 0 <. Ka < Substi- 

tLition of equation (7.6) into (7.5) results in 

^ S" (x) 


M (L) = SisfJl p0:+1 / 

-1 


[ 


d x / o 7 

0+9 / s P w • f ; 


where 


F (x) = = Rj,(Ka,L) 


(7.3) 


discretization of integral terra in equation (7.7) by n- point 
bobat to quadrature 5 3 , 1 p6 


leads to the follovri.ng expres 


sion 


M (L) - < 


! 2n(n-l) p 


a +2 


n- 1 j' 

+ 1 i ' 

1:= 2 ' . L 


.1 


I 

■.11^4 ., ^ e 1“'^^ 

,;i+xps » M ^ 


■ n ■ . 

O * W, 


(7.9) 


"i' ~ -j 

vdiere x^ and are i-th abscissa and weight, respectively, in 
n-point bobatto quadrature formula. For determination of 
f"(x) (i.e. Rj,(Ka,I-)), the follovdng least squares objective 


function is found to be quite satisfactory. 
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N f 

Z ^ M (L) 
2=1 I ^’3 


, a . +1 

'2( ^ ' 


n~i 

E 

i=2 




fc+2 

2n(n-l)p.3 

J 


V?. 

1 


i 2 


1~z, ) 


i ;ri +i"ni j 


a^+2 ( 1+x^. )- I 


(7. 10) 


where subscript j denotes the values of recovery and residence 
time parameters at j-th slurry flow rate, j = 1,2,,.., IT. The 
problem is nov/ reduced to minimizing the objective function 


with appropriate choice of P (x) . The constraints on Rj,(Ka,L) , 
which are identical to those defined by equations (5* 22) , (5. 23) 
and ( 5 . 24 ) and are reproduced below for convenience j considerabl 
simplixji" the evaluation of P '(x). 


(i) 

Rj.,(Ka ,L) > 0, 

0 < Ka < “ (7. -11) 

(ii) 

-R-giCKa,!) is a 

monotonically non-increasing function 

of Ka . 

Therefore s 



1 2 R27(Ka-, .L) ^ 

„d Rp(Ra 2 . 2 Rji( KSp,, 5 E) •••2. ^ ( 7 * 12 ) 

where 

0 < Ka, < ICao 

rn 

(iii) 

Finally , 



RpCO,!) = 1 

' : 'i' (7.13) 

and 

Rp(°°,L) = 0 

(7.14) 
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CO&iPUT'AT lOKlIi ASPECTS *. 


Again a 8 point (n = 8) Lobatto quadrature was chosen 
to aiscretize tee integral in equation (7.7). As pointed out 
earlier, the choice of this foramla is dictated by the fact 
th£;t it is exact at the limits of the integral and the cons- 
traints in equations (7.13) and (7.14) are automatically sati^ 
fiea. Apatite and gangue data for steady state concentrate 
output as a j. unction of the feed slurrj' flow rates were taken 
from Table 6 of vtoodburn et al®^. The number of data points, 
vAich for different particle sizes, ranged from 6 to 9, v/as 
deemed co be insufficient. In order to smoothen out any 
scatter in the data and over determiie the unlcoown variables, 


a third order polynomial was fitted to the data set and 
additional points were iiterpolated from this curve. Tables 7.1 
and 7.2 give the Woodbum et al data in terms of recovery of 
apatite and gangue, respectively, as a function of throughput 
for different particle size groups after sTioothening and inter- 
polation. Total number of smnothened and interpolated data 
points ranged frora 11 to 17 for diffe::cnt particle size groups. 

Erom data given in Table 2 of Ifcodburn et al®^, plots of 
— 1 —0 3 

mean and variance” of particle residence tire distribution 
as a function of feed slurry flow rate for various particle size 
groups were obtained and are shown in Figs. 7.1 and 7.2, respec- 
tively. Using these two figures and eqt^tions (2*24) and (2.25) 
Values of parameters a and p in residenoe time distribution 
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Table ?. 


Particle 
range ^ 
.micron' 


37 to 63 


67 to 88 


i: Practional Recoveries of Apatite as a function 
Ql Particle Size Range and Slurry Plow Rate. 


size Slurrj/ flow 
rate 

GU f 


5 . 70 
6. 10 
6»50 
6.95 

7.40 

7.80 

8.20 

8.80 

9.40 
10.05 
10.70 
11.30 
11.90 
12.40 
.12.90 
■1 3. 90 
14.90 

•5.70 
' 6 . 10 
'6.50 
6.95 

7.40 
7.80 
8.20 
8.80 


9.40 
10.05 
10.70 
11.30 
1 1 . 90 
12.40 

12.90 

13.90 

14.90 


Experi- 
mental ■ 

Polynomial 

Pit 

ecp_yej?X- 

Predicted 

0.585 • 

0.570 - 

0.528 

— 

0.540 

_ 

0.486 .. 

0.513 

0.493 

— 

0.436 


0.475 . 

0.462 

0.458 

0.415. 

0 . 443 


0.426 

0.432 

— 

0.404 

0.394 

0.336 

0.397 

- 

. 0.370 


0.357 ■ 

0. 356 

0 . 365 

- 

0.345 

G.3S0 

0.334 

; 0.339 

- 

0.325 


0.284 

0.316 

0.321 

- 

0.296 


0.276 

0. 269 

0.2S9 

•0.613 

0.605 

0.597 

- 

0.593 


0.558 

0.579 

0.562 

- 

0.563 


0.553 

0.547 

0.537 

- 

0.532 


O’. 534 

0.517 

0.500 

- 

0.494 


0.453 

0.471 

0*464 ' 

■_ 

0.447 

. 

0.415 

0.423 

0.429 

- 

0.402 


0.410 

0.382 

0.402 


0.367 


0.336 • ' 

0.353 ' 

0.362 

■ ■- . 

0.329 

: 

0.314 

0.312 

. 0.332 


Gontd 
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Contd. . . 

Particle size 
Range , 
micron 


88 to 125 


125 to 177 


Slurry flow 
rate 

cu ft/min 


.-™_™,„£rajC t ioiial_^e_c__q y e ry 

iiixperi- Polynomial Predicted 
mental Pit 


5. 

70 

0 

.590 

6 . 

10 



6 . 

50. 

0 

.543 

6 . 

95 



7. 

40 

0 

.497 

7. 

80 



8. 

20 

0 

. 458 

8. 

80 


9. 

40 

r\ 

( 

.408 

10. 

05 



10. 

70 

0, 

■ 363 

11. 

30 



11. 

90 

0, 

.330 

12. 

40 



12. 

90 

0, 

.307 

13. 

90 



14. 

90 

0. 

.270 

5. 

70 

0. 

.490 

6. 

10 


— 

6 . 

50 

0. 

. 448 

6 . 

95 



7. 

40 

0. 

422 

7. 

80 


— 

8. 

20 

0. 

385 

O 

u* 

80 



9. 

40 

0. 

329 

10. 

05 



10. 

70 

0. 

271 

11. 

30 



11. 

90 

0. 

257 

12. 

40 


— 

12. 

90 

0. 

211 

13. 

90 



14. 

90 

0. 

107 


0.576 

0.546 

0.552 


0.530 

0.511 

0.507 


0.486 

0.477 

0.468 


0.451 

0.450 

0.429 


0.408 

0.414 

0.388 


0.370 

0.381 

0.355 


0.341 

0.356 

0.331 


0.320 

0.337 

0.301 

.... 

0.284 

0.301 

0.488 

0.473 

0.472 


0,455 

0.430 

0.436 


0.416 

0.406 

0.399 


0.382 

u. 369 

0 . 356 

— 

0.331 . 

0.332 

0.305 


0,281 

0.290 

0.260 

... 

0.242 

0.266 

0.229 


0.218 

0,230 

0.202 

— 

0.197 

0.210 


Contd 
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Gont cl. . . 

Particle size 
range, 
micron 

177 to 250 


> 250 


Slurry flow 
rat e , , 

cu ft/min 


Fractional 1 

le CO very 

Experi- 

mental 

- Polynomial 
Pit 

Predicted 

5.70 

0.422 

0.409 

. 0.375 

6. 10 


0.369 


6.50 

0 , 295 

0.337 

0.320 

6.95 

- 

0.309 


7.40 

0.339 

0.288 

0.291 

7.30 


0.274 

— 

8.20 

0.241 

0.262 

0.269 

S.80 

— 

0.245 

— 

9.40 

0.223 

0,226 

0.243 

10.05 

— 

0.196 


10.70 

0.155 

0. 153 

0.199 

5.70 

0.219 

0.218 

0. 177 

6. 10 

— 

0. 189 


6.50 

0. 163 

0. 164 

CO 

• 

o 

6.35 

— 

0,139 

— 

7.40 

0.112 

0. 117 

0.117 

7. 80 

— 

0. 100 

— 

8. 20 

0.094 

0.085 

0. 107 

8.80 

— 

0.067 

— 

9.40 

0.047 

0.052 

0.067 

10.05 

— 

0.038 

- 

11.70 

0.028 

0.027 

0.050 
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Table 7. 


Particle 
range , 

micron 


57 to 63 


63 to 38 


Fractional Recover-^os of Gangue as a Function 
of Parcicle Size Range and Slurry Plow Rate- 


size Slurry flow 

rciXo g 


Ixa ct iona 1 Re co v ery 


1 f' 

c/m in 

Ex perl- 

Pclynomial 

Predicted 



mental 

Pit 



5. 

6 . 

70 

10 

0.361 

0.345 

0. 322 

0 

, 30 s 

6 , 

50 

. 0.277 

0.501 ^ 

0 

,279 

6 . 

95 

- 

0.280 


'7. 

40 

0.258 

0.261 

0 

^253 

7. 

80 

- 

0.245 


8 . 

20 

0 . 233 

0.231 

0 

234 

8 , 

SO 

— 

0.213 


- 

9. 

40 

0.206 

0. 198 

0 

209 

10 . 

05 

- 

0.183 


. 

10 , 

70 

0 . 187 

0. 171 

0 

1S3 

11 . 

30 

— 

0. 162 



11 . 

90 

0. 142 

0. 154 

0 

172 

12 . 

40 


0. 148 



12 . 

90 

0.135 

0. 142 

0 

161 

13. 

90 

- 

0. 132 



14. 

90 

0.125 

0. 121 

0 

142 

5« 

70 

0.369 

0.363 

0 

354 

6 • 

10 


0.343 



6. 

50 

0.313 

0.325 

0 

314 

6* 

95 

- 

0.306 



n 

( • 

40 

0,292 

0.28S 

0 

277 

7. 

80 


0.273 



8. 

20 

0.259 

0.259 

0 

256 

3. 

SO 


' U. 24 O 


■ 

9. 

40 

0 . 22 S ' 

0.225 

0 

229 

10. 

05 


0.206 


•«*, . 

10. 

70 

0.191 

0. 192 

0 

205 

11. 

30 


0. 180 



11. 

90 

0. 164 

0.170 

0 

180 

12. 

40 


0.162 


— 

12. 

90 

0.159 

0. 156 

0 

170 

13. 

SO 


0. 145 


_ 

14. 

90 ■ 

0.136 

0. 136 

0. 

147 


CoBtci 
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Contd o . . 


Particle size Slurry flow 
range, • rate, 

micron cu ft/min 


5.70 
. 6. 10 
6.50 
6.95 

7.40 

7.30 
8.20 
8. 80 

88 to 125 9.40 

10.05 

10.70 

11.30 

11.90 

12.40 

12.90 

13.90 

14.90 

5.70 
6.10 
6 . 50 
6.95 

7.40 
7. 80 
8.20 
3.80 

125 to 177 9.40 

10.05 

10.70 

1 1 . 30 

1 1 . 90 

12.40 

12.90 

13.90 

14.90 





fractional Recovery 

Experi- 

Polynomial 

Predicted 

mental 

Pit 


0.315 

0.315 

0 . 285 

- 

0.294 


0.273 

0.275 

0. 258 

- 

0.256 


0.242 

0.238 

CM 

• 

o 

" 

0.224 


0.215 

0.211 

0.215 

- 

0. 194 


0. 171 

0. 180 

0. 192 

- 

0. 167 


0. 164 

0. 156 

0. 172 

- 

0.148 


0. 134 

0.141 

0.157 

~ 

0. 136 

— 

0.137 

0.132 

0. 146 

- 

0. 125 


'o.n? 

0. 118 

0. 129 

0.256 

0.252 

0.219 


0.232 


0.21C 

0.214 

0.196 

-- 

0.196 

— 

0.175 

0. 179 

0.169 

- 

0. 165 

— 

0.153 

0.153 

0.155 


0.137 

■ — 

0.127 

0. 123 

0.134 

— 

0.111 


0.104 

0.101 

0.117 


0.093 

— 

0,087 

: 0.037’ 

C . 104 


0.083 


0.075 

0.030 

0.094 

— 

0.075 

- ^ 

60*075:.': 

0.073 

0.0901 


Contd 
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Contd. . 


range , 
micron' 


177 to 250 


> 250 


Slurry flow 
rate, 
cu ft/min 

fractional Reco' 

very 

Experi- 

mental 

Polynomial 

fit 

Predicted 





■ * ■ 

5.70 

0/181 

0.177 

0. 146 

6.10 


0.157 


6 . 50 

0- 128 

0.141 

0. 130 

6 . 95 


0.126 


7 . 40 

0. 129 

0.115 

0. 116 

7. SO 


0.107 


8. 20 

0.096 

0. 100 

0, 106 

3.80 


0.902 


9.40 

0.073 

. 0.080 

0.094 

10.05 

— 

0.066 


10.70 

00 

o 

0 

o 

0.047 

0.073 

5.70 ■ 

0.064 

0.064 - 

ICS 

6 

d 

6. 10 

— 

0.056 l 


6.50 

0.049 

0.049 '■ 

0.040 

6.95 

~ 

0.041 


7.40 

0.032 

0.035 

0.035 

7 . 80 

— 

0.030 

— 

8. 20 

0.030 

0.026 

o 

o 

3.80 

— 

0.02t 

■ 

9.40 

0.013 

0.015 

0.025 

10.05 


0.010, 


10 .70 

0.006 

0.0C6 . 

0.021 











0.125 - 0.250 
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could be obtained at any desired slurry flow rate. The values 
of aj bnua obtaineds^ were found to be quite small. This shows 


that tus cell nearly behaves like a perfect mixer. Three 
optimization algorithms^ univariant search using ribonacci 

technique"^- Rosenbrock''59 and Rardom Search were 

tested. Althoughs Random search turned out to be the most 


sa ti siac b orj ' 5 it is possible that a still more efficient 
method might exist. In order to check against possible reali™ 
zatioii of local minima j widely different initial guesses were 
employed. The scaling factor was varied from 1 to 5 . 


FLOTATION RATE DISTRIBUTION: 

Figs. 7*3 and . 7. 4 shov^r plots of computed cumulative flota 
tion rate distributions, RyCKa,!) for apatite and gangue for 
six particle size groups. These distributions were employed to 
back calculate fractional recovery as a function of throughput* 
Figs. 1 , 3 , 7 . 6 , 7 . 7 , 7.8 and 7.9 show the plots of fractional 
recovery as o function of throughput for apatite and gangue for 
different particle size groups. It can be seen in these figures 
that the overall agreement between the calculated recovery 
Values and experimental values is satisfactory. 


The Validity of the throughput invariant distributions 
in Figs. 7.3 and 7.4 has been, tested in a somewhat wore , 
stringent manner by predicting the recovery when the rougher 
concentrate is refloated in a. steady .state cleaner cell. 
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following the argument used to derive equation (7.2), it can 
le shown that the mass flow rate in concentrate streara of cleaner 
cell is 

a^+l 

cx) R ^ 

Mo od) = I lie i(Ka,I) [l- Toa+f ] «<a (7-15) 

0 ’ (Ka+p2) 

where subscript 1 and 2 refer to rougher and cleaner cells, 


respectively. With changed notations the recovery of particles 
?jith attribute , from equation (7.2), can be vo’itten as 


M^^^(Ka,l) = lipdCa,!) [i 


P 


a^+1 


.) 


a.j+1 


1 


(7. 16) 


Substitution of equation (7. 16) follo’i'ved by substitution of 
equation (7.5) into (7.15) results in 


= f 


8Bj,(Ka,l) 


eKa 


[1 


a.+1 

3 

11 


a.+1 

(Ka+fi,) ■ 


J 


ag+l 


"'1 - 


Ug+I 


1 1 




(7.17) 


. (Ka+fjg) 

Integration by parts leads to 

(a.+l)P-'' (ap+l)p? 

Ej,(ra.i)[— ' 

(Katp^) (Ka+pg) 


a.+1 , . 


2 (^) 


I 

0 




ll 


. a.42 ^ ' ' 

"" -(Ka+Pg) 




(7.18) 
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CtiaBge of variabl© froa Ka "bo tising equation 1.S v/ith s = 1 
re stilts ill following eq-ua,tion 5 ' ^ 

, CC^-rl tto+l 

. s J (a<5+l)pQ'^ 

lU JL) = 2 f F (x) [ — I 

,) ^ (S+pO 2 


-1 


1 


^ 1 O'r) “^1 j 


'-] p2 ■ \ ( + ( a2’^1 ) 

- «-* a ^^2 . , (x^i^2 * J . 

(x+p^) (x+pg) 


dx 

( 1-fX)' 


(7.19) 


I 1 

where S' (x) is same as defined in equation (7.8) and 


X = 


IxS 

'^^^x 


(7.20) 


Discretization of the integral term in equation (7. 19) by 
8-point Lobatto quadrature results in 


K 


a .+1 

(a^+l) '' '-^2' 


ao+1 

^ ^ (o:2+l)p2" 

0 (1)=2 2 F(xp[--f TTs 

' (Xj^+pj) ® 


Pg’ { (“i+n('q+P2)_+ (gyH)(X j^.p^)| - 


0^+2 _ 

(Xi+p^) (x^+pg) 




w. 


'( 1+x^)' 


(7.21) 


Substitution of the numeric values of F (x), computed earliei 
in estiiiation of the flotation rate distribution, in the above 
equation gives the predicted fiactioml recovery 
in the cleaner cell. Experimental and predicted fractional 
recovery values as a functiou-Of ^utry flow rate for apatite 



for different particle size groups are shown in Fig. 7.10. 
Similar plots for gangue are shown in IP ig. 7.11. Por reasons 
unknown 5 predicted recovery values in cleaner for particle 
size > 0.25 mm in the case of apatite and for particle size 
0. 177-0.25 mm in the case of gangue show considerable deviation 
from’ experimental values. Comparison of Tfoodburn et al model, 
based on throughput dependent Ka., mth the present throughput 
invariant first order distributed rate constant approach sho?/s 
that, on the wiole, both lead to comparable results as far as 
the hack calculated rougher recovery is concerned, but in the 
cleaner float the latter model turns out to be somew/hat more 
accurate. We may therefore conclude that based on the data 
that has been scrutinized, the simpler kinetic representation 
is puite adeq’uate atleast for modelling p’urposes, as compared 
to the model of Woodbum et a 1 where in order to obtain an 
acceptable fit, it was feund necessary to ii’itroduce slurry flo^. 

rate in the expression for the, flow rate constant. 


EPPEC-1 OP PAR’TIGIiE SIZE: 

The computed flotation rate distributions provide a 

good opportunity to test the validity of equation (6.14) 
which was originally proposed for semi-hatch flotation operation 
When RjCKa.l) is equal to 0 . 5 , substitution of equation (6.14) 













Particle Size 



v.Ak!.* 


fV/VK* 




WM?B 




■'4;-,. \^Vi 6 ,i',''' 'V'.v'.. V •;'(;/'••, ''i'Vv 


IWM 

TET?./ 











140 


where iva.(0.5 5 lj) is the median apparent flotation rate constant 
for particles of size L. As discussed in the last chapter 
(equation > 0 . 22 ), n(l^) is proportional to Ka(0.5jl) and the 

constant of proportionality K°(0.5), is^ by definition, 
independent of particle size. Therefore, equation (6.14) 
becomes 

_ ■ K° (7.?3) 

Ka(0.5,l) K°(0.5) 

This implies that plots of Rp(Ka,L) as a function of dimen- 
sionless flotation rate constant Ka(L)/I[aX0. 5,1) should be 
independent of particle size, provided the hydrophobic it, y of 
the species (i.e. , S°) remains unchanged with particle size. 

A mathematically rigorous derivation of this result is given 
in previous chapter. Pigs. 7*12 and 7. 13 show the collapsed 
distributions for apatite and- gangue. The values of the median 
flotation rate constants v;ere read off from smooth curves in 
Figs. 7.3 and 7 . 4 . Apatite results are renia.rkably good, seme 
deviation in gangue data points pertaining to particle size 
> 0.25 miii could le due uo either experimental and computational 
errors or change in hydrophobe city of large particles. The 
collapsibility of the flotation rate distributions has a 
significant implication, namely, other things equal, the 
apparent rate constants in a distribution are modifieu to the 
same extent with change in, particle size. 
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to the 


Sinc.^ tti'^ median rate constant KaC.O.SjL) is proportional 
j-UT-etioii j(li) 5 it is now possitle to deteiTnine the 


relationship between 


the lat uer function and the mean particle 


size wii-h an uncertainty of a constant of proportionality, 
^io* 7*14 snows na(0.5si') as a function of the niem particle 
size for both apatite ana gcngue. As expected, in both the 


cases, tnere exist characteristic peaks at inteimediate size 
range where the flotation rate is maximum. Ihese curves, 
riioreov'cr , show bhat optimal separation between apatite and gangue 
would occur when the particle size is about 0.075 ram and negli- 
gible flotation is to be expected when it exceeds about 0*45 mm. 
From equation (6.14), it follows that the ratio of the median 
flotation rate constants of apatite and gangue is invariant of 
particle size. Consequently, if the ordinate of the gangue curve 
in Fig. 7. 14 is approxiriately rescaled, the resulting curve 
should match the apatite curve. -The broken cur/e in this figure, 
scaled by a factor of ,3, is in fair agreement vdth that of 
apatite except in fine particle size range. Due to possible 
flocculation and agglomeration of fines the actual situation is 
perhaps more complicated than what is allowed for in an admittedly 
highly simple expression in equation (6.14). 

OOKPUIATIOI'I OF FLOIAIIOI MIS DISIRIBUIIOHS EiFlTIjS: GASS: 

The method used here, for 'determining the flotation rate 
distribution from steady state ■ flotation kinetics _per tain ing to 
different slurry flow rate.s,'^® dictated' by availability of 






145 


extensive data in the literature. l?ith slight ^'modification, 
this technique may oe readily applied for response data to a 
tracer impulse in a continuous cell which would be more conve- 
nient provided a suitable tracer is available. As far as the 
finr.l o oU 1 1 o Cij-C Guiicemed Woodbum et al have alread"''^ decions- 
tra. teo. the: e c^cii le nee of stead 3 ^ state and dvnaiHic tests. 
Briefly / 5 the computation of algorithm in the impulse case is 


as folloYmrs: In response to a tracer impulse, the mass flow rate 
in tailings stream time t later is 

lVi2,(Ka,L,t) = ECt) Mj,(Ka,L) exp [-Ka t] (7.24) 

or 

oo, 

M,-(L,t) = E(t) / My(Ka,l) exp [-Ka t] dICa (7.25) 

0 

Substitution of equation (7.3) leads to 

== 0R-(lCa,l) 

Mj(L,t) = E(t) / - — — 1 ^-— . exp [-Ka t] dKa (7.26) 

.Integrating by carets ■ 

CO 

M^(L,t) = iKt) [l-t/ dKa] (7.27) 

o 

Change in .-variable from Ka to x leads to 

M..n(L,t) = S(t) [l- / r(x) (x^.l)®'^"’’ dx] (7.28) 

h r _ ■■■ ■: ' . ■ g - 

where the transformation variable x is same as defined in 
:btmtion r(5.9,) , viz. /ggg'grgh!:;- ;yghgr/g^g 
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X > 2 

exp [_^ 

f- ] -1 

(7.29) 

and 




b(x) = 

R( s In 


(7.30) 

Niscre t isat ion 

of the 

integral term in equati 

on (7.23) by 

n-point Lobatto 

quadra 

turc fomiula gives 




It is now possible to set up the least squares objective function 
using data at different times vAiich Eiay then be optimized vfith 
suitable choice of unknovra I'(x), Although, the impulse technique 
is experimentally more efficient, steady state approach may have 
to be used if suitable radioactive tracer is not available. 


flSCUSSIOil AITD CONCLUSIONS: 


Given a flotation cell, suitable experiments can be designei 
to determine th; reoidence time distribution of particles and the 
apparent flotation rate dirjfcribution by either d:/n8mic testing 
with a tracer napulse or by monitoring the steady state output 
£t different slurry flow rates. The apparent flotation rate 
distribution turns out to be independent of the slurry flow i*ate, 
results, which to a first approxination and at least for modelling 
purposes', are intuitively valid provided, of course, the hydrody- 
namic regime is not drastically altered with the slurry flow rate. 
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The effect; of parficle size on the apparent flotation rate 
constant is similar to that in semi-batch case and hence 
the analysis of notation kinetics in semi-batch case can bo 
extended to a' continuous configuration. Altho-ngh suitable data 
is not available, the effect of the aeration rate in contin-'-’”’- 
cell may also turn out to be the same as in sem5--batch cell, 
in view of the results shown in Chapter 6. More significantly, 
the kinetic model for continuous cell has been demonstrated to 
possess sufficiently accurate predictive poi;rer and is deemed 
suitable for modelling, control and optimization purposes. 'Tho 
latter aspect is discussed in the next chapter. 




PART ..II 


CHAPTER 3 

OPTIMAL AHL SUBDPTEIAL SYIlTHBSIS A.^I) DESIST OP 
PLOTATIOn CIRCUITS 

'THe iraportance of synthesis and design of chemical 
reactor units for optimization of the efficiencj?- of the 
process has been touched upon and some of the generalized 
techniques, developed so far, have been discussed in the 
literature reviev/ in Chapter 2. Whenever, two or more 
reactor units are to be interconnected to perfoin a task, 
the problem of synthesis of reactor units comes into play. 

Por simple systems trial and error approach to synthesis 
can be applied, but it becomes progressively cumbersome 
and eventually impossible as the number of reactor units in 
the assembly increases. 

In this ch.apter wc have carried out tentative investi- 
gation of the foasibility cf sjcnt bo sizing and designing 
optimal and subop -cimal multicell flotation circuits, employing 
a number of appropriate objective functions of recovery and 
grade, and profit, by a direct optimization approach due to 
Umeda , ■ Hirai and Ichikawa First, a flotation circuit 
with t?/o cells and feed comprising of two species of known 
overall rate constants' is synthesized alid designed, and the 
results are then compared with the real optinuss circuit 
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obtained by using the trial and error technique. Second, 
the a.pproa oil is extended go a flotation circuit v/ith four 
cells and for batteries of cells for a feed comprising of 
three species. It v/ill be seen that the direct optijrization 
techniciue seems to be quite promising for 3:ynthssis and 
design of the industrial flotation circuits with large number 
of interconnected cells and feed comprising of a number of 
species which ccn lead to better performance and efficiency of 
the flotation process, 

SlAlMIiMl 01' SHE PROBIIM: 

I'or illustration purposes, consider two interconnected 
cells. Pig, 8. 1 shov/s six possible intuitively meaningful 
configurations in which, for the time being, split streams 
are excluded. In the trial and error approach the problem is 
broken into t?/o parts. First, all possible configurations 
are enumerated, then the best design parameters are determined 
for each structure, and the optimum configuration is selected 
by comparing tho values of the objective functions. In general, 
a flotation circuit cemprisos of many cells, performing 
rougher, cleaner and scavanger functions and it is quite 
tedious to enumerate all meaningful structures, koreover, 
there is no a priori justification for ignoring split 
streams and the enumeration problem, thus becomes even more 
complex. In the integral approach* proposed by bmeda, Hirax 

and Ichikawa'^^^ , enumerations- ar'e totally circumvented and 
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the optiiiial s uTucture and dasign parameters are extracted 
simultaneouniy frora a generalized circuit by direct search. 

This technioue has been applied to solve linkage (including 
split streams) problem and to determine the optimal residence 
time distriDution of particles in each cell. The residence 
time distributions as a design parameter, is an automatic 

b 1 

engineering choice since, as pointed out by Woodburn et al , 
only the bolding time is at present capable of being included 
into a quantitative description of the flotation kinetics. 

These authors had also earlier determined the optimal residence 
times in perfectly mixed cells but for an assumed circuit confi- 
guration. 

PLOTATIOli KDiETIG MOML: ; 

Erom the literature reviewed in Chapter 2, it is evident 
that for all practical purposes the kinetics of flotation can 
bo approximately described by a first order expression with 
distributed flotation rate constant. For industrial modelling 
purposes it is also knov/n from data of ¥foodbum et al that the 
cell behaves like, more or less, a perfect mixer ana the mean 
residence time is a function of particle size and density. 

Thus, if Mj.(ICa) is mass flow rate of particles of overall 
flotation rate constant iva to the cell, as defined earlier, 
and t is the mean residence time in the pwlp, the mass flow 
rates in concentrate and tailings , stream at steady state are 
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IVinU' 

U 


(Ka) 


Mj,(Ka) 


Ka t -j 
l+Ea t 


(S. 1) 


M (Ka) = IU(Ka)[ 


1+ICa t 


( 8 . 2 ) 


Using these relationships it is now possible to construct a 
stead^f state mathematical model for n^ particulate species 

O 

(discretised apparent flotation rate distribution) of same size 
and density in a flotation circuit ?/ith n cells, for simpli- 

C/ 

city, the treatment is restricted to single size particles, or 
more correctly to particles of same value of t. Relaxation 
of this constraint is discussed later. It will he noted that 
water balance of circuit, while important, is not really germane 
to the problem, as each cell is normally eq.uipped with arrange- 
ment for addition of water (and surfactant) to bring the pulp 
to required consistency and level. A constraint, given in 
the seq.uel, will assure that the pulp density never exceeds 
a limiting value, say 20 percent. 

Let the specific flotation rate constant of the j-th 
species be Ka. and f. be its fractional valuable mineral 
content, and t^ be the mean residence time of solid particles 
in .i~th cell. In the generalized configuration, the feed is 
split into n streams and fed to each cell. Similarly* eeeh 
concentrate and tailings stream is split into (n^+l) streams 
and connected as feed to every othef cell as well as short 
circuited back to the cell fr<HB/W|i^fe, it originated 


split into n_ streams and 

■ V/ ■ . ii' ' ' 





(n^+l)-th split leaves the circuit as concentrate or tailings 
output, as the case may be. I'he generalised configuration for 
2 cells 13 sho\¥n m Ilg. 8.2. Consider now the i~th cell. The 
total feed flow rate of j-th species to this cell is 


'Di 


= E 


‘>fl 


■ 

+ 2. M 


n, 


ki ^-Cok hi 


(8.3) 


where is flow rate of j-th component in new feed, 6^^ is 
fracbion of this feed going to i-th cell, is tailings flov^ 

rate of ,j~th species from k-th cell, 6^. is split fraction of 
this flow going to i-th cell, is concentrate flow rate of 

0 -th species originating from k-th cell and is split fractio 
of concentrate flow from k to i-th cell, let in equations (8.1 
and (8.2) 


a. 


1 


l+Ka^ t. 

J J- 


(8.4) 


Therefore, equations for flow rates of j-th species from i-th 
cell in concentrate and tailings streams can be rewritten as 


and 


"‘Cji 


Mm . . 

Tor 


ji 


(8.5) 


zr ' P .( fY . . ■ 

01 ^01 


Combining equations (8.5) and (8.6) 


( 8 . 6 ) 
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Substitution in equation (8.3) gi 


give; 


^’fi 

I’-c 
‘ £ 
k=1 

^^310 

Ir 

1 

""dk ■ 


^ki| 

(8.8) 

e qua t ion 

(8. 

6) result 

s in 



i 




nc 


1-a.T 



8fi 


“31 

£ 

k=1 

^S3k \ 

r 

. «3lc 




i 

= 1 

9 

j <1 5 • * • 






1 

= 1 

$ 2 . ^ • • • 

’ ^c 


(8.9) 


These (n +n^) linear algebraic equations can be solved for Mn- - 

-1.31 

in beriiis of linkage and design parameters. The formal solution 
can be va? it ten as 


M^. . 


Si 


i's, a's) 


( 8 . 10) 


Similarly from equation (8.7) 

M... = az 5 .. 

C 31 ' or L J 


(8.11) 


let 5. and 6 . be the fractions of tailings and concentrate 
froc; i-th ce.ll, respectively, leaving the circuit, then, by 
definition (iraetion of concentrate/feed flowrate.)} 'fcba 
percentage overall recovery is, 


“c ^s 
£ £ 
1=1 3=1 

’ : n„ ; 


'Cji ho 


■'s 

£ liiT 






( 8 . 12 ) 
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of overloading of bubbles and loss of flotation selectivity, 
the pulp deiisicy 7)^ (Percent mass of solid in slurry) should, 
no V exceed a prescrioed value d^ in any cell. By definition of 
the mean residence time the hold up of solids in i-th cell is 


^i *^I';ii 

D=1 ^ 


(8.19) 


Therefore 5 the lower bound on volume of the i-th cell is 


V. > H. r '4- + 

1 “■ 1 L. j3 


100-d^ 


( 8 . 20 ) 


where P is solid density. Moreover, the overall size of the 

s 

plant is also restricted to some reasonable volume Vj, of all 


cells combined, hence 


c 

£ V, 


( 8 . 21 ) 


Following the suggestion of Woodbum et al, for convenience in 
compiitation 5 personal judgement is used to impose a q.uite 
liberal upper constraint t on t, i.e. , 


0 < t^ < t 


i — 1 } 2,.»«5 n _ 


( 8 , 22 ) 


To be acceptable, it is necessary that the grade of the concen- 
trate does not drop below a minimum level which gives rise to 


G > g 


(8.23) 


An obvious objective is to maximize recovery given in equation 
(8.12), subject to above constraints.. .In.view of the fact that 
as Ka^ > Kag ..., > Sa^ ^ 

V ..vf" dd': r/ Vi" n, .o'': ^ ; p-'' , 
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vre may infer fronj physical nafure of fhe prohle^n fhaf the 
optlm® lies on boundary of constant S = g. ihls. however, may 
not be nocessar 5.1y true for more complicated objective functions 

as illustrated subsequently, 

COMrUlAlIOPii ASPECTS: 

The integral approach used here requires search through 
a lather large parameter space. This gives rise to well Icnown 
difficulties in mattei. of slov/ rate of convergence and false 
optima. An allied problem is of deciding ?/hether a split 
stream iSj in factj missing or not when linkage parameters 6's 
and P's are near one or zero. On the otherhand, realization 
of a large number of suboptima, is not necessarily a disadvantage, 
since some configurations may turn out to be simpler to design 
with lower capital investment and more economical to operate 
than the global optimum structure with only a marginal loss in 
performance. 

The set of algebraic equations in (8,9) were solved 
by Gauss Jordan technique A systematic evaluation of the 
best optimization algorithm was beyond the scope of tills 
investigation. However, modified complex''^-' aiid Random 
search"^^*^ methods were employed (See Appendix B). Success with 
the first algorithm, used by Umeda et al^^^ also, was found to 
be highly sensitive to initial guess,, the latter method 
although slow, was quite reliable,-- every case a number of 
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dixfsrGii’t ini’tidl guGSSGs wgx’g upisd in opden "to assixcSj as fax 
as posoiblc. 5 a ijiutj opiiin'unio ins opciniizat/ion wqs caxxiGd oui 
in TWO ox liiore stages. In course of computation it was observed 
that eventually the improvement in the objective function 
became rather slow. The search vas stopped at this stage and 
linkage ratios which were less than or greater than pre specified 
values were taken as zero and one, respectively. Starting 
vdth the reduced generalised configuration, the search was 
renewed in the smaller parameter space. This stage wise search 
procedure was continued as long as the successive reduced confi- 
guration resulted in an improvement in the objective function. 

As pointed out by Hendry, Rudd and Seader , this integral 
approach Yiith direct stage wise search may not necessarily lead 
to the global optimum although for two cell case treated in the 
sequel the convergence has been verified by independent means. 

TVa SPECIES, TWO CELL PROBLEM; 

I'or domonst ration, consider first the case n^ = n^ = 2. 

This v/ill peinnit evaluation of the efficacy of the integral 
method by comparison with 'the enumerated structures, in Pig. 8.1, 
optimised individually with jpespect to the mean residence time 
in cells I and II. Table 8,1 lists the numeric values of feed 
characteristics and the constraints on the process design 
Variables' employed for this exercise, S&pansion of equation (S,9y 
for this particular case, and rearrangement of terns leads to 
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Table 8.1; Numeric Values Used iu Computation of Design 
X a j- amet ex s in Two Cell Probleii!, 


feed composition Constraints on Values of the 
percent variables parameters 


Valuable = 22 

Vj, < 20 cuft 

Ka T 
val 

- 1 

= 1 min 

Uaiigue =78 

dp < 20 percent 

^^gang 

= 0.1 min" 


0 < t ^ ,^2 < 20 

II 

25 Ibs/rnin 



% = 

10.0 units 



Pg = 

1.00 units 
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following two. equations. 




-t M 


&21 


Dl 

1-a. 

2l'- «32 


a^1 


and 


q = 1,2s. .. 5 Hg (8.24) 

-%j1 ^*^12 ^2 ^32 

’t ■*•**3]^ 2 [■'^22 *^32 ^ ^22 ^*%3 ^f2 ^12 ~ ^ 

3 = 1 ,2,. . . , r.g (8.25) 

Eq.uations (8.24) and' (8.25) can be solved simultaneously fo^ 

and in tenns of , 5's, p’s and a's to give 

vfhich, then can be substituted in eq_uation (8.11) to give 

M„... Once M^.. are knoTO, explicit expression for recovery 
031 031 

R and grade G in terms of linkage parameters 6 ’ s and p's and 
c 

design parameters a's can be obtained from eq^uations (8.12) and 
(8. 15 ) 5 a s follows 


R 


(:..^11 + ■^''^u_2l!3p.Jli!kllJlfc22 


100 (8.26) 




and 


G = 




X 1C\ 


Pgi 1 “csf ^10 * '“oi 2 “022’ P20 


(8.27) 
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where 


(1-a, J 


'Cjl 


6 

“d1^-3 


3 l®fl“j1 D-5 


22®.j2 " ^22^ ''"“o2^'J 


+ ^2 1^21 P21 -;ji; 


a. .W j = 1,2 (8.28) 


f- 

‘■‘Cj2 = ■■ a.2T%3 1 “f2“32 I1-Sll“3l - 

J • 


I'-a • 


+ [6^2 + P 12 (“— -^"yj a^-2| , .1 = 1,2 (8.29) 


and 


^3 ° 


[jS ^ ^ ocjj ^ ^ ( 1— ^ ) — lj L52^0Cj-| + P 2 -] f '•‘- j ■] 


'ii“ji ^ hr '“hi 

1-a. 


1-a^ , 
:i2 


i-a. . 

L^12h2 h2^’ a.. '^h2^ t^22 h2 h2^'’"h2^“^^ 


3 = 1,2 (8.50) 

Two objective fixnctions used rteTe, masiomize recovery subject to 
constraint that grade does not fall below (A) 35 and (B) 75 percent 
respectively. In addition, a profit objective function (c), given 
below, was also formulated. 

(8.31) 


(C) 


-+ (L- - G-^,) p,^ 


v/here p^, is price of feed ore of grade G^i and Pj^ is the inci*ease 
in price of the concentrate for every one percent iopiovenient 
in the grade. 

Results are presented in fables 8.2(A), 8.2(B) and 
8.2 (C) , where initial guess, optipal eohfiguration at first 



'Table 8.2(A) 


Values of Designed Parariicters for Stage 
Wise Optimization for 2 Cell 2 Species 
Problem in Case of Objective A. 


■^araeters 

i-i'iitial 

Stage 

3 t :? ee 

Opt 

ima 

'v^alus 

T 

X 

TT 

j... -i- 

IX. 

VlO 

^10 

0.525 

Gm 764 ■ 

1 . oc 

1 . 

00 

^11 

0. 150 

0.067 

0.00 

0 . 

00 

§12 

0.325 

0 . 16 S 

0.00 

0 . 

00 

^'20 

0.717 

0 . 00 

O.GC 

0 . 

00 ■ 

C 

""21 

0.101 

n / o 4 

W * ‘‘•f ,j ~ 

1.00 

1 . 

00 

^22 

0.182 

0.506 

0.00 

0 . 

00 

ho 

0.345 

0.184 

CnOO 

0 . 

00 

hi 

0. 135 

0 • 1 2 2 ■ 

0.00 

0 . 

JO 

h 2 

0.469 

0.694 

1.00 

1 . 

00 

ho 

0.661 

0.900 

1.00 

1 . 

00 

^21 

0.152 

0,00 

0.00 

0 . 

00 

h 2 

0. 187 

0. 100 

0.00 

0 . 

GO 


0 . 254 

0.00 

O.DC 

0 . 

00 

^ ° f 2 

0 . 74 -S 

i.oo 

1 .00 

1 . 

00 


6.04 

9. 100 

1 A 545 

13 . 

700 

'h 

5.238 

2.015 : 

3 . 9 - r 5 

3 . 

7 30 

h 

_ 

8.06 

:■ 11.90 

12 . 

51 

V r: 

- 

'.6 a; 6:. 59 ■ 

7.43 

■■n. 

' f * 

48 

lOVQXy 


■ 661449 ' 

6 61 ,99 

: 6 2 • 

01 

.de 


35.00 

35.00 

35 . 

00 



Table 


Table 8,2( 

B ) : Value s 

of Designed 

Parameters 

for Sta 


Optimiz 

ation for 2 

Cell 2 Species Pr 


ii'i Cass 

of Objective B. 


Paraiiieter 

Initial 

value 

‘^tage i 

Stage II 

Optimal 

val"-'3 

^10 

0.352 

0.00 

0.00 

c.oo 

^11 

0 . 320 

0.522 ' 

0.00 

0.00 

f ''2 

0.328 

0.478 

1.00 

1 . 00 

&20 

0 . 840 

0.973 

1.00 

1.00 

^21 

0.038 

0.012 

0.00 

0.00 

^22 

0.072 

0.015 

0.00 

0.00 

*^10 

0.643 

0.783 

1.00 

1.00 

i°1 1 

0.144 

0. 141 

0.00 . 

0.00 

Pl2 

0.213 

0.076 

0.00 

c.oo 

^20 

0.025 

0.002 

0.00 

0.00 

^21 

0 . 485 

0.619 

1.00 

1.00 

^22 

0.490 

0.379 

0.00 

0,00 

X 1 

0.078 

0.00 

0.00 

0.00 

^f2 

0.912 

1.00 

1.00 

1 . 00 

^1 

1.095 

0.996 

0.658 

0.65-3 

tg 

^•527 

3.6:50. 

14.032 

14.032 

^1 

"’2 


1.612 

1.420 

1 . 420 

4 . 

12.170 

18.500 

18. 5 

Recovery 


22.350 

24.867 

24. 6'.67 

Grade 

p.-. " . 

35.00 

35.00 

35.00 




fable 8.2(C); Values of Designed Jrarameters for Stage- 
tese Optimization for 2 Cell 2 Species 
Problem in case of Objective G. 


Parameter s 

Initial 

Value 

stage I 

Stage II 

Optimal 

value 

“io 

0.518 

0.908 

1.00 

1.00 

®ii 

0.377 

0.088 

0.00 

0.00 

®12 

0. 105 

0.004 

O.OC 

0.00 

o 

CM 

0.235 

0.006 

0.00 

0.00 

^21 

0.237 

0.443 

1.00 

1.00 

^22 

0.478 

0.551 

0.00 

0.00 

^'10 

0.651 

0.001 

0,00 

0.00 


0.236 

0.325 

0.00 

0.00 

^12 

0.113 

0.674 

1.00 

1.00 

^20 

0.670 

0.880 

1.00 

1.00 

P21 

0.049 

0.002 

0.00 

0.00 

P 22 

0.281 

0.113 

0.00 

O.Ov) 

^ f 1 

0.2M 

0.271 

0,00 

0.00 

^f2 

0.75! 

0.729 

.1.0C' 

1.00 

‘^1 

8.943 

13.460 

14,815 

14.806 

C\J 

4 • 353 

! 0.715 

0.822 

0.822 


— 

16.65 ’ 

17.31 

17.31 

^2 

“ 

y ;3.;33; 

■.; 2.69 ' 

2,69 

Recovery 


34.69 

33.630 

35.645 

Grade 


’ ■;:5?.-60:'. !i 

60.70 

60.70 

Price 


595.70 

409.76 

409.78 



stage and "the final optimal structvnes sjTithesised from the 
reduced configuration for three objective functions {k) , ( 3 ) 
and (0)5 respectively, are shown. The reduction was carried 
out by equating large and small linkage para-ueters (> 0.3 and 
< 0.2) with one and zero, respectively. As expected; in first 
two cases optimum occurs at grade constraint boimd::ry, but, 
interestingly, depending on the grade stipulation, the struc- 
tures of the optimal circuits are quite different, for objectiv^ 
(a) and (B) the resulting configurations are equivalent to 
configuration ITo. (liE) and (V) in Big. 8. l(Cell number 1 in the 
computed configuration refers to cell number 2 in Big. 8. 1 and 
vice-versa). The last column under each of the three objective 
presents results of trial and error approach obtained by 
optimising individually the six configurations in eig. B. 1 
with respect to residence times only. The agreement provides 
a partial check that it is possible, at least in this case, to 
derive an optimal circuit starting from the generalised 
configuration in Big. 8 . 2 . 

The twe optimal circuits are intuitively reasonable. 
Thus, when recovery is to be increased at some accrifice of 
the grade, as in objective (A), it is obvious that one of 
the cells must act as a scavanger (Cell Us 5 ‘ig* 8, l(vi)). On 
the other hand in objective (B) high grade restriction majces 
it necessary that the rougher concentrate be refloated 
cleaner cell (Coll I, tn-lig. 8.1 (V))» 
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i'HtlSB SPECIES FOUR CELL PROEM; 


From industrial flotation point of viev;, this is a more 
meaningful problem. Ihe additional particulate species. 


considered in this caee, is middlings in v/hica valuable mineral 
is locked in gangue matrix. Table 8.3 gives the rrorueric values 
used for maximizing the recovery subject to a minimum of 55 per- 
cent grade. As shoura in Fig. 8. 3? the generalized structure has 
a rather large number of parameters and attainneiit of convergence 
in a reasonable time is doubtful. Therefore, optirnlzaticn v/as 
carried out in four stages. At the end of a stage, all split 
stream ratios <_ 0.05 v/ere set zero and those >_ 0.95 were 
eq_uated with one. Each stage resulted in improvement in 
recovery. The optimal structure is shown in Fig. 8.4, the linko/'i 
parameters are given in the paranthesis. Table 8.4 gives stage- 
wise progress of the linkages and design parameters and the 
objective; function. 


It will be noted in Fig. 8.4 that in the cptimal confi- 
guration there is one split stresin ol value 0.17 only. In order 
to simplify the structure, it was decided to remove this stream 
and seek ne\' parameters for best recovery. Results of this 
stage of computation, in column 7 in Table 8.4, shew that the 
recovery has now fallen somewhat. Fin^ther reduction was carried 
out in stage 6 by eliminating cell 4 altogether Prom Aie 
circuit 5 since its computed volume was comparatively sjni'ill. <i/> 
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Table 8.32 iJuiseric Values Used in Computation of 
Ifesign Parameters in 4 Cell 3 S-pecies 
Problem. 


Peed Composition, 
percent 


Con at ra inti 


Values of the 
paraiixrL'ers 


Valuable 

= 15 

g = 35 


. — 1 

= 1 miri 

Middling 

= 20 

0 < t . < 20 

0 

"■^®mid 

= 0.1 rnin"'^ 



D ' ~ 1 5 k . 5 4 



Gangue 

= 65 

0 

CM 

tl 

^'^gang 

= 0.01 rair”'® 


Vr,, = SO cuft i'ij, = 25 Ibs/min. 














cootimiraticn for three species 
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■Tat'le 3,4: 'Values of Designed Parameters for Stage. 

Wise Op biiiiization for 4 Cell 3 Species 
Problem in Case of Objecti've A,. 


Parame- 

ters 

Initial 
va lue 

Stage 

I 

Stage 

II 

stage 

T j T 

Stage 

IV 

Stage 

stare 

Vi 

®10 

0.484 

0.031 

0.002 

0.000 

0.000 

0 . uOu 

0.000 

b1 

0.046 

0.068 

0.031 

0,000 

0.000 

0.000 

0.000 

b2 

0.180 

0.389 

0.352 

0.000 

0.000 

0.000 

0.000 

b3 

0.166 

0.403 

0.615 

1.000 

1.000 

1 . 000 

1.000 


0.124 

0.109 

0.000 

0.000 

0.000 

0.000 

_ 

®20 

0,541 

0.049 

D .016 

0.002 

0.000 

0.000 

0.298 


0.250 

0.251 

0.014 

0.000 

0.000 

c.ooo 

0.000 

622 

0.1 30 

0.367 

0.238 

0.743 

0.654 

0.680 

0.702 

^23 

0.036 

0. 165 

0.423 

0.000 

0.000 

0.000 

0.000 

^24 

0.043 • 

0. 168 

0.309 

0.255 

0.346 

0.320 

_ 

^30 

0.285 

0.015 

0.020 

0.000 

0.000 

0.000 

0.000 

^31 

0. 183 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

^32 

0.207 

0.800 

0.945 

1.000 

1.000 

1.000 

1.000 


0.099 

0.143 

0.035 

O.OOG 

0.000 

0.000 

0 , 000 

S34 

0.226 

0.03& 

0,000 

0-000 

0.000 

0.000 


^40 

0.507 

0.403 

0.543 

0.962 

1.000 

1 . 000 

** 

^41 

0.220 

0.045 

0.000 

O.OCO 

0.000 

0.000 


^2 

0,053 

0.072 

0.000 

o.coo 

0.000 

0.000 


^43 

0.219 

0.407 

0.457 

0.C18 

0.000 

0.000 

■ 

^44 

0.021 

0.073 

0.000 

0.000 

0.000 

0.000 

-- 


CO nil’ ci »> * • 
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contd. • . 


Parame- 

ters 

Initial 

value 

S ta 

I 

Stage 

II 

stage 

III 

IV“ 

stage 

V 

stage 

VI 

Pio 

0.371 

0 . 756 

0.653 

1.000 

1 . 000 

1.000 

1.000 

^11 

0. 156 

0.003 

0.000 

0.000 

0.000 

0,000 

0.000 

Pi 2 

0.244 

0. G 18 

0.000 

0.000 

0.00 c 

0.000 

0.000 


0.061 

0.239 

0. 342 

o.oco 

O.OOQ 

0,000 

0.000 

^14 

0.163 

0.004 

0.000 

0.000 

0.000 

0 . 00c 

- 

ho 

0.414 

0.990 

1.000 

1.000 

1.000 

1 . occ 

1.000 

^21 

0. 210 

0.007 

0.000 

0.000 

0.000 

0.000 

0.000 

P22 

0.245 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

P23 

0.084 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

^24 

0.047 

0.003 

0.000 

0.000 

0.000 

0.000 

- 

ho 

0.551 

0.582 

1.000 

1.000 

1 . 000 

1.000 

1.000 

P31 

0.021 

0.077 

0.000 

0.000 

0.000 

0.000 

0.000 

P32 

0.205 

0. 194 

0.000 

0.000 

0.000 

0.000 

0.000 

P53 

, O.4O6 

0.014 

0.000 

0.000 

0.000 

0.000 

0.000 

1^34 

0.177 

0. 133 

0.000 

0.000 

0.000 

0.000 


P40 

0.323 

0. 105 

0.811 

0.030 

o . ca ; 

0.000 


P /, -1 

0. 150 

0.559 

0.090 

0.702 

0. 170 

0 . coo 


P42 

G.O58 

0.227 

^ 0.089 

0.268 

0.830 

1,000 


P/3 

0.24 S 

O.3O8 

0.000 

o.oai 

0.000 

0.00'’ 

' 4 ... ■ 

^■14 

0 . 2 20 

0.001 : 

0.010 

0.000 

0.000 

0,000 

■ ... ■ 


Oontd 
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Contd . . . 


Parame- 

ters 

111 it ial 
V alue 

‘I'to.ge 

I 

Stage 

II- 

Stc ge 
I.lI 

Sta 

lY" 

Stage 

Y 

Stade- 

VI 


0.030 

0.279 

0.637 

0.577 

1.000 

1.000 

1.000 

^f2 

0.004 

0. 550 

0.362 

0.423 

o.oco 

0.000 

0 . 000 


0.045 

0 . 170 

0.000 

o.ooc 

0.000 

0.000 

0.000 

6 f 4 

0.871 

0 . 001 

0.000 

O.GOO 

o.oco 

0.000 



4.580 

15.015 

11.707 

14.38 

18.090 

19.890 

20.000 

^2 

4.390 

15.028 

14.206 

19.98 

19.420 

18.810 

16.290 


0. 196 

5.077 

7. 107 

15.19 

19.020 

16.540 

16 . 8 S 0 


1.610 

4.332 

14.803 

14.83 

5.410 

5.920 

_ 


- 

20.42 

10.23 

10.67 

20.18 

21.60 

21.69 

''^2 

- 

45.73 

42.96 

47.95 

51 . 14 

31.84 

31.42 



9.78 

14.99 

8.85 

16.45 

14.31 

14.54 


... 

3.41 

11.80 

7.77 

4.28 

J * i-/ 

- 

Recovery 

12.80 

69.08 

70.27 

70.72 

71.13 

71 . 0 Q 

70.95 

G-rade 

79.40 

55.00 

35.01 

35.00 

35.00 

35.00 

35.00 
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compared to other cells. It will be seen from coluiTin 8 of 
Table 8.4 that the recovery has fircther decreased somev.’hat 
but the final circuit is now considerably simpler. The last 
two structures are shown in fig. 3.55 stage 5 structure in 
broken and stage 6 in solid lines. ■ This example illustrate Ss 
given the necessary data, it should be possible to iiicorporate 
capital investment and operation cost also in a more comprehen- 
sive profit function for an objective decision on choice of the 
most efficient circuit. It is reiterated that the final struc- 
ture shovm here may not necessarily be optimal in the global 
sense. Nontheless, in view of the very large nuraber of possi- 
bilities, even the suboptimal or near optimal configuration is 
very difficult to synthesize based on only plant experience, 
intuitive feeling or even test results of a single flotation 
cell. 

THREE SPECIES ilTE BaTTEHISS Of fLOTAflON CELLS: 

Batteries of flotation colls, frequently used in industry, 
can also be synthesized and designed using the integral approacu. 
A battery consists of flotation cells arranged. Ei scriec so 
■that the tailings of one cell goes to tbo next cell in series. 
Vfeile the p-arallel concentrate streams’ from each cell combine 
together to form concentrate stream fro© the hatteiy, the • 
tailings of the last cell in. series constitute .the tailings 
stream. A typical batte|!y with four flotation cells is sbow-n 
in Eig. 8 . 6 . A lAth four ■batteries will be 



* 1 1**/ 1 ; ; 
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identical to the one given in lig. 8 . 3 . 


Eowevers it should 


he recognized that each block ( 1 to 4) in tnis figure now 
represents a battery sirtiilar to one shov/n in lig. 8.6 rather 
than individual flototion coll. If a single battery consists 
of n^ flotation ceils each of viiicb is a perfect mixtui'o, and 
if the mean residence time of each cell is equal to t, it can 
be shovm (See Appendix C) that the mass flow rate ox ; 3 -th 
species in tailings is 




1+Ka, t 


n . 


( 8 . 32 ) 


If a. . is redefined as 

a,, = (8.33) 

l+Ka j t^ 

it follows that the analysis of the circuit is identical to tlx. 
case discussed above. For illustration purpose a flotation 
circuit with four batteries and a feed consisting ol throe 
species with same characteristics as discussed in the preced,irig 
section, is considered. Each battery in, the circuit corjsisto 
of four flotation cells. In this case, instead of defining the 
upper bound for the total volume the volume of each coll 
was predefined (25 cu ft). This problem was analyzed for the 
same objective, as in the preceding section, for a mass flow 
rate of solid feed equal to 100 Ibs/oin. Table 8.5 gives 
stage wise progress of linicag© and design parameters and 
oTcjective function. The structure is shown In lig.'. 
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Table 8.5: Values oi Designed Parameters for Stage-Wise 
Optimiziatiou for flotation Cell Batteries 
^ Spacies Problem in Case of Objective A. 


Parameters 

Initial 

Value 

stage 

I' 

stage 

If 

stage 

III 

Sta go 
IV * 







°10 

0. 67 s 

0.098 

0.043 

0,781 

1.C00 

6^1 

0.000 

0.000 

0.000 

0.000 

0.000 

^12 

0. 159 

0.442 

0.393 

0.092 

0.000 

^13 

0. 156 

0.269 

0.300 

0.070 

0.000 

^ 1 4 

0.007 

0. 191 

0.258 

0,057 

0.000 

^20 

0. 660 

0.630 

0.719 

0,709 

1.000 

^21 

0. 140 

0.341 

0.281 

0.291 

0,000 

^22 

0.031 

0.000 

0.000 

0.000 

0.000 

“23 

0. 168 

0.029 

0.000 

0.000 

0.000 


0.000 

0.000 

0.000 

O.OCO 

o.coo 

^30 

0.468 

0.036 

0.037 

0.036 

0.020 

6^1 

0. 100 

0. 242 

0.333 

0.216 

0. 277 

“32 

0.000 

0.000 

0.000 

0.000 

0.000 

633 

0.209 

0.114 

C.OlO 

0.000 

0.000 

“34 

0.223 

0.608 

0.619 

0.74s 

0.702 

^40 

0.219 

0.059 

0.001 

0.000 

0.000 

^ ,1?, 

c . 203 

0.488 

0.741 

0. 34 1 

0. 141 

6 3j p 

C . 2^9 

0. 175 

0. 110 

0.370 

0.545 


0. 179 

0.097 

0.000 

0.000 

0.000 

» 4-1 

0. 190 

0, 180 

0 . 14 S 

0 . 2 B 9 

0.313 

^’10 

0.78-i- 

0.789 

0.98* 

1,000 

1.000 

1 

^12 

^3 

^14 

0.000 

0.000 

0.000 

0.000 

0.000 

0.212 

0.291 

0.016 

o.ot^o 

0.000 

0.001 

0.000 

0 .<X )0 

0.000 

0.000 

0.003 

0.000 

0-000 

0.000 

Contd. . . 

0.000 



Contd. . . 


Pa garnet ers 

Initial 

Value 

Stage 

I 

ote.ge 

II 

otage 

III 

Sta se 
IV'' 







i ^20 

0.765 

1.000 

1.000 

1.000 

1 . 000 

P21 

0.003 

0.000 

0.000 

0.000 

0.000 

P22 

0.008 

0.000 

0.000 

0.000 

0.00c 

P23 

0,224 

0.000 

0.000 

0,000 

0.000 

P24 

0.000 

0.000 

0.000 

0.000 

0.000 

P30 

0.718 

0.741 

0.638 

0.825 

1.000 

hi 

0.012 

0.007 

O.GOO 

0.000 

0.000 

P32 

0,000 

0.000 

0.000 

0.000 

0.000 

!^33 

. 0.079 

0.051 

0.000 

0.000 

0.000 

P34 

0. 191 

0 . 201 

0.312 

0. 175' 

0.000 

^40 

0.352 

1.000 

1.000 

1.000 

1.000 

P41 

0.005 

0 . 000 

0.000 

0.000 

0.000 

P /12 

0.010 

0.000 

0.000 

0.000 

0.000 

P43 

0.010 

0.000 

0,000 

0.000 

0.000 

P44 

0.023 

0.000 

0.000 

0.000 

0.000 

"f1 

0. 106 

0.031 

0.000 

0.000 

0.000 

2 

0. 181 

0. 154 

0.017 

0.000 

0.000 

0 

0.129 

C. 117 

0.129 

0.448 

0.495 

^f4 

0.534 

0.63S 

0.853 

0.552 

0.504 

."1 

■ , t’ 2 

2. 92 1 

3.613 

•;-.203 

9.681 

8. 5’:' 4 

4.274 

8. 190 

3. 137 

10.624 

12.732 

t.. 

4 . 3 f 0 

12.365 

12.660 

9.220 

9.54 

t, 

, ‘"T' ....'. 

2.69^- 

4.012 

3.000 

4.587 

4.376 

Recovery 

45.35 

70.57 

70.76 

71.23 

71.27 

Grade 

51.22 

55.00 

55.00 

' 35.00 

35-00 
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is iiidependsnt oi particulate environment in the suspension. ■ 
Jiq_a8,oion (8.2) for particles of size L can then he written as 


(L) [ -^ '1 
1+Ea(l) t(L) ^ 


(S.35) 


Now the design parameters to he optimised are tbo ccn-tants 

r "2 


C.,s Cq 5 ••• in an empirical expression relating t to size 1 


t(L) = td; 


'1 


Cq , • - • ) 


( 8 . 36 ) 


8 -'. 


Experimental data of Woodburn et sho?/s that the following 

relationship is perhaps sufficient for this purpose 

t(l) = c,, + cd + ♦•• (8.37) 


The main drawback in the integral approach, as indicated 
earlier, is in the optimisation step which he co .mo s more acute 
as the Parameter space increases. After this work v/as completed 
the author has come across a more efficient hypereonical random 
search method ?/hicli may prove to ho better for this kind of 
problem. 

Using the algorithm developed in Chapter 7, the apparent 
rate distribution can he determined and then discretisod in tv« 
or more groups of species. This may he incorporated in a genera- 
lized circuit in order to synthesize either the opttoal or a 
number of near optimal circuits with residence times as the 
design xjaraneters which play a central role in the kinetics oi 
flotation and separability of species* The optimal circuits 



were found bo be strongly dependent on the objective function 
which in industrial practice changes freq^uently vatn availabilit 
of ore 5 market price and deinand trend. Consequently 5 the 
procedure enumerated here provides a rational and objective 
method for design of new and adaptation of existing circuits. 
Even a marginal Improvement in the performance can make sijgni- 
ficant economical impact in vie’r of the very large tonnage 
treated in typical flotation circuits. 



CHAPTER 9 


SIlffiARY 


In tne xirst part of this investigation sorae quantita- 
tive aspects of flotation kinetics have been studied. These 
inclucie (l) a phenjnienological roodel for flot'j.tioii kinetics 
in distributed flotation rate constants (2) an algoritl-im 
for nuiiierical inversion of the Laplace transform for esti- 
mation of the apparent flotation rate distribution; (3) the 
effects of aeration rate ^ particle size and pulp density on 
the flotation rate distributions; (4) the effects of particle 
size and feed rate on the flotation rate distributions in a 
continuous cell. The second part deals with (5) the synthe- 
sis and design of a . flotation circuit. 


1 . A j?he>^enc3a3_gj^al Mo del for flo tation Ki_ne t icj 


Equations (4.18), (4.23) and (4.25) describe the 
kinetics of flotation under assumption of a rate distributed 
first order particle-bubble attachment sub-process. The uuasi- 
Bteady state conditions have been , stipulated, which are 
necessary to reduce these equations to a first order distii- 
buted apparent flotation rate kinetics, form, which has been 
postulated by many previous workers. She relationship between 
the apparent flotation rate constant and the i^rticle-bubblc 
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attacixcnent. detachme-nt rate constants is given in equation ( 4.5 
The apparen-c ilotation rate and the attachment rate distributi- 
ons arc; related to each other in equation ( 4 . 41 ). Wheii the 
detachment probability of the particle is high, such as tor 
coarse iparbicles, the first order distributed apparent flota- 
tion rate kinetics is valid for single sise particles only. 
Ho'vvever, v/hen A predominates, the relationship betv/een the 
apparent and atta-ichment rates is given in equation ( 4 . 43 ) and 
the apparent first order kinetics prevails even for a dis- 


persed particle size spectrum. I'or e given size, a distribu- 
tion of valuable mineral in attachment rate constant is 


introduced starting' from a tri variate density function in 
surface reactivity, particle size and -grade. It has been 
pointed out that the flotation kinetics can be formulated 
concurrently for distribution in total aass and valuable 
miner::;! using the apparent flotation rate constant as the 
cormuon attribute. These two distriUitions are related to each 
other through the cuiaulative grade as shown in equation (4.49). 


2 • Estimation of Apparent ElpJcrtipn Rate. JDist r.ibut ioii^ bj 
kuTiperi Cci 1__ Irrye r^iengjo f_ t h_e * • 


The problem of estimation of the distribution function 
of the apparent rate constant is equiv; lent to the inversiori of 
the Laplace transform from M^Ct) to lig{Ka)» It v«s found that 
the estimation of the rate constant distribution in cumulativo 
mode was easier and more reliable as compared to that in the 
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frequency xunction because the foimer is smoother and less 
prOiie xo influence oi errors. Moreover, the cuxiulative distri- 
bubiOii pxOVj.dcs a iiumber of constraints which not only constrict 
the region of numerical search but also prevent possible v.lld 
fluctuations in the computed distribution. Approximation of 


the integrcii in equation (5*12) by the lobatto quadrature for- 
mula automat icall^;" satisfies the constraints at the end points 
of the cumulative distribution. The cumulative distribution 
has been evaluated at discrete points by opt imizsttion of a 
least squares objective function with inequality constraints 
using univariant -search with Fibonacci search technique, and 
B. 0 senbrock ’ s method of optimization. The proposed algorithm 
has been tested for synthetic data generated from different 
distribution functions. It has been shown that by using this 
algorithm it is possible to extract back the distribution which 
is used to generate the data. Data accurate to two decimal 
places in range 0.01 to 1 has'been inverted satisfactorily by 
the proposed method. The success of the method does not seem 
to depend on the initial guess values in the optiiisation step. 
Off starting points also give accurate results but at the 
expense; of efficiency with respect to computation time. For a 
8-point Lobatto quadrature foimula, as small as seven menber of 
well spaced data point& lead to reasonable results. With 
appropriate modifiCc'.tion, data pertaining to feed with non- 
floatable component can also be Inverted satisfactorily. 
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3 . M and 

Rate D istr ibii'tions : 


PulolishGd flotation kinetic datg has boon eraployed for 


computa fcion of the flotation rate distribution by numerical 
inversion ox the Laplace transforiijj using the algoritlim developed 
in the preceding section. These distributions and corresponding 
apparent flotation rate constants, have been analysed as a 
function of aeration rate, particle size- and initial pulp 
density. Thus, it has been possible, for the first tjiae, to 
scrutinize the distribution as a v/hole for the effects of these 
variables. The principal conclusion is that the nature of the 
function R(I^) does not change with variation iJ3 aorc.tion 
rate and particle size. Only the argument of the function, 

Ea,is altered or rescaled as in equation (4.43) and (6.14), 
and the scaling parameters may be determined from the appropriate 
data. The scaling factor increases rapidly with the rate of 
aeration and exhibits a characteristic maximum for inteimediata 
particle size.. The evaluated distribxit ions are independent of 
initial pulp density when dilute, but show marked divergence 
for concentrated pulp. In this manlier it has been possible to 
develop useful correlation for the structure of the apparent 
flotation rate constant Ka and to check the valluity of the 
conclusions derived from model in part (l)* 
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4 . . of Particle Size and Feed Rat e f he I 'lot ctic n 

Rate Distributi on in a Continuou s Cell : 

The approach used lor analysis of flotation kinetic data 
for seuii-batch flotation operation has been eztendad to the 
analysis of the flotation kinetic data for continuous steady 
state operation. The cumulative rate constant distribution, 
RjiCKa,!), in equation (7.5) has been evaluated by I'lSing a method 
similar to numerical inversion of the Laplace transform. The 
integral term in equation (7* 5)? after appropria I'e transformations 
has been discretized by a Lobatto quadrature fomiula and the dis- 
tribution at discrete points is evaluated by optimizing a least 
squares objective function by univariant, Rosenbrock and Random 
search methods of optimization. The method is again not restric- 
ted to any a priori assumption concerning the mathematical form 
■of the distribution function. Using tha tabulated data of 
Woodburn et al. the effects of particle size and feed rate on 
rate distribution have been examined. It is observed tliat these 
distributions are independent of the slurry ilori rates and the 
effect of piarticle size on the distribution was sinmlrr to that 
in the semi— batch process. The estimabed rtste disi-ribution has 
been employed to predict quife accurately the Icinetics of a 
cleaner cell ?;hose feed comprises of the rougher concentiote. 
Similar results have been obtained for gangue component. An 
algorithm, given in equation • (7. 5l) , for analyzing' the kinetic 
data in the impulse case has 'also been outlined. 
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5 . 

Tbe optimal otructure and design parameters are extracted 

simultaneously from a generalized circuit. This technique has 

been applied to solve the linkage (including split streams) 

problem and to determine thp n-n+imoi 

lixxie une optimal mean residence time of 

particles m each cell, first, a flotation circuit with two 

cells and feed comprising of tw species of Iraovm overall rate 

constants is synthesized and designed for three objective 

functions: (i) maximize recovery such that concentrate grade 

> 35 percent; (2) maximize recovery such: that concentrate 

grade > 75 percent; and (3) maximize profit, which is defined 


as tne difference between the prices of concentrate and ore. 
Lodified complex method and Random search techniques v^re used 
for optimizing the objective function but Random search seems 
to be more reliable. The results obtained by using this approach 
have been compared with the real optimum circuit obtained by 
using the trial and error technique in v/hich all possible 
meaningful configurations are enumerated (fig. 8. 1) and each 
configuration is then examined separately. The optimum configu- 
rations for above three objective functions are identical to 
those shown in figs. 8. l(vi) ,8. l(v) and 8. l(vi), respc'ctively. 
This approach has been extended to a flotation circuit with 
four cells and for batteries of cells for a feed -jeoBorising 
of Valuable, gangue and middlings. The optimal configuration 



for these two 'cases for objective (i) are shov/n in Figs. 8.4 
and 8. i, respectively. Simplification in the configurations 
which result rn somewhat smaJlar values of tho objective functic 


have been obtained by using engineering guesses, for four cells, 
and resulting siraplified configurations are shown in Fig. 8.5. 
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SITGOSSTIONS ROR FUTURE YORK 

ihe continuaGion of the present work may be carried out lo 
the following directions. 

1. Ihe importance of the study of the effect oi process 
variables on flotation kinetics has been duly emphasised in 
the previous chapters. In Chapter 6, the effects of aeration 
rate, particle size and initial pulp densities on the flotation 
rate distributions have been studied. The effects of seme of the 
other important variables such as chemical reagents, e.g. colle- 
ctor and fro ther etc. , agitation (which is related to the 
turbulence in the cell) and froth level etc. may be studied. 

This will give further insight into the structure of the flota- 
tion rate constant ?/hich may be important from the view point 
of the control of the prooes.s and for eventual derivation of the 
flotation rate constant from first principles. 

2. As discussed in Chapter 6, the assumed structure of 
apparent rate constant Ra in equation (6« 14) is sonievAiat similar 
to the model proposed by Woodbum et al and given by onuation 
(6.25). Yftien only one variable is changed, either aeiation 
rate or particle size, both models require sel.v.— similar 
distribution in scaled apparent flotation rates. However, in 
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eq.uatio]a (6.25) parameter 0, which is a functioii of particle 
size and a hydrodynamic parameter, does not permit separation 
of par bide size effect as in the model proposed in the present 
study. To discriminate between these two structures of the 
apparent flotation rate constant, experiment s ? at the laboratory 
scale, with particles of different sizes and at different 
aeration rates for each of these particle • sizes ms-j be carried 
out. The data thus obtained can be analyzed as follows: 

I'irst a plot of R(K) versus ICa/Ka(0«5) for different 
particle sizes is obtained, which is expected to result in numbe. 
of curves of the type shovwi in Eig.6.2, each curve correspondin ': 
to a particular particle size. If these self similar forms of 
distribution collapse on to a single curve, similar to the 
curves shown in Fig. 6.8 to Fig. 6.11, the validity of equation 
(6.14) is established. 

3. The problem of synthesis and design of optimal - subopt imal 
mutlticell flotation circuits, discussed in Ghap'cer 8, is rather 
simple. A more comprehensive pro fir function can be formulated 
incorporating the capital investment, oper'^.tion cost and price 
of the products etc. The analysis can be extended, to particles 
of various sizes and densities as Bug.gested in Chapter 8. The 
random search technique oi optimization is perhaps not the best 
technique of optimization. Faster and more reliable technique- 
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of optimization ?rould be more desirable. Directional random 
search techniques seem to be more appropriate for this type 
of optimization problems involving large number of constrained 
variables. One such technique is the hyperconical roridom search. 

4. Recently; a number of techniques of synthesis of chemical 
reactors have been proposed which take core of some of the draw 
backs of direct search technique of synthesis, used in the pre- 
sent study, A new method of synthesis proposed by Ichikav/a 
and Pan and a multilevel technique proposed by Osakado and Pan 
seem to be encouraging and feasibility of these techniques for 
synthesis and design of 0 multicell flotation circuit may be 
investigated. 

5. In formulating the kinetic model in Chapter 8, it is 
assumed that the mean residence time 01 water and solid partic ■ 
in the. cell is same, A -more sophisticated mode. L, incorporati.ii. 
different values of mean residence time for different species, 
can be formulated and this model then .can be used for synthesis 
and design problem, discussed in Chapter 8 . 
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JkJiGOillTa^S POR NIMERICili BITERSIOil 01 


THE LiPLACE TRAIsSi'OR.. 


AS pointed out in the text, if 

Mg(t) = / Mg(lCa) exp [-ICa tj dKa (A-l) 

o 

Ivl (t) is tooTOi as the Laplace transfom of function 1C (Ka) , 

hj 

Mg(Ka) is the inverse of the Laplace transforiTi. For equa- 
■tion (A-1) to he valid Mg(Ka) must he defined for all positive 
■values of and in our problem it can not take negative 
values. 

One of the commonly used methods for inversion of the 
Laplace transform consists of approximating M^Cla) by Jacobi 
polynorflial"^'^^ after carrying out suitable transformation to 
change the limits of the integral, in equation (A-1), to -1 to 
1 (see Chapter 5). Thus Mg(Ka) , which after transiorrairig ICa 

to X is denoted by M’(x) , is given by 

A P °° ( a , p ) . 

K*(x) = (1-x) (1^) I (x) (A-2) 

n=0 ^ 

wbere 4. j a'? and p are the unknovm parameters, C^(n =0,1, 
ar-e coefficients in series and P^ is the Jacobi ix)lynoCTial. 
Coefficients b^ can be analytically dateriained once the series 
is truncated using some suitable criterion and values of 
parameters A , V , a’ and.p* are' suitably chosen, mfferent 

values of these parameters slid ttchis_Hues for evaluating the 
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C 0 © xf iC i_ 011 b S t) llQ-VS X6SDli~Prl 1 n *t7 -VO n ^ n - » « 

■ Cl in Various algorit«niiiS| doeuiBerited 

in literature. Srdelyi''^*= ’ has • ■ • - x- \ 

uas given expncii- ior;:i 3 for 

but in the method proposed by Papoulis'* and Miller and 

” "^n embodced. in a S 3 j'stein of linear equations, whose 
coerficit^nt matrmc is triangular. Iiiaerical examples based on 
these approaches have been given by Pelix, Sajalsli and 
Kuntzman , and Nugeyre 1 Erdelyi obtained a general inver- 
sion formula analogous to equation (A-2) by use of an orthogonal 
set of functions over the interval (0,“). 


Bellman, Kalaba and Lockett-^'^^' have discussed the 
numerical methods of inversion in which integral in equa- 
tion (A~i) can be approximated by either Laguerre, Legendre 
or shifted Legendre quadrature formula after carrying out 
proper transformation of the variables. Phis leads to an 
equation of the following form 

Mg(t) = L f(x.,t) (A-3) 

i=i ' “ , 

where the quantities l'„(x.)? i ~ 1,2,..., n are unknoTOi 
function of x,. and f(xu ,t) is .lome knov/n function of and t, 
the form of which will depend on the quadrature formula used 
for discretizinv the Integral in equation (l-l), and wp, 
as defined earlier, are the i~th abscissa and W 2 ight, respective! 
of the quadrature formula. Now allowing the ^ variable t to 
assume n different values, i - tjS,..., n, one can obtain a 



system of x. equations in n unknc^s. Recentlv Pi... 
has^ modified Beliman^e method to avoid the loss of 
cant figures inherent in earlier methods. He has she 
his,, as well as, Bellman- e method are theoretically 
to a more generalized method proposed by lanczos^^'b 
the methods are suitable for noise corrupted d^'t,-'. 


111 


oirr.ifi- 
■’VO' that 
equivalent 
Hone of 



B 


O^TIialZATIOB TEC Hi" 


Ghe j-ast; "fewo der"'ii<^>-. - 

^ec,:.aafa, a large iiueIg- of algorithms 

ha'^G 1 e 011 proposed for rr i-t-rt. 

^r ma..ira^sxiig or miriimisin 

fu.K'.^'biOii s with or ivitlioiit pTni'^i -i-i-tr - 

-out eo^slity an::, inequality coastraints, 

using numerical techninuep nniT. ps, 

e-. Onlj ^.nosQ optimzction tochni- 

ques which have been emni mr^a ■ i 

oen employed m rhe are sent study are dis- 
cussed in this Appendix. These are univariant , Fibonac .'i , 

Rosenbrock, modified complex and random search met’uods of 
optimization. ' 


UlilVARIillT SEARCH: 

ihis method, also known as one-at-a-time , or sectioning 
meohod, has been basically developed for unconstrained 
objective functions. However , due to the special feature of 
our problem in Chapter 5 this method v/os found to be quite 
suitable for optimization of the objective function in equation 
(5.2l). lor the sake of clarity the objective function can be 
rewritten as, 

XI = E( X 2 ) , ... .| ) , ( B— « 1 ) 

subject to constraint 
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i/vhere no?/ corresponds to ) in e:,paation (3.15). iiie 
univarxant search used in this work is siEilar to I'rie juan ml 
Savage method and is as follows: 

Choose ^;n°l^ as initial lase point 

such that tine constraints in enuation (B-2) are co.ti’sfied. 


_Step__ 2 : 


The Cb.jective function is no?/ minimized wite rtspoct 


to one variable onljr, say, Xp, keeping all other variables 
constant. Fibonacci search technique (to be discussed in the 
folloudng section) in ranee 0 to is usei for this search. 

/.N ■ ■ ~ 5 

If X 


( 1 ) 

2 is the optimal point, the new base point new becomes 


/y(l) j(0) 

\ ^2 ^ ^ 


llext, the objective function is mini- 
mized with respect to X^ keeping all other variables consto^nt. 
The search is restricted in range Xg"^^ to rather than 


y(0) 

^4 


and X';;^^. This procedure is continued until the objective func- 


tion is optimized with relspect to all the variables once, to give 
a new vector (Xg'^^, ’**’ completes tlie first 

stage of computations. 

St ep 3- The v/hole cycle ics repeated unril there is no signifi- 
cant change in the objective functiOxi, or until vhc cnongc is 
within a preselected value. 

A flow chart used for computer programming is given in 


Fig. (B-l). 




w£W^m 




Read parameters 

in ' 


Choose .initial, 

such that constraints in- 

Eq. CB-2) are satisfiaa 


Evalwats functional 










Cilitcift''' for' ■’ 
^stoppinf critorioft^ 


Fig. B-1 Flow chart J 
t' ' Qt'O 









PIBOIAGGI SEARCH: 


l!37s 1 5- 


This technique is used for optimizing objective fuze 
tions with a Single variable. The steps are as follows 

points c.,nd (x^ > x^) are placed at a distance l^,froiTi 
each end of the original search interval ly The distance 
1^ is defined as 


(.InnS £ s - 

^ i' ” tv ) 

n n ‘ 


(3-3) 


where b is arbitrarily predefined aaall number (say, 0.0C5) 


which defines the distance between the last two trial points, 


'E ' defines 
n 


the 


h“th number in the 


sequence of Fibonacci 


nmnbers and n is the total number of trials finally required. 
The sequence of Fibonacci numbers is defined as follows; 




= F + F 

3-1 , 1-2 



(3-4) 

n' 


The values of the oblective function at the tw iriul prints 
are compared and one of the three sub intervals contatjiin,3 the: 
worse rrial point is relucted such that the rejected trial point 
forms one of the two end points of the remaining oubinterval 
and the other trial point lies in this siibinterval. 
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The new reduced interval of search 
after first two trials, is 




rsfiE ^ L-.l)fL.i; 

1 V 1 



Nov/, tv/0 trials in this region are placed at a 
from each end of the interval such tnat 


i *1 r-i -L 




v/here is defined hy following equation 

■‘' d +2 " ^2 ~ ^* 0+1 . ( B - 7 ) 

It can be shov/n that one of the tvvo iiev/ trial points virin 

exactly coincide with the trial point remaining from the 

previous search cycle. Thus, in second cycle, only Qri,- 

■' new 

trial point is necessary to determine the third scorch 

interval. The new trial point is placed such that tho v 

nev; 

point and the old point are syixane trie: xly placed iii th^ 
subinterval of search. This holcte for subsequent dVelo- ^^30 
Thus, only one trial point is to be added in each Again 

the Values of objective function at these tv;o i-oiritp 
compared and the subinterval v;ith the worse point 13 

Step 3; The search is carried out in this manner, 
trial point symmetrically with respset to the trial 
remaining from the previous cycle in the subinterval 
until all but one trial point' are placed. 



last trial point is placeu, in ths siibii.terval 
accoreiing to dichotomous search plan, that is, riev; point is 
placed syEiiDtitrically with respect to old point nt ai^ aioitroi'i!.'' 
small discance £' from it. Again the sub interval with v.orco 

point is rejected end the subiriterval L / reirainin;: after 

. ■ n ' 

this trial, becomes the final region of uncertainty. 


In Fibonacci search, the number of trial points ia to 
be predefined, larger the number of trials, ae^iller is the 
region of uncertainty. In order to ;-ucelerate the iinivwriant 
search, in Laplace inversion problem, in first stage only three 
trials were prescribed. After every additional stage the 
number of trials was increased by one to maximur' 11 trials 
in 10 stages. This results in a maximum possible final regio’'' 
of uncertainty of only 0.5 percent. This procedure, which is 
computationally more efficient, begins v/ith coarse adjustment 
initially with progressively fine 'tuning' of the vector 

^ at the final stages. Flow chart for fibonacci 

search is given in Fig. (B-2). 


ROSmiBROCK I..STnOJ)i 

This method incorporates a .direct search otrategy ii-hich 
gives acceleration in both direction and distance. The compu- 
tational algoritto is illustrated for miiiimizing the folio wi. 4 , 
objective function 


:5ji 

























Altnough, basically the original nethod v,aB Inclc -c:i for air-.: 

"tiQ-iiiGd. 'v3X'i.QDl0s oijIyj j_T Ip 

the constraints also. 


beeri rnodified .n crate 

Ihe tiottiflsd mctbod la siBwarlsel as 


Select a set of s sirtaelly oi-:-a.o.::cr.;,l dircctioa 
coinciding initially i^ith the coordit 
not necessary^, let this set be 
:( 0 ) 


late 1 i r «'' c o i o ri n ( n 1 1 h a \ 3 ; h 


= (1 


p ■ 


c(o) 

'•52 


0 ) 

= (0,15...0) 


Ct-9) 


= (0,0,...1) 


Let ®2’*‘*’®n the set of step sizes and 
be che starting point vdiich satisfies all constra iirts. 


.§is.P.._2: A new trial point (X 


(C) . e v(0) y(0)\ V. 


obtaiiieci anu the ■value of the objective function at tliis point 
is compared with the functionaj value at the storting base 
point. If bhe trial is a success, i.e., the functional value 
at the new trial point io araller than that at the.’ bane noiirt, 
the step size is increased by a factor a (o > 1). However, in 
the case of failure _ the ' step sise is reduced by a factor 
p(0 < p < 1) and the oM base point is retained. This procedure 


is repeated until one success and one failure io recorded to 


(c) 
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175 

direction . Rosen'brock recornitienaed 
If the f-uTictional values at the trial ' 
are sairiej the trial is considered as a 
that there v/ill be atieast one success 
because if there is ic success, the re 
will ultimately bring tne trial point ' 


G = 3 and f, = 0.5. 
point ana the base poin 
sue ce 3 3 . Th i s cn s’’. -.re s 
and cue failure 
due tier, of step sene 
to orijin-'l base point. 


4 ..^ 

y 


dhe search is now carried cut alone 




..( 0 ) 
* 't'n 


directions sequentially using toe Sana criterion of one success 
and one failure in each direction. Tos net of trials iraie with 
one set of directions and the subsequent change of these dii-ec- 
tions is called a stage. 


Step 3^ Once a stage of computations is over, a nev/ set of 


orthogonal directions is generated by using Grain Sclisidt ortho- 
gonalization ruetbod which is as follows. 


Let sdg, . . . be the algecraic sian of all the success 

( id f ■) f id 

ful steps in direction. » da"' ’ • " ’ ^^n ^ ” 

independent vector be defined as 



- i r ( .1 ) X -1 ^ + + 

d r; ^ 3 ) 
n n 

Ur) 

jfC i1 ) j. j. 

A f ( .1 i 

Si ’ ! h 



d r ^3^ 



n ^ 

Orthogonal 

unit vectors 

-(3+1) 
• *-** ^ n 


obtained in the follovdng tBaaner 


(B-IO) 


are now 



223 


c(d+i) 


f ( D + 1) 

■'2 






u. 


•11,, 


-(■3 

31 



,- ( ui“-'i ) 
1 1 


% 

t ( D +1 ' 
’n 


= u 


M 


ii— J 

E 

i=1 




■:. (;i+i ) c (i+i) 

■^1 


i a. 


(B-11) 




where ...j refer to a set of ortbogorial vectors 

which is norraalized to give a set of orthonomal vectors 


t(3+l) c(d+1) c C3'^1) 

'T^ , ;2 s --'j '-Jri 

Stj_,p_4; The process of generatirig the set of ortnonomal 
directions and then searching along t’;. sc directions is continued 
until a predefined convergence criterion is satisfied. 

If there are seme constraints to be satisfied, nt each 
trial point variables are tested for constraints. If a particular 
trial point violates the constraint, the trial is conbidered 
as a failure. A flow chart describing the computer program 
used in the present study is given In Fig. (B-3). 

























226 


COkJ^bEX MEiTHOl): 


-he coxnplsz method of Eox”''^® 
iQio'Tn Simplex method in two respec-ce 
not regular and. secondlj, it can be 
and/or explicitljr constrained optiiai 


d j.. f f 3 1* 'f 0'^ " 

. firstly, t 
used to solve 
sat ion problc 


the v/oll 
he complex is 
implicitly 
::.s. iTlxe method 


is more useful if the variables are bounded, i.e., 

^iL - - '^'i 21^- and refer to Lov/er and upper 

limits of Xi, respectively. Even if 'che bounds of variables 
are not defined in the problem, some guess can always be made. 


The algorithm, for complex method for an optimisation problem 


?/ith n Variables is as follov/s: 


S tep 1 ; Generate a complex with k vertices where k _> ri+i. box 
suggested that for n upto 5? k = 2n but for ri > hi k can be less 
than 2n. To generate a complex, first a base point which satis- 
fies all the constraints is obtained. The remaining of the 
k-1 vertices are obtained by using the iollo?.'iii;j epueitioh, 



=^11. •- =^.13 - *11.^ " 

- b' 1."^ i 


1 ,2, . . . ,n 

1 j 2 j . . . 5 .v» 


(B-12) 


?/h6r6 X- . rsfers to i— dh coordr/fioto of vert ox oikI .rjj 

is c random number usually betweeii 0 to 1, but this rouge 
can be modified depending on the type of constmints on the 


variables. 



227 


of fbc Vertices, so determined, violates any 
of oue cons ..raint s, move half way in towards the ccntj^oid 
of the already chosen vertices, fhe base point wnich satlsfie 
all constraints ensures ttie generatj.cn of the coriulex wfioso 
vertices satisfy all constraints. 


Ste£j;, Comparing the fimcticnal value, 


lli . -4 ,1 “» i 5 C 1 • • # I 


the vertices find the worst point and its position i 


L!i si.'ace, 


A new complex is formed hy replacing tne worst ■■>cirit by a 
point a times as far from the centre id of the j-'emaining points 
as the reflection of the v/orst point in the centroid. The coorc 
nate of this new point, X-. are 


. k 

f £ 

(k-1) 3=1 




j, — • 21 


(B-13) 


where X.-o is the i-th coordinate of tne worst potrit to be 
rejected. 

Steji 4: If the new point is still a ’vorst point or violates 
some constraint, a new trial point is ohtaine • by noviiw' half 

Vv’ay. toY/ards the cGiitroad... 

step This cycle is repeated until s precpecified stopping 

criterion is satisfied. 

To improve the reliability of the method £3oae n»diliC8"- 

, 162 

tions have been proposed by Ifeeda and IchitoTO f snd Gu.n 
and are discussed belows ' 
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1 . 


While evaluat: 


ing the ceritroid cf the renaL'ir;:: eohit.: 
the functional values at vertices of the compleJ. should also 

be given wei^^htags. Thus, the modified coordinates of the 
centroid are 


'•iC 


S X.. 

3:.1 3 ID 




k 

A 


> i = 1,2, 


(3-14) 


where m is a positive number, an: 


,£%. 3 . 


3 S . 

S - E . 

max min 


- Ej^r, 

E 


for minimisation problem 


for maximizat ion problem 


(B-15) 


max min 

where E, „ and E^. are respectively, maximum and minimum 

' mill , 

functional values. Eow the nev/ trial point become:. 


k 1' 

£ (-’d^S.)® X.. 

,1=1 


^i'i'T ~ ^ 1+0;) 


k . „ 1 

1=1 


- a X.p; i = 1.2 


jt,. 


,n (B~l6) 


2, Above modification should be used as long as 

W/W >:°-5 

When this inequality does not hold, original complex method 

f, , 

should be applied by taking' ns » <3» 



3. In the case of violation of inequality constraints, insieo ^ 
of moving half way towards the centroid of the i-oimeinin, points 
a more general mj of moving is defined as 




a'Xi-, + (i-a ) x.j^. i ^ i,2,...,n (l-lR) 


v/neri,. u j thi^ rts-ileccion factor, isj ODts.ined iii such a v/a’/ th-.-rt 


the functional value is optimum with re, 3 pect to X. 

in 

4.' The method may fail if all points lying on the line joining 
the worst point and the centroid of the remaining points are 
vvorse theai the original point because eventually the projected 
point will coincide with the original point. To avoid this 
situation 'duin suggested that if a falls below a certain value 
(Say, < 10“^) v/ithout obtaining a better functional value this 


point should be replaced by its original unproj acted position 
and the second worst point be rejected. 

5. In the case of non-convex curves, centroid itself may violat,- 


the constraints. If this happens at some stage, all but the best 
point of the complex are discarded and a nev; cojaplex is generated 


according to 


X. 


3-3 


■^ib 


t >1. (X^^, 


:X, 


ib^ 

3 


^ 1 2 j , n 

1,2,... ,k-1 


(P-19) 


where and refer to i-th coordinates of the unfeasible 
centroid and the best vertex in- the. old complex. 



V’ -i- 
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G-eriera.t_i o Ti 


ial^a s,e„ : 


bLoug »-5 8 fsi/v iii©vliod.s axs aVcilal)!© ir; lil'arp'turfi to 
evaluate the initial base point, it v/as foimd that the notijro 
o± constraints in the present study -vis giacb that etiia-niori (?.~li 
itself could be used to gen-rsie the base point. A set of n 
randoiii numbers was generated and initial base point ttis obtained 
by using this equation. If this point violated some of the 


cons fcra3.nt s 5 a new set of random numbers was goneroted and the 
procedure was repeated until a point, which satisfied all the 
constraints, was obtained. Although this method of obt-jining 
the base point is not the mo st efficient, it is conceptually 
very simple and easy to apply. 

A flow chart of the complex method used in tlic present 
study is given in fig. (B-4-). 

RANIX)M SEARCH: 

One of' the main drav/backs of Rosanorock and complex 
methods is t^iat the solution obtained may not be the global 
optimum. Sven large number of initial guesses donot guarantee 
the global solution. Ibis problem becomes even more acute 
with increasing number of variables. However, this oroblera 
is taken care of to a great extent by using a random search 
technique. Basically j the method consists of specifying 
each of the variables, to be investigated, in a random iitaiu'jer 
within the region of interest, Itech variable i« assumed to be 
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independent of others. Some n-p 

other attractive fectnres of 

this technique are: (i) ease iv, 

-^n computer proara!m'.:irr' : (2) free for 

from mathematical restriction p 

on functional form, and (3) corver-. 

gence in noisy environment. Detpii^ - ^ ^ , 

xtciiio ui tnese features and 

various random search schemes a,-e extensively locumcnloa in 

literature^^^’ T'bp . 

--- citeiis 3.nvolveo. in tne algorithiE 

used in the present study ore listed belo,.-''®. ihe proble- 
considered here is same as discussed above, reniely 
Minimize, 

li = Jl/(iL , Y \ 

' 1 ’ 2’’ * * 

Subject to constraint 


( 1 - 20 ) 


X-T < X. < X-^ 
iL - 1 "• 


i -- 1,2 


g £0 I 


'■11 


(B-21) 


S:Lej?_l: Denote the initial base point Avhich satisfies ail 

constraints by Xj^jOld ~ 111,...,!!) and the functional value 
at this point by 


Step 2; 


G-enerate a sot of n random numbers, (l = 1,2,...,n) 
by using some suitable subroutine and then evaluote a new trial 
point by using following ©QUation. 


^i,new 


^i,old i = 


Vvhere 


AX- 

1 




^ 2si/2 


w'( ^ r.) 

i=sl 


(Iw2?) 


(I'-23) 


where w is 3 weigbt^^ fftotor. 
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Step g: Test the lew trial T)nin+ -f-.- 4 . • ^ 

foi- constraints. If al.l eons-- 

traints are satisfied, fu-o+irr s’ - . 

.-’J-Lca, iU::etiori corresponrli:::, to nev/ 

trial point, is obtained and is compared witr 1 ^^,. if 
®new < ^old’ ca-iginal ard X, . are discarded and S. 

.XjWj«'JL II 0W 

and are stored end, respectively, become E , , .-.^nd 

OKi 

^ijOld* point is obtained by using equation (B-22) 


usin 


■g the same value of as; in the previous stage, 

,±: The whole procedure is repeated, if B > E . 

nev.' old.' 


oi* 


any of the constraints, other than those defined by (B~ 21 ), 
io viola oed. The old trial point is retained and the new trial 
point is discarded ana a new trial point is obtaiiied by using 
a new set of random numbers and equations (B-22) and (B-23). 
However, if constraints given by equation (B- 2 l) are violated, 
is replaced by corresponding limiting value, i.e., if 
Xf > X^, it is replaced by X^ and if < X^, it is replaced 



■Step 5: Ihe search is cDUtinued until a pre specified 

convergence criterion is satisfied. 

Detailed computational procedure is given in flow chart 
in Fig. (B-3). The number Ig in the flow diagram refers to the 
maxiiaum number of successful trials allowed with same /^X^. 

This is used to prevent' situation auoh as only a very sivall 
improvement is made on E in each successful trial •’ I.| is the 
number of successive unsuccessfhl trials after wliiob computation 

















ck ^ 
•&qml 
:Strai, 
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is stopped. Tbe assumption here is that successive 

unsuccessful trials mean that S has attained the optimum value. 

The number defines the total number of trials that are allo’.ved 
3 

whether they are successful or not. 



ippeithix c 


MASS BAIMGE EQUATIOB' H CASE 


OE A BATTERI OF FLOTATIOH 


CEILS 


Consider a bsttery of n iderMical flotatiou cells 
arranged in series in such a my tnai tailingg fron one 
cell Toecomes feed for the next cell in series whereas 


concentrate from each cell is collected in parallel. A 


batterjA of four cells is shov/n in Pig. 8.6. Consider cell 1 
in the br:ttery using the same notations as in Ghciptcr 8, the 
mass flow rate of o-th species in tailings and concentrate 
streams, respectively, are 

Mm.CKa) = Kp(la) [ — ^:; ] (C-1) 

^ 1+Ka t 

and 

M^.(i:a) = lilp(Ka) [ 1 (C~2) 

^ 1+1^ t 

It should be noted here that subscript j used in Chapter 8 
to denote 3 ~th species has been removed for stoplicity in 
notation but 7 /ith the understanding tint these arc appiicabl<? 
for any particular species, say j-th. Bov; for cell 2 in the 
battery becoiaes the mass flow rcte in feed, therefore, 
tailings and concentrate flow rates for cell 2 are, 
respectively, 
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and 


(iCa) l 

1-fKa t ' 


lip^ClCa) = lin, (1(a) 

1+ife t 


(G-r) 


(C-4) 


Substitution of equations (C-1) in (C-5) and (C»a) results in 

2 


M,_(Ka) = M^(Ka) [ — i— 1 

1+ICa t 


and 


M^gdva) = M^(Ka) [ 


Ea t 


TJ-- - L el L— -G] 


(C-.5) 


(C-6) 


This argument can be extended to n -th cell to give 


I%_. (Ka) = Mj^CKa) [ ] 


n 


and 


M (Ka) = M,,(Ka) [ — 1-~ ] 
g : 1+Ka t 1+Ka t 


n „-1 

C. ' 


(C~7) 


( C~8) 


It should be reccgnized that tailings flow rat'^ for the battery 
aa a wtiole (Mrp(Ea)) is same as t*iat ior n^.-tli cell, whereas 
the concentrate flow rate (if^(Ka)) is the iTumm-itio.j of concen- 
trate flow rate of individual cells in battery. Thus, 


.MgiCha) = (Kc’' 


. H,(Ka) 


IJ, 


l+Ka t 


(C-9) 
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and 


M.(Ka) 


M 


i=1 


Ci 


11 


i=1 ’ " 14-Ka t 


1+Kci t 


i-1 




II . 

i... ) .j 


1+Ka t 


Thus, if ?7e redefine a^- in equation (8.4) as 


n. 


a 


( — ) 

1+Ka t 


(C-10) 


(C-11) 


substitution of equation (C-11) in (C-9) and (C-IO) leads to 


M^(Ka) = a 


and 


Mp(Ka) 5= M-r,!(l'^) (l -r) 


(C-12) 


(C-13) 


Thus, the basic forrr: of equation (C-12) and (C-U) is same as 
of equations for a cinple cell (oqur'clons (B.h) and (8.6))* 




5^3- 1>- M'tsM 



